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STATEMENT  OF  THE  PROBLEM.  - - — 

Photoluminescent  porous  silicon  samples  were  synthesized  and  their  physical  properties  studied  and 
correlated  with  their  electrical  and  optical  properties.  Samples  were  prepared  by  electrochemical  anodization. 
Characterization  was  carried  out  using  SEM,  TEM,  AFM,  Auger,  Raman,  time  resolved  spectroscopy  and 
standard  optical  and  electrical  characterization  techniques.  Electronic  surface-states  densities  were  determined 
using  a  novel  electrolyte-semiconductor  blocking  array.  Phototransport  properties  were  measured  using  the 
photocarrier  grating  technique. 

The  electrode-PSi  formation  was  investigated,  for  electroluminescence  applications,  by  promoting 
electropolymerization  in  the  inner  walls  of  the  pores.  Applications  of  porous  silicon  to  the  silicon-based  UV 
detector  technology  were  explored. 

Nanocrystalline  silicon/Si02  composites  were  synthesized  and  their  optical  properties  studied  and  compared 
with  PSi.  The  study  of  the  effects  of  the  nanostructure  on  the  rear  earths-doped  silicon  emissions  was  initiated 
with  Erbium  and  Europium. 


SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS. 

•  Luminescent  porous  silicon  films  synthesized  from  crystalline  silicon  wafers  with  a  large  variety  of 
doping  levels  both,  n-  and  p-type,  were  obtained.  The  samples  had  different  pore  morphologies  ranging 
from  porous  sizes  of  the  order  of  microns  in  width  and  columnar  orientation  to  nanosize  porous  with 
random  orientation.  Self-supporting  flat  nanoporous  films  with  optical  quality  were  prepared  and  studied. 

•  Combined  studies  of  the  photoluminescent  spectrum  and  the  electronic  surface-state  density  and  effective 
surface  area  if  the  porous  silicon  surface  prepared  from  Si  p-type  wafers  shown  that  intensity  and  spectral 
position  of  the  photoluminescence  spectrum  maximum  changes  as  a  function  of  anodizing  time  attaining  a 
maximum  after  ~4min  anodization.  The  effective  surface  area  and  the  density  of  surface  states,  on  the 
other  hand  increases  monotonously,  the  surface  state  density  attaining  a  value  of -101 3cm'2  suggesting 
that  the  surface  states  near  the  valence  band  are  not  involved  in  the  luminescence  process. 

•  For  samples  for  which  surface  photovoltage  spectroscopy  revealed  a  band  gap  of  the  order  of  2.0eV,  the 
bandtails  widths  were  of  the  order  of  0.3eV.  The  majority  carriers  in  strongly  luminescent  films  were 
holes,  while  in  the  weakly  luminescent  films  they  were  electrons  even  though  the  starting  material  in  both 
cases  was  p-type. 

•  In  macroporous  Si  we  measured  an  electron  diffusion  length  of  2pm  and  a  diffusion  length  for  holes 
(minority  carriers)  of  0. 1  pm. 

•  Decay  time  analysis  of  the  photoluminesce  of  PSi  brought  to  the  development  of  a  new  approach  to 
extract  the  decay  time  distribution.  We  could  obtain  the  distribution  based  on  the  fact  that  the  non¬ 
exponential  luminescence  decay  behavior  of  the  sample  is  multi-exponential. 

•  Electron  Loss  Spectroscopy  analysis  showed  that  the  C02  involved  in  the  critical  point  drying  method 
contaminates  the  pores’  surface  with  carbon. 
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•  The  study  of  the  growth  process  of  conducting  polypyrrol  by  electropolymerization  inside  the  pores  of 
PSi  showed  coverage  of  the  inner  surfaces  according  to  the  curvature  and  consistent  with  the  electric 
fields  distribution.  PSi  samples  with  reduced  size  of  the  pores  showed  limited  filling  and  coverage. 

•  Optical  absorption  and  Raman  measurements  of  highly  porous,  optical  quality,  self  supporting  films 
made  from  p-type  Si  wafers,  permitted  us  to  clarify  the  characteristics  of  the  nanostructure  of  the 
samples.  This  nanostructure  appeared  to  be  an  array  of  spherical  nanocrystallites  with  a  lognormal  size 
distribution.  The  maximum  of  the  distribution  is  dependent  on  the  preparation  parameters  but  is 
typically  around  ~2.5nm  for  the  optical  quality  films. 

•  The  analysis  of  picosecond  nonlinear  optical  response  of  PSi  samples  showed  a  broad  and  intense 
absorption  band  in  the  UV-visible  region  with  a  cut-off  edge  at  ~400nm.  The  non-linear  optical 
response  consisted  of  an  instantaneous  component  followed  by  a  long-lived,  slowly  decaying 
component.  The  former  is  associated  with  the  third-order  susceptibility  of  the  material  whereas  the  later 
originated  from  the  contribution  of  surface  states.  The  response  from  surface  states  was  significant  and 
with  potential  in  practical  applications. 

•  Silicon  Based  UV  detector  prototypes  using  luminescent  PSi  were  prepared  by  growing  PSi  overlayers 
on  top  of  a  n-p  junction.  The  UV  photoresponse  at  366nm  was  enhanced  as  compared  with  control 
prototypes. 

•  Photoluminescent  nanocrystalline  silicon  particles  embedded  in  Si02  were  synthesized  by  RF  co¬ 
sputtering  of  Si  and  Si02  targets.  Stable  photoluminescence  ranging  from  green  to  deep  red  was 
obtained  from  samples  exposed  to  air. 

•  Optical  absorption  and  Raman  analysis  of  Si/Si02  composites  revealed  nanoparticles  size  distributions 
peaking  in  the  range  between  l.lnm  to  3.0nm.  The  measured  photoluminescence  related  with  each  size 
distribution  was  in  agreement  with  the  corrected  local  density  approximation  calculation. 

•  Si/Si02  nanocomposites  were  doped  with  Erbium  and  Europium  and  their  photoluminescence 
properties  were  analyzed  in  relation  with  the  nanostructure. 

•  The  intensity  of  the  characteristic  1 .54pm  emission  from  Er+3  was  found  to  be  dependent  on  the 
nanoparticles  size  and  density.  In  particular,  we  observed  that,  in  the  range  of  the  Si-particles  sizes  we 
studied,  larger  particles  produced  larger  IR  emission.  It  was  also  observed  that  the  use  of  Er203  for  the 
synthesis  of  the  doped  composite  could  introduce  large  amounts  of  oxygen  and  prevent  the  formation  of 
the  silicon  nanocrystals. 

•  We  were  able  to  incorporate  Eu+3  into  the  composite.  We  also  found  that  we  obtained  Eu+2  in  the 
process.  The  Eu+3  related  photoluminescence  peaks  were  only  observed  where  the  size  distribution 
maxima  of  the  nanoparticles  were  lower  than  1 ,3nm. 
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Abstract 

Measurements  of  photoluminescence,  electronic  surface  states  and  effective  surface  area,  at  various  stages  of  the 
anodization  process,  are  presented.  The  results  show  a  behavior  different  from  that  observed  on  n-type  porous  Si. 

Keywords :  Porous  silicon;  Luminescence;  Surface  states 


1.  Introduction 

Porous  silicon  [1]  (PS)  exhibits  high  lumines¬ 
cence  efficiencies  in  the  visible  range.  The  reasons 
for  the  high-efficiency  luminescence  are  still  under 
debate  [1,2].  It  was  suggested  that  the  phenom¬ 
enon  is  due  to  the  amorphous  or  microcrystalline 
nature  of  the  PS  or  that  the  formation  of  silicon 
compounds  are  involved  in  the  luminescence.  To 
gain  some  insight  into  this  process  we  employed 
combined  studies  of  the  luminescence  spectrum, 
and  the  electronic  surface-state  density  and  effec¬ 
tive  surface  area  of  the  porous  surface. 


2.  Experimental 

The  starting  material  was  high-grade  p-type  sili¬ 
con  wafers  of  resistivity  in  the  range  30-50  £2  cm. 
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A  p+  layer  was  formed  by  diffusing  A1  into  one  of 
the  faces  to  attain  an  ohmic  contact.  The  sample 
was  etched  in  20%  HF  and  then  anodized  in  a  solu¬ 
tion  of  HF,  ethanol  and  water  (1:1:2)  with 
a  current  density  of  100  mA/cm2.  The  PS  was 
studied  at  different  stages  of  the  anodization  pro¬ 
cess.  The  luminescence  was  excited  by  a  10  mW 
He-Cd  laser,  X  =  442  nm.  The  anodizing  solution 
was  replaced  after  each  anodization  stage  by  an 
indifferent  electrolyte,  an  aqueous  solution  of  KC1. 
The  electronic  measurements  were  based  on  the 
‘pulsed  capacitance’  method  at  the  semiconduc¬ 
tor-electrolyte  interface  and  have  been  described 
previously  [2]. 


3.  Results  and  discussion 

In  Fig.  1  we  present  typical  luminescence  spectra 
measured  after  different  anodization  times,  as 
marked  on  the  spectra.  We  notice  that  both  the 
amplitude  and  the  wavelength  of  the  maxima 
change  as  a  function  of  anodizing  time.  This  is 
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Fig.  1.  Luminescence  spectra  for  various  anodization  times  as 
marked. 


illustrated  in  Fig.  2  where  we  plotted  the  amplitude 
and  the  wavelength  of  the  maxima  as  a  function  of 
anodizing  time.  As  can  be  seen,  the  luminescence 
spectrum  attains  a  maximum  after  ^  4  min  anodiz¬ 
ation,  after  which  it  decreases.  This  is  similar  to  the 
results  found  [3]  on  n-type  PS.  We  measured  also 
the  effective  surface  area  and  the  density  of  surface 
states.  We  found  that  the  effective  surface  area 
increases  upon  anodization  fairly  monotonously  by 
a  factor  of  15^10,  while  the  surface-state  density 
increases  from  -  2  x  1012  cm"2  for  the  ‘virgin’  sur¬ 
face  to  -  1013  cm'2  for  the  anodized  surface.  This 
increase  occurs  already  after  2  min  anodization  and 
the  density  stays  fairly  constant  upon  further  anod¬ 
ization.  This  is  quite  contrary  to  n-type  PS,  where 
we  have  found  [3]  that  the  surface-state  density 
attains  a  minimum  when  the  luminescence  reaches 
its  maximum.  The  discrepancy  between  the  two 
results  can  be  easily  reconciled  once  we  realize  that 
in  both  cases  the  surface  states  are  measured  under 
accumulation  condition.  The  surface  states  in¬ 
volved  are  then  those  near  the  majority  carrier 
band  edge;  namely,  near  the  conduction  band  for 


Fig.  2.  Peak-luminescence  intensity  and  wavelength  as  func¬ 
tions  of  anodization  time. 


n-type  and  valence  band  for  p-type.  Thus  our  pres¬ 
ent  results  suggest  that  the  surface  states  near  the 
valence-band  edge  are  not  involved  in  the  lumines¬ 
cence  process. 
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Abstract 


We  present  a  summary  of  a  comprehensive  study  of  transport,  phototransport,  surface  photovoltage  spectroscopy  and 
luminescence  of  self-standing  films  of  mesoporous  and  nanoporous  silicon.  We  conclude  that  the  visible  luminescence 
arises  from  optical  bandtail  transitions  in  the  silicon  crystallites,  while  the  transport  is  via  a  disordered  silicon  shell  that 
wraps  the  crystallites. 

Keywords:  Porous  silicon;  Transport;  Phototransport;  Surface  photovoltage 


It  has  been  recently  shown  that  visible  photo¬ 
luminescence  (PL)  in  silicon  can  result  from  a  quan¬ 
tum  confinement  (QC)  effect.  This  and  the  findings 
in  porous  silicon  (PSi)  provide  convincing  evid¬ 
ence  that  the  strong  red  PL  in  PSi  is  due  to  QC 
effects.  The  next  issue  that  needs  to  be  resolved 
is  the  electroluminescence,  i.e.  the  understanding 
of  the  transport  and  the  phototransport  (PT)  in 
PSi  and  their  relation  to  the  luminescence.  Rela¬ 
tively,  this  issue  has  not  been  intensively  studied 
thus  far. 

Accordingly,  we  have  carried  out  a  comprehens¬ 
ive  study  of  the  transport,  PT  [1,  2],  surface  photo¬ 
voltage  [3],  and  PL  [4]  in  ‘free-standing’  PSi  films. 
In  particular,  we  have  studied  mesoporous  films. 
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where  QC  effects  are  negligible,  and  nanoporous 
films,  where  QC  effects  are  dominant.  In  this  paper 
we  summarize  our  results  [1-4]  and  derive  a  self- 
consistent  picture  that  accounts  for  them. 

•  All  the  samples  were  fabricated  from  boron- 
doped  c-Si  (100)  wafers.  In  order  to  eliminate  the 
contribution  of  the  c-Si  substrate  to  the  surface 
photovoltage  and  the  electrical  conductivity, 
‘free-standing’  films  were  required.  Mesoporous 
films  were  prepared  from  wafers  of  0.005  fi  cm  [1]. 
Using  the  Raman  spectra,  the  average  crystal¬ 
lite  size  d  was  established  to  be  10  nm.  This  and 
the  weak  PL  that  is  peaked  at  1.05  eV,  as  in 
c-Si,  determine  the  ‘mesoporous  nature’  of  these 
samples. 

Three  sets  of  nanoporous  samples  were  prepared 
from  c-Si  wafers  of  higher  resistivity  by  using  higher 
current  densities  [2-4].  The  Raman  spectra  re¬ 
vealed  that  d  %  3-4  nm  [1],  i.e.  the  ‘nanostructure 
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nature’  of  the  films.  The  Raman  spectra  also 
contain  features  which  indicate  the  presence  of  dis¬ 
order,  as  in  hydrogenated  amorphous  silicon  (a- 
Si:H).  In  the  surface  photovoltage  spectroscopy 
(SPS)  measurements,  it  was  found  that  the  films 
of  one  set  (set  A)  were  of  p-type,  the  second  set 
(set  B)  of  n-type,  and  the  third  set  (set  C)  of  ‘weak’ 
p-type  [3]. 

Our  main  experimental  results  were: 

(1)  Efficient  (rj  ~  1%)  PL  in  films  of  set  A  [4]. 
This  PL  peaks  at  a  photon  energy  hv  of  1.7  eV.  The 
peak  of  the  emission  of  sets  B  and  C  was  at  the 
same  hv,  but  r\  was  ten  times  smaller  for  set  C,  and 
was  hardly  detectable  in  set  B.  On  the  other  hand, 
the  photoconductivity  in  set  B  was  ten  times  larger 
than  that  of  set  C  and  it  was  undetectable  in  set  A. 

(2)  The  dark  conductivity  in  all  samples  was  orders 
of  magnitude  smaller  than  that  of  the  wafer  sub¬ 
strates.  The  dark  conductivity  activation  energy 
is  larger  than  half  the  optical  bandgap  or  half 
the  energy  of  the  PL  photons.  If  we  use  the  ‘micro¬ 
scopic  prefactor’  of  a-Si:H,  i.e.  150  Q-1  cm”1,  the 
energy  separation  between  the  conduction  ban- 
dedge  and  the  Fermi  level  is  found  to  be  0.5  eV,  in 
both  mesoporous  and  nanoporous  materials.  The 
current-voltage  (I-V)  characteristics  were  all 
linear.  It  was  found  that  the  presence  of  adsor¬ 
bed  oxygen  and/or  water  vapour  reduces  the 
luminescence,  enhances  the  conductivity  and  the 
photoconductivity  and  renders  the  films  n-type. 

(3)  The  diffusion  length  of  the  electrons  in  the 
nanoporous  n-type  flakes  is  of  the  order  of  1  pm, 
while  that  of  the  holes  is  about  0.2  pm.  (4) 
The  temperature  dependences  of  the  dark  and 
photoconductivities  roughly  resemble  those  found 
in  a-Si:H.  (5)  The  width  of  the  photovoltage  band- 
to-band  absorption  edge  is  0.3  eV  for  both  n-type 
and  p-type  samples. 

In  view  of  the  limited  scope  of  this  paper  we  only 
highlight  the  conclusions  drawn  from  the  above 
new  experimental  findings.  Based  on  these  results 
we  derive  the  following  self-consistent  picture. 
From  the  findings  (1)  and  (2)  it  appears  that  the 


conduction  is  due  to  transport  in  a  disordered  Si 
network,  and  is  not  due  to  tunneling  between  crys¬ 
tallites.  From  finding  (3)  it  appears  that  the  con¬ 
cepts  that  apply  to  bulk  semiconductors  also  apply 
to  PSi,  in  spite  of  its  coral  structure.  From  findings 
(3)  and  (4),  it  appears  that  the  current  flows  through 
the  a-Si :  H-Iike  shell  that  wraps  the  crystallites. 
From  findings  (1)  and  (5)  it  appears  that  the  domi¬ 
nant  optical  absorbtion  takes  place  in  the  crystal¬ 
lites.  Hence,  there  are  both  conduction  and  valence 
bandtails  in  the  crystallites,  which  are  about  0.3  eV 
wide.  The  origin  of  these  tails  and  of  the  lumines¬ 
cence  is  mainly  due  to  the  relaxation  of  the  lattice 
momentum  selection  rules.  The  PL  is  due  to  radi¬ 
ative  transitions  between  the  crystallites’  bandtails. 
The  high  efficiency  of  the  PL  is  due  to  the  combina¬ 
tion  of  this  relaxation  and  the  geometrical-poten¬ 
tial  confinement  of  the  carriers  in  the  crystalliltes. 
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ABSTRACT 

Pulse  measurements  on  the  porous-Si/electrolyte  system  are  employed  to  determine  the 
surface  effective  area  and  the  surface-state  density  at  various  stages  of  the  anodization  process 
used  to  produce  the  porous  material.  Such  measurements  were  combined  with  studies  of  the 
photoluminescence  spectra.  These  spectra  were  found  to  shift  progressively  to  the  blue  as  a 
function  of  anodization  time.  The  luminescence  intensity  increases  initially  with  anodization 
time,  reaches  a  maximum  and  then  decreases  with  further  anodization.  The  surface  state  density, 
on  the  other  hand,  increases  with  anodization  time  from  an  initial  value  of -2x10 12  cm-2  for  the 
virgin  surface  to  ~1013  cm*2  for  the  anodized  surface.  This  value  is  attained  already  after  -2  min 
anodization  and  upon  further  anodization  remains  fairly  constant.  In  parallel,  the  effective 
surface  area  increases  by  a  factor  of  10  -  30.  This  behavior  is  markedly  different  from  the  one 
observed  previously  for  n-type  porous  Si. 


INTRODUCTION 

Porous  silicon,1*4  (PS)  obtained  by  electrochemical  etching  procedures  applied  to  crystalline 
Si  surfaces,  exhibits  high  luminescence  efficiencies  in  the  visible  range.  It  is  quite  clear  now  that 
the  visible  luminescence  originates  from  the  band-gap  enlargement  due  to  quantum 
confinement 4  At  the  same  time,  the  reasons  for  the  high-efficiency  luminescence  are  still 
somewhat  under  debate.4-5  It  was  suggested  that  it  is  the  amorphous6  or  microcrystalline7  nature 
of  the  porous  Si  that  is  responsible  for  the  phenomenon,  or  that  the  formation  of  silicon 
compounds  such  as  siloxene  (SigO^Hg)  or  species  of  Si-H,  Si-0  and  Si-F  bonds  are  involved  in 
the  luminescence.8  One  way  of  gaining  further  insight  into  the  luminescence  process  is  to  carry 
out  a  variety  of  measurements  on  samples  of  different  porosity.  To  that  end  we  have  employed 
combined  studies  of  the  luminescence  spectrum,  the  surface-state  density  and  the  effective 
surface  area  of  the  porous  surface.  Such  studies  were  carried  out  at  different  stages  of  the 
anodization  process  and  thus  for  different  morphologies  of  the  porous  surface.  The  luminescence 
spectra  were  measured  by  conventional  methods.  The  surface  state  characteristics  and  the 
effective  surface  area  were  determined  by  pulse  measurements9  on  the  PS/electrolyte  system. 
This  system  is  particularly  suitable  since  a  capacitative  contact  to  the  terrain  of  the  porous 
surface  is  best  achieved  by  an  electrolyte,  and  it  was  successfully  used10  to  investigate  n-type  PS. 
There  we  found 1(1  a  strong  correlation  between  the  surface-state  density  near  the  conduction  band 
edge  and  the  luminescence  intensity.  In  this  paper  we  present  similar  measurements  on  p-type 
porous  Si  and  we  compare  the  results  with  those  obtained  on  n-type  material. 


EXPERIMENTAL 

The  starting  material  was  high-grade  p-type  silicon  wafers  of  resistivity  in  the  range  20  -  50 
Qcm.  A  p+  layer  was  formed  by  diffusing  metallic  A1  into  one  of  the  faces  to  obtain  an  ohmic 

461 

Mat  Res.  Soc.  Symp.  Proc.  Vol.  452  °  1997  Materials  Research  Society 


^sso  occupied  surface  states.  Vcq  and  P^  can  be  determined  quite  accurately  from 
measurements  in  the  depletion  range.1*  The  entire  plots  of  Ps  and  P^  vs.  Vs  can  then  be 
constructed  by  using  the  relations  Vs  =  Vso  +  8VS,  Ps  =  Pso  +  8PS,  and  P^  =  P^  +  5PSS.  So 
much  so  m  the  absence  of  an  insulating  buffer  layer  (such  as  an  oxide)  at  the  semiconductor 
surface.  If  such  a  layer  is  present,  the  as-measured  barrier  height,  i.  e.,  the  measured  voltage  drop 
between  the  Pt  electrode  and  the  p+  contact  just  after  the  pulse  termination,  yields  Vs+Vc,  where 
Vg  is  the  voltage  drop  across  the  insulating  layer.  * 


RESULTS  AND  DISCUSSION 

In  Fig.  1  we  present  typical  photoluminescence  spectra  of  p-type  PS  surfaces  prepared  by 
anodization.  The  different  anodization  times  are  marked  on  the  spectra.  We  notice  that  the 
luminescence  intensity  at  the  beginning  increases  with  anodization  time,  attains  a  maximum  and 
then  decreases.  This  behavior  is  illustrated  by  the  higher  curve  in  Fig.  2  and  is  quite  similar  to 
that  observed10  previously  for  n-type  PS.  The  lower  curve  in  Fig.  2  shows  the  appreciable  blue 
shift  of  the  spectra,  suggesting  that,  on  the  average,  the  porous  structure  gets  finer  with 
anodization  time. 


500550600650700  750  000850 
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Fig.l.  Photoluminescence  spectra  for  various 
anodization  times,  as  marked. 


Fig.  2.  Peak-photoluminescence  intensity 
and  wavelength  as  functions  of  anodization 
time  for  the  sample  in  Fig.  1 . 


Typical  results  of  the  free  surface-hole  density  Ps  (diamonds)  and  the  density  of  occupied 
surface  states  Pss  (stars)  against  the  as-measured  barrier  height  Vs+Vg  as  obtained  for  an  etched 
virgin  silicon  surface,  are  displayed  in  the  semilog  plot  of  Fig.  3.  In  the  depletion  range,  Ps  is 
negative  but,  because  of  the  logarithmic  scale  used,  the  plot  is  that  of  -Ps.  The  light  curve,  in  the 
accumulation  range,  labeled  Cg  =  oo,  represents  the  theoretical  dependence  of  Ps  on  Vs  for  a 
buffer-free  surface  (Cg  =  oo,  Vg  =  0),  as  derived  from  a  solution  of  Poisson’s  equation  for  the 
value  of  the  hole  bulk  concentration  pj,  marked  in  the  figure.  It  is  seen  that  this  curve  does  not 
account  well  for  the  data  in  the  accumulation  range.  The  best  fit,  represented  by  the  bold  curve, 
labeled  Cg  -  4.5  pF/cm2,  was  obtained  by  assuming  the  presence  of  an  insulating  buffer  layer  of 
Cg  =  4.5  pF/cm2,  corresponding  to  a  thickness  of  --4  A.  Most  probably  this  layer  consists  of  an 
oxide,  but  a  small  contribution  of  the  series  Helmholtz  layer  cannot  be  ruled  out.  Again,  this 
behavior  is  similar  to  that  found  on  n-type  PS.  Turning  now  to  the  surface-state  hole  occupancy 
Pss,  it  is  seen  to  rise  slowly  from  a  low  value  of -101 1  cm"2  at  Vs+Vg  « -0.5  V  up  to  ~1014  cnr2 
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Surface  Hole  Density  (cm -2) 


Fig.  3.  Free  surface  hole  density  Ps  (diamonds) 
and  density  of  Qccupied  surface  states  (stars) 
vs.  the  as-measured  barrier  height  Vs+Vg  for  a 
virgin  Si  surface.  The  light  and  bold  curves  are 
theoretical  plots  of  Ps  as  explained  in  the  text. 


Fig.  4.  Ps  (diamonds)  and  Pss  (stars)  vs.  Vs+Vg 
for  four  anodized,  porous  Si  surfaces.  The  light 
and  bold  curves  are  theoretical  plots  of  Ps. 


at  Vs+Vg  »  -2  V.  No  saturation  value  for 
the  surface-state  density  could  be 
reached.  This  is  quite  different  than  the 
behavior  we  found10  on  n-type  PS,  where 
the  surface-state  density  saturated  at 
IO12  cm’2.  Because  of  this  lack  of 
saturation,  we  shall  use  the  values  of  P^  at 
Vs+Vg  «  -1  V  to  compare  surface  state 
occupancies  for  different  anodization  times. 
For  the  case  of  the  virgin  surface,  this  value 
is  ~2xlOl2cm’2. 

Results  of  Ps  and  P^,  similar  to  those  in 
Fig.  3,  for  four  porous  surfaces  are 
presented  in  Fig.  4.  These  results  were 
obtained  after  the  sample  of  Fig.  3  has  been 
anodized  for  different  times,  as  marked  in 
the  figure.  Since  our  aim  is  to  compare  the 
surface-state  densities,  we  show  only  the 
accumulation  range  here.  Because  of  the 
increase  of  the  effective  surface  area,  the 
highest  surface  potential  barriers  (for  holes) 
attained  were  around  -1  volt.  The  curves  in 
the  figure  are  theoretical,  calculated  for  the 
same  Cg  values  as  in  Fig.  3.  We  notice 
again  that  for  all  four  porous  surfaces,  the 
surface-state  density  (stars)  increases 
monotonously  with  the  potential  barrier 
through  the  whole  region  shown  and  does 
not  exhibit  signs  of  saturation.  As 
mentioned  above,  we  choose  for 
comparison  the  values  of  Pss  at  —1  volt. 
These  values  in  the  figure  are  scattered 
around  2x10 13  cm’2,  about  an  order  of 
magnitude  higher  than  on  the  virgin 
surface.  However,  the  interesting  thing  is 
that  this  surface-state  density  remains  fairly 
constant  with  anodization  time,  very  much 
different  from  the  results  found10  for  n- 
type  PS.  A  behavior  similar  to  that  of  the 
surface  states  was  observed  also  for  the 
effective  surface  area.  In  Fig.  5  the  surface 
area  factor,  i.  e.  the  ratio  of  the  effective 
area  to  the  area  of  the  virgin  sample,  for 
two  typical  samples  is  plotted  against  the 
anodization  time.  The  area  factors  were 
derived  from  measurements  of  the  type 
shown  in  Figs.  3  and  4  in  the  strong 
accumulation  range.  The  results  for  the 
two  samples  are  qualitatively  the  same;  the 
area  factor  increases  with  anodization  time 
till  it  reaches  a  saturation.  The  saturation 
values  observed  varied  from  sample  to 
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Fig.  5.  Surface  area  factor  vs.  anodization 
time  for  two  samples. 


Fig.  6.  Photoluminescence-peak  intensity  vs. 
anodization  time  for  two  samples. 


sample  between  10  and  30 

Jr  *«»  & 

PS/electrolyte  system.  In  Fig  6  the  lumine^cnc^01”!^-  after.the  measurements  on  th 
■notation  dJ  to  to  «>  L  tnp”'";/?  SIT'S  *<=*""'  “> 

luminescence  peak  obtains  already  after  ?  min  *•  , e  no*Ice  that  the  maximun 

upon  further  anodization.  In  Fig.  7  we  plot  the  valulfnf  ft!  T  ^en  ^  ^um^nescence  decrease: 
function  of  anodization  time.  The  peak  shifts  to  the  hfi  *  Iununescence-peak  wavelength  as  i 
its  value  at  2  min  anodization.  *  blue  Up°n  anodl2ation  clos*  to  200  nm  frou 

com^rL^S^  (here  fs  TdiIc!ePMcy  £  A ®'  2  7  ”?tice-that  theV  b°*  show 
maximum  luminescence.  To  check  the  influence  of  fi°n  tlme  needed  *°  attain 

measurements,  we  measured  the  spectra  on  PS  samnlcftt.  ^C*  electroIyte  used  in  the  electronic 
10  min  into  a  KCI  electrolyte  iffSte  S,  L  p  7  T?,Zation  were  Versed  for 
plotted  the  corresponding  curve  from  Fig.  1.  We  notice  that  ' the  KCI We  re' 
maximum  to  a  different  anodization  time  however  it  Uli  K.u,nt  does  not  shift  the 
This  latter  is  probably  due  to  adsorotion’of  some  ’  *  d  ■  ?Wer  ^  luminescence  intensity, 
ascribe  the  adsorption  of  some  species,  from  the  electrolyte.  Thus  we 


Fig.  7.  Photoluminescence-peak  wavelength 
vs.  anodization  time  for  two  samples. 


Fig.  8.  Photoluminescence-peak  intensity  for 
an  untreated  and  for  a  KCI  treated  sample  vs. 
anodization  time. 
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shift  of  the  anodization  time  for  maximum  luminescence,  to  possible  surface  damage  due  to  the 
application  of  the  voltage  pulses. 


CONCLUSION 

P-type  PS  behaves  quite  differently  from  that  of  n-type  PS.  While  the  photoluminescence 
intensity  exhibits  a  pronounced  maximum  with  anodization  time,  the  effective  surface  area  and 
the  surface  state  density  appear  to  reach  a  more  or  less  constant  value  as  a  function  of 
anodization  time.  This  is  in  contrast  to  n-type  PS,  where  a  close  correlation  between  the 
effective  surface  area,  the  luminescence  intensity  and  the  surface-state  density  was  found.10  The 
discrepancy  between  the  two  results  can  be  reconciled  once  we  realize  that  in  both  cases  the 
surface  states  are  measured  under  accumulation  conditions.  The  surface  states  involved  are  then 
those  near  the  majority  carrier  band  edge;  namely,  near  the  conduction  band  for  the  n-type 
material  and  near  the  valence  band  for  the  p-type  material.  Thus  our  present  results  suggest  that 
the  surface  states  near  the  valence-band  edge  are  not  involved  in  the  luminescence  process. 
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ABSTRACT 

We  have  studied  the  four  phototransport  properties  as  a  function  of  temperature  in 
undoped  a-Si:H  films  deposited  with  different  substrate  temperatures  in  the  range  150  -  225°C. 
The  analysis  of  the  results  indicates  how  T,  determines  the  densities  of  the  various  defects.  The 
general  trend  of  decrease  of  both  the  density  of  dangling  bonds  and  the  Urbach  energy  is  in 
agreement  with  the  weak  bond  breaking  model.  However,  we  conclude  that  a  sbght  modification 
of  this  model  is  required. 


INTRODUCTION 

One  of  the  most  important  parameters  in  the  deposition  of  a-Si:H  films  and  devices  is  the 
substrate  temperature,  T,.  Previous  works  have  shown  that  the  variation  of  T,  effects  the 
Uibach  energy  (Eu)  and  the  density  of  dangling  bonds  [1],  degree  of  disorder  12],  concentration 
of  Hydrogen  in  the  film  [3]  and  the  phototransport  properties  14].  The  latter,  and  in  particular 
the  mobility-lifetime,  (IT,  product  of  the  minority  carriers,  are  of  particular  importance  for  the 
production  of  high  quality  bipolar  devices,  such  as  solar  cells  [5],  Despite  that  importance  and 
die  basic  interest  in  the  properties  of  the  material,  the  understanding  of  the  exact  relation 
between  the  phototransport  properties  and  the  nature  and  concentration  of  the  various  defects  is 
Kill  at  a  rudimentary  stage.  It  is  thus  still  unclear  to  which  extent  the  minimization  of  material 
properties  such  as  Eu  and  the  density  of  the  dangling  bonds  really  leads  to  the  optimization  of 
the  phototransport  properties. 

Recent  works  [6-8]  have  proposed  that  the  recombination  processes  in  a-Si.H  are 
determined  by  an  interplay  between  the  valence  band  tail  states  and  the  dangling  bonds. 
Specifically,  it  was  suggested  that  at  room  temperature  the  phototransport  of  the  majority 
carriers  (i.e.  the  electrons)  is  determined  by  the  dangling  bonds,  while  the  phototransport  of  the 
minority  earners  (the  holes)  is  determined  by  the  valence  band  tail  states  [7].  In  this  paper  we 
examine  the  validity  of  this  suggestion  by  directly  measuring  the  phototransport  properties  as  a 
function  of  the  temperature  in  a  series  of  a-Si:H  films.  The  only  difference  in  the  way  these  films 
were  prepared  was  the  T,  under  which  the  deposition  took  place.  The  properties  which  were 
measured  include  the  jiT  products  of  the  two  carriers  and  the  exponents  that  characterize  their 
yu  intensity  dependence  (y  and  S  for  the  electrons  and  the  holes,  respectively).  As  we  have 
previously  shown,  this  procedure  provides  quite  a  complete  characterization  of  the 
recorribination  kinetics  [9].  The  determination  of  these  four  phototransport  properties  is  done 
tlrough  the  measurement  of  the  photoconductivity,  and  ambipolar  diffusion  length,  L,  and 

|  the  application  of  the  relations  [5]: 

;  (Hf)«  =  Opt,  /  qG 
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(2) 

(3) 

(4) 


(HX)k  =  qLJ  /  2kT 

Of*  oe  GT 

L«  G"5 

Here  q  is  the  electronic  charge,  G  is  the  photogeneration  rate  and  kT  is  the  thermal  energy. 

As  we  show  below,  the  results  we  obtained  from  the  temperature  dependence  of  these 
quantities  confirm  the  current  understanding  of  the  recombination  kinetics,  yielding  an  excellent 
agreement  with  the  conclusions  of  our  recent  Raman  scattering  study  [2]  of  these  films.  We 
further  show  that  the  results  are  consistent  in.  their  gross  features  with  the  weak  bonds  breaking 
model  [11,  although  further  considerations  should  be  made  in  order  to  account  for  the  slight 
deviations  of  the  experimental  results  from  this  theory. 

EXPERIMENTAL 

The  a-Si:H  films  used  in  this  study  were  prepared  using  the  rf  glow  discharge 
decomposition  of  silane  technique.  The  substrates  were  made  of  Coming  7059  glass.  The 
substrate  temperatures  during  deposition  were  150,  175,  200  and  225°C  (samples  1-4, 
respectively),  and  the  deposition  rate  was  about  4  A/sec.  The  thickness  of  all  the  films  was 
about  0.8  pm.  The  a-Si:H  films  were  deposited  on  top  of  predeposited  NiCr  layers,  which 
provided  coplanar  ohmic  contacts,  with  a  0.4  mm  separation. 

The  phototransport  properties  were  determined  using  the  photocarrier  grating  (PCG) 
technique  [5,10],  This  technique  yields  the  phototransport  properties  of  the  two  charge  carriers 
by  the  simultaneous  measurement  of  Op*  and  L,  followed  by  the  application  of  equations  (1)  • 
(4).  The  light  source  in  our  measurements  was  a  He-Ne  laser  with  a  wavelength  of  633  nm.  The 
maximum  illumination  flux  used  was  10  mW/cml,  corresponding  to  a  photogeneration  rate  of  0 
-  2*10"  cm'W).  For  the  determination  of  the  light  intensity  exponents  we  have  reduced  the 
above  photogeneration  rate  in  steps,  down  to  101*  cmJsecl.  The  samples  were  placed  ioi 
cryostat,  which  enabled  us  to  cany  out  the  measurements  in  the  temperature  range  of  180  • 
350°  K. 

EXPERIMENTAL  RESULTS 

Following  the  above  outlined  procedure  we  have  measured  the  four  phototranspoit 
properties,  for  all  four  samples,  as  a  function  of  the  temperature.  All  the  results  are  presented  a 
Fig.  1.  Starting  with  the  majority  carriers  we  present  in  Fig.  1(a)  the  temperature  dependence 
of  (px)«.  It  can  be  seen  that  for  sample  1  (T,  =  150°C)  (pt),  increases  monotonically  with 
temperature,  and  its  values  are  considerably  lower  (by  more  than  an  order  of  magnitude)  that  ! 
those  obtained  for  all  other  samples.  The  temperature  dependence  of  (HT)e  in  the  other  films  is 
more  complicated:  At  first  (pt),  increases  sharply  with  temperature,  and  then  at  a  much  stows 
pace.  This  is  a  manifestation  of  the  well  known  thermal  quenching  (TQ)  of  Of*  [6,11].  He 
temperatures  at  which  the  TQ  occurs  (Ttq)  are  about  200,  220  and  250  K  for  samples  2  *  4 

respectively.  . 

It  is  well  known  that  the  TQ  is  accompanied  by  a  cusp  in  the  temperature  dependences 1 
corresponding  exponent,  y  [6.8,11],  In  Fig.  1(b)  it  can  be  seen  that  indeed  such  peaks « 
observed  for  samples  2-4,  thus  confirming  the  above  explanation  of  the  results  of  Fig.  I(i) ' 
We  can  also  see  in  Fig.  1(b)  that  for  sample  1  y  is  monotonically  decreasing.  This,  and  Ik 

.t 
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tndofTn),  suggest  that  for  sample  1  the  TQ  occurs  at  a  temperature  which  is  lower  than  the 
operimenta]  range  accessible  in  our  measurements. 


fig.  1:  The  measured  temperature 
dependence  of  the  phototransport 
properties  for  samples  1(T,  =  I50°C, 
circles).  2  (T,  =  175°C,  triangles).  3  (T, 
*200°C.  diamonds)  and  4  (T,  =  225°C. 
l^uies).  The  smooth  lines  are  drawn  to 
guide  the  eye. 


Turning  to  the  minority  carriers,  we  present 
in  Fig.  1(c)  the  temperature  dependence  of  (px)*  for 
our  four  samples.  We  can  see  th^t  for  ail  the 
samples  (px)h  increases  monotonically  with  T. 
However,  the  T,  dependence  of  (px)h  is  not 
monotonic  as  that  of  (px)e:  For  samples  1  -  3  we 
observe  an  increase  of  (ju)h  with  T„  while  for 
sample  4  the  (px)h  values  are  tower  than  those  of 
sample  3,  indicating  that  T,  =  200°C  is  an  optimal 
value  for  obtaining  a  high  value  of  (pxV  Finally,  we 
present  in  Fig.  1(d)  the  temperature  dependence  of 
the  light  intensity  exponent  of  the  minority  carriers, 
S.  There  are  two  characteristics  to  be  noted  in  the 
behavior  of  S:  First,  S  has  a  non-zero  value 
throughout  the  temperature  range  of  this  study, 
indicating  that  at  least  two  different  types  of  defects 
are  involved  in  the  recombination  process  [9]. 
Second,  for  all  samples  S  is  monotonically 
decreasing  with  temperature. 

In  order  to  summarize  the  influence  of  T*  on 
the  phototransport  properties  as  described  above, 
we  plot  the  T,  dependence  of  the  two  jix  products 
(as  measured  at  250°K)  in  Fig.  2(a),  and  of  the  Ttq 
in  Fig.  2(b).  Since  we  did  not  find  experimentally 
the  value  of  Ttq  for  sample  1  we  have  plotted  an 
extrapolated  dotted  line  for  T,  values  below  175°C. 

DISCUSSION 

For  the  interpretation  of  the  experimental 
results,  let  us  first  review  the  recombination  model 
which  is  commonly  used  for  the  analysis  of 
phototransport  properties.  We  consider  a  recent 
version  of  this  model  which  was  analyzed  in  detail 
by  Tran  [6],  The  model  consists  of  a  valence  band 
tail  which  is  roughly  twice  wider  than  the 
conduction  band  tail,  and  a  single  type  of  dangling 
bond  centers,  which  are  assumed  to  be  neutral  in 
equilibrium.  Accordingly,  there  are  a-priori  three 
possible  recombination  channels.  However,  since 
the  capture  cross  sections  of  the  states  in  the  two 
band  tails  arc  assumed  to  be  symmetric,  the  higher 
total  density  of  valence  band  tail  states  causes  the 


211 


recombination  rate  via  this  taiJ  to  dominate  the  recombination  rate  via  the  conduction  band  tail  at 
all  temperatures.  The  temperature  dependence  of  the  phototransport  properties  is  determined 
then  by  the  recombination  interplay  between  the  valence  band  tail  states  and  the  dangling  bonds. 

Recent  publications  have  shown  quite  convincingly  that  at  low  temperatures  the  majority 
carriers  lifetime  is  determined  by  the  valance  band  tail  [6,8].  As  the  temperature  increases,  the 
thermal  release  of  trapped  holes  causes  a  decrease  of  the  recombination  rate  through  this 
channel,  and  an  increase  in  the  recombination  rate  through  the  dangling  bonds.  This  process 
continues  until  at  high  enough  temperatures  the  electrons  lifetime  is  determined  solely  by  the 
dangling  bonds.  Naturally,  there  is  a  temperature  at  which  the  recombination  rate  through  the 
two  channels  is  about  equal,  and  Tran  has  shown  |6]  that  this  is  the  temperature  at  which  the 
thermal  quenching  takes  place  (i.e.  Ttq).  The  TQ  is  caused  then  by  a  transition  between 
temperature  regions  where  different  channels  dominate  the  recombination  process.  For  lower 
light  intensity  a  higher  temperature  is  required  for  such  a  transition,  leading  to  the  high 
sensitivity  of  (px)e  to  the  light  intensity  at  temperatures  around  TTq.  This  is  manifested  by  the 
cusp  in  y(T).  Furthermore,  the  existence  of  two  compatible  recombination  channels  enables  that 
cusp  to  exceed  unity,  contrary  to  a  situation  where  only  one  recombination  channel  is  operative 
(a  situation  which  leads  to  y  <  1  [9]). 


Ts  (K) 


Fig.  2:  The  effect  of  T,  on  (a)  the  values  of  the  two 
px  products  (taken  at  250°K)  and  on  (b)  the  thermal 
quenching  temperature. 


Turning  to  the  minority 
carriers,  it  is  expected  that  their 
lifetime  will  increase  monotonically 
with  temperature  [6].  Based  on  the 
measurement  of  (pt)b  as  a  function  of 
light  soaking  time  Wang  and  Schwarc 
[7]  concluded  that  at  room 
temperature  the  holes  lifetime  is 
determined  by  the  valance  tail  states, 
while  the  electrons  lifetime  is 
determined  by  the  dangling  bonds  (Le. 
that  Ttq  is  lower  than  room 
temperature).  According  to  their 
analysis  the  dominance  of  the  holes 
lifetime  by  the  dangling  bonds  will  be 
achieved  only  above  350°1C  Recently 
we  have  shown  [8]  that 
correspondingly  the  exponent  S  should 
decrease  with  increasing  T  from  a  non 
zero  value,  which  is  related  to  the 
recombination  process  in  the  valance 
band  tail,  to  a  value  of  zero,  which  is  { 
related  to  the  recombination  process  io  j 
the  dangling  bonds. 

With  all  this  in  mind,  we  can 
now  see  that  all  the  experimental  : 
results  in  Fig.  1  are  in  agreement  with  ; 
this  model  The  TQ  of  is  indeed 
accompanied  by  a  peak  in  y  with 
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values  higher  than  unity.  The  minority  carrier  jit  product  increases  and  the  exponent  S 
decreases  with  increasing  temperature.  It  should  be  emphasized,  however,  that  since  S  does  not 
reach  a  zero  value  for  any  of  the  samples,  we  must  conclude  that  there  is  no  range  of  the 
temperature  interval  under  study  in  which  the  dangling  bonds  control  solely  the  holes 
recombination  process.  Hence,  the  current  results  provide  further  support  to  Tran’s 
recombination  model. 

In  view  of  that  we  can  obtain  spectroscopic  information  regarding  the  influence  of  T,  on 
the  density  of  the  various  defects.  We  saw  in  both  Fig.  1(a)  and  Fig.  2(a)  that  (px)e  increases 
with  the  increase  of  T,.  This  means,  in  accordance  with  previous  suggestions  (1],  that  there  is  a 
reduction  of  the  total  defect  density  as  T,  increases  within  the  T,  interval  under  study. 
Furthermore,  following  the  analysis  of  Tran  [6],  we  conclude  that  the  increase  of  the  Tn?  with  T, 
(Fig.  2(b))  indicates  that  the  ratio  of  the  concentrations  of  dangling  bonds  to  the  valence  band 
tail  states  also  decreases  with  T,.  Hence,  the  new  conclusion  of  this  study  is  that  the  most 
important  effect  of  the  increase  of  T,  is  a  reduction  in  the  concentration  of  the  dangling  bonds. 
While  the  electrons  recombination  is  determined  by  both  the  dangling  bonds  and  the  tail  states, 
the  holes  recombination,  as  discussed  above,  is  dominated  by  the  valence  band  tail  in  the 
temperature  interval  studied  here.  Hence,  the  conclusion  from  Figs.  1(c)  and  2(a)  is  that  as  T, 
increases  from  150°C  to  200°C  the  valence  band  tail  becomes  narrower,  giving  rise  to  the 
increase  of  (nx)*,.  However,  a  further  increase  of  Ts  causes  a  slight  decrease  in  (px)h,  which 
must  be  attributed  then  to  a  rewidening  of  the  valance  band  tail.  It  thus  seems  that  while  the 
reduction  of  the  concentration  of  the  dangling  bonds  is  monotonic  with  T,  this  is  not  the  case  for 
the  valance  band  tail  width,  for  which  a  minimal  value  is  obtained  for  an  optimal  substrate 
temperature.  Indeed  it  is  possible  [1]  that  the  minimum  of  the  Eu  found  here  for  T,  s  200°C  is 
only  a  local  minimum,  and  that  a  further  increase  of  T,  will  results  in  a  further  reduction  of  E0. 
The  current  results  show  then  that  the  relation  between  Eu  and  the  concentration  of  the  dangling 
bond  density  is  not  trivial  (see  below). 

Since  this  discussion  is  based  only  on  the  analysis  of  the  photo  transport  data,  we  now 
turn  to  examine  previous  works  in  which  the  effect  of  Tt  was  studied  by  other  means.  We  have 
recently  presented  a  study  of  these  films  by  Raman  scattering  [2],  and  found  that  the  degree  of 
disorder  in  our  a-Si:H  films  decreases  from  sample  1  to  3,  and  then  increases  slightly  for  sample 
4.  This  is  true  both  for  the  short  range  disorder  (i.e.  width  of  the  bond  angle  distribution)  which 
effects  the  location  and  width  of  the  spectral  band  of  the  transverse  optical  (TO)  mode,  and  for 
the  intermediate  range  disorder  (i.e.  dihedral  angle  distribution)  which  determines  the  ratio  of  the 
intensities  of  the  transverse  acoustic,  TA,  spectral  band  to  the  TO  spectral  band.  Since  the 
disorder  in  the  a-Si:H  is  related  to  Eu  [12]  the  results  of  the  Raman  scattering  are  consistent 
with  the  interpretation  given  above  for  the  phototransport  data.  Furthermore,  this  shows  that 
the  decrease  of  (px)h  with  the  increase  of  T,  from  200  to  225°C  is  not  an  experimental  artifact 
but  rather  a  true  property  of  the  material. 

The  dependencies  of  the  defect  density  in  the  tails  and  in  the  dangling  bonds  on  Tt  have 
been  analyzed  in  detail  by  Stulzmann  [1].  Correspondingly,  he  suggested  the  week  bonds 
breaking  (WBB)  model,  in  which  the  dangling  bonds  originate  from  spontaneous  break  of  weak 
silicon-silicon  bonds  whose  energy  exceeds  a  certain  critical  limit.  Hence,  a  higher  Eu  means 
more  states  above  this  critical  energy,  and  thus  a  higher  density  of  dangling  bonds.  Stulzmann 
has  farther  showed,  by  examining  a  collection  of  experimental  results  (obtained  by  different 
experimental  methods),  that  minimal  values  for  both  the  density  of  dangling  bonds  and  the  Eu 
are  obtained  for  T,  of  about  300°C,  It  is  obvious  then  that  our  finding  of  decreasing  density  of 
dangling  bonds  is  consistent  with  the  collection  of  experimental  data  presented  by  Stutzmann. 
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We  note,  however,  that  the  slight  increase  concluded  here  for  the  Eu  in  sample  4  seems  to  be  in 
disagicemcnt  with  the  WBB  model.  It  seems  then  that  our  results,  while  generally  in  agreement 
with  this  model,  suggest  that  more  refined  details  which  depend  on  T,  may  be  important. 
Examples  for  such  details  are  the  density  of  states  at  the  mobility  edge,  the  concentration  of 
hydrogen,  or  the  energy  of  the  gap.  The  effect  of  these  paranwers  will  be  considered 
elsewhere. 

In  conclusion,  we  have  studied  the  phototransport  properties  as  a  function  of 
temperature  for  a-Si:H  samples  deposited  under  different  substrate  temperatures.  We  have 
shown  that  the  results  are  consistent  with  Tran’s  recombination  model,  and  that  the  changes 
with  T,  are  generally  in  agreement  with  the  weak  bond  breaking  model.  However,  a  refinerwnt 
of  this  model  should  be  considered. 
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The  effects  of  a  controlled  aquaregia  treatment  of  indium-tin-oxide  (ITO)  substrates  on  the 
behavior  of  highly  efficient  vacuum  evaporated  double  layered  8-tris-hydroxyquinoline  A1 
(Alq3) -based  light-emitting  diodes  are  described.  It  is  found  that  in  suitably  treated  devices,  both 
current  injection  and  the  electroluminescence  (EL)  are  significantly  enhanced.  The  enhancement  is 
believed  to  result  from  the  greater  ITO/hole  transporting  layer  contact  areas  and  the  contact 
conditions.  The  observed  dependence  of  /(V),  the  EL  output,  and_the  EL  efficiency  on  the  ITO 
surface  morphology  indicates  that  space-charge-limited  currents  dominate  the  behavior  of  the 
devices.  ©  1997  American  Institute  of  Physics.  [S0003-695 1(97)04020-5] 


Transparent  conducting  indium  tin  oxide  (ITO)  is  used 
widely  as  the  anode  material  in  organic  light-emitting  diodes 
(OLEDs)  due  to  its  efficient  hole  injection  properties.1"5  It  is 
also  generally  recognized  that  the  cleanliness  of  the  ITO  sur¬ 
face  is  an  important  factor  governing  the  performance  of  the 
devices.  Consequently,  many  cleaning  methods  such  as 
scrubbing,  ultrasonication,  chemical  washing,  vapor  degreas¬ 
ing,  irradiation  in  a  UV-ozone  chamber,  etc.,  have  been  de¬ 
veloped  to  treat  its  surface.6-8  Yet  no  studies  of  the  effect  of 
the  ITO  surface  morphology  on  the  electroluminesce  (EL)  of 
OLEDs  have  been  reported.  In  order  to  study  this  effect,  we 
developed  a  simple  method  to  treat  ITO  by  etching  it  with 
diluted  aquaregia.  This  letter  describes  its  effects  on  double¬ 
layered  OLEDs  with  an  8-tris-hydroxyquinoline  A1  (Alq3). 
We  show  that  their  performance  is  strongly  dependent  on  the 
ITO  surface  morphology  and  that  optimized  treatment  yields 
highly  efficient  devices. 

Figure  1  shows  the  structure  of  the  OLEDs,  which  con¬ 
sist  of  the  aquaregia-treated  Applied  Films  Corp.  ITO  sub¬ 
strate  (untreated  sheet  resistance  of  —20  fi/sq),  a  500  A  thick 
hole  transporting  layer  (HTL),  a  500  A  thick  Alq3  layer,  and 
a  1800  A  thick  Mg09Ag01/Ag  layer.  The  HTL  was  either 
a  “three-armed  star”  (3-AS)  4,4',4"-tris[A-(3- 

methoxyphenyl)-Af-phenylamine-triphenylamine]  or  triph¬ 
enyl  diamine  (TPD).  The  molecular  structures  of  Alq3,  3-AS, 
and  TPD  are  shown  in  Fig.  2. 

The  aquaregia  was  prepared  by  adding  HN03  to  HC1 
very  slowly,  to  a  ratio  of  1:3,  and  letting  the  mixture  brew 
for  15  min.  The  aquaregia  was  then  diluted  with  distilled 
water  to  a  1:5  ratio.  The  extent  of  the  aquaregia  etching  of 
the  ITO  was  monitored  via  its  sheet  resistance  /?,  which 
increased  with  etching.  After  thorough  cleaning  with  an  ul¬ 
trasonic  surfactant  wash  and  rinsing  with  distilled  water  fol¬ 
lowed  by  acetone,  the  ITO  substrate  was  immersed  in  the 
diluted  aquaregia  in  an  ultrasonic  washer.  When  the  sub¬ 
strate’s  sheet  resistance  reached  the  desired  value,  it  was 
rinsed  in  distilled  water,  followed  by  acetone,  isopropanol, 
and  blown  dry  with  pure  N2. 

Both  the  HTL  and  Alq3  were  deposited  by  thermal 
evaporation  of  the  source  materials  in  a  vacuum  chamber 


with  a  base  pressure  of  10~6  Torr  located  in  a  high-purity 
dry  N2  glove  box,  in  which  water  and  oxygen  levels  were 
typically  less  than  1  ppm.  The  deposition  rate  of  the  organic 
layers  was  3-5  A/s.  A  1500  A  Mg0.9Ag0 A  metal  mixture  was 
coevaporated  on  the  Alq3  from  separate  sources  at  a  rate  of  5 
A/s  Mg  and  0.5  A/s  Ag,  in  a  separate  vacuum  chamber  with 
a  base  pressure  of  10-6  Torr.  A  300  A  thick  Ag  layer  was 
deposited  on  the  MgAg  layer  to  inhibit  oxidation.  The  broad¬ 
band  EL  intensity  IEL  was  measured  with  a  calibrated  silicon 
photodiode.  I(V)  and  /EL(V)  measurements  were  all  per¬ 
formed  in  air. 

Figures  3  and  4  display  I(V)  and  /E L(V),  respectively, 
of  LEDs  with  3-AS  and  TPD  deposited  on  treated  ITO  with 
varying  R.  To  minimize  the  effects  due  to  variations  in  the 
fabrication  procedure,  the  3-AS  and  TPD  devices  were  fab¬ 
ricated  simultaneously.  As  clearly  seen,  at  any  given  bias,  all 
devices  with  treated  ITO  had  moderately  higher  /  and  much 
higher  /EL  than  those  with  untreated  ITO.  Indeed,  relative  to 
devices  with  untreated  ITO,  optimal  treatment  enhanced 
/ EL  more  than  100-fold  and  1  more  than  tenfold  over  a  broad 
bias  range,  which  translates  into  a  considerably  higher  lumi¬ 
nescence  efficiency  77.  The  treated  ITO  devices  also  exhib¬ 
ited  a  lower  EL  threshold  voltage  Vih. 

We  note  that  the  optimal  treatment  depended  on  the 
HTL.  For  3-AS,  treated  ITO  with  R** 30  fi/sq  exhibited  the 
lowest  Vth(  —  4.2  V)  and  highest  /  and  /EL  at  ail  bias.  For 
TPD,  however,  the  devices  with  R^10  fi/sq  exhibited  the 
lowest  Vth(— 5  V)  and  highest  /  and  /EL  at  all  bias  voltages. 


FIG.  1 .  The  structure  of  a  typical  organic  LED  described  in  this  work. 


App!.  Phys.  Lett.  70  (20),  19  May  1997 


0003-6951/97/70(20)/2741/3/$10.00 


©  1997  American  Institute  of  Physics  2741 


FIG.  3.  /(V)  behaviors  of  LEDs  with  (a)  3-AS  as  the  HTL,  and  (b)  TPD  as 
the  HTL. 


some  aspects  of  the  I(V)  and  / EE(  V)  behaviors.10,11  It  is  also 
generally  accepted  that  tunneling  dominates  for  large  barriers 
while  SCL  current  injection  is  dominant  in  smaller  barriers. 
Based  on  the  latter,  Burrows  and  Forrest  suggested  that  in 
ITO/TPD/Alq/MgAg  devices  the  /( V)  curve  is  controlled  by 
trapped-charge-limited  (TCL)  electron  transport  through 


FIG.  2.  Molecular  structures  of  (a)  8-tris-hydroxyquinoline  A1  (Alq3), 
(b)  “three-armed  star”  (3-AS)  4,4\4"-tris[/V-(3-methoxyphenyl)- 
yV-phenylamine-triphenylamine],  and  (c)  triphenyl  diamine  (TPD). 

Note  that  the  ionization  potentials  of  3-AS  and  TPD  are  4.9 
and  5.2  eV,  respectively.  At  high  bias,  the  current  injected 
into  3-AS  was  much  higher  than  into  TPD,  but  /EL  of  the 
former  saturated  at  lower  V  than  that  of  the  latter.  In  the  high 
bias  region,  it  was  also  observed  that  the  decay  rate  of  /EL 
was  proportional  to  the  injection  current  density,  consistent 
with  other  studies.8 

Figures  3  and  4  also  suggest  that  in  LEDs  with  treated 
ITO  and  two  HTLs  such  as  3-AS/TPD,  lower  Vth,  higher 
/EL,  and  high  rj  should  be  observed  over  a  broad  bias  region, 
consistent  with  other  observations.9  We  note  that  in  all  de¬ 
vices,  the  I(V)  curves  followed  a  power-law  behavior  above 
Vth ,  but  below  Vth  they  clearly  deviated  from  such  behavior. 

In  order  to  understand  the  behavior  of  these  LEDs,  we 
performed  both  x-ray  photoelectron  spectroscopy  (XPS)  and 
atomic  force  microscopy  (AFM)  measurements  on  the 
aquaregia-treated  ITO.  XPS  indicated  no  chemical  difference 
between  the  treated  and  untreated  ITO.  However,  the  AFM 
images  presented  in  Fig.  5  show  that  the  surface  morphology 
of  the  treated  ITO  is  more  porous  than  that  of  the  untreated 
ITO.  In  OLEDs,  both  space-charge-limited  (SCL)  injection 
effects  and  tunneling  processes  have  been  used  to  explain 


Voltage  (V) 

FIG.  4.  /Fi  ( V)  behaviors  of  LEDs  with  (a)  3-AS  as  the  HTL,  and  (b)  TPD 
as  the  HTL. 
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FIG.  5.  AFM  images  of  (a)  untreated  ITO  (surface  sheet  resistance 
—  20  H/sq)  and  aquaregia-treated  ITO  with  surface  sheet  resistances  of  (b) 
-30  ft/sq,  (c)  -70  H/sq,  and  (d)  -90  H/sq. 


Alq3.  10  Yet,  the  porous  surface  morphology  of  aquaregia- 
treated  ITO  provides  more  contact  area  with  the  organic 
layer,  which  induces  greater  hole  injection  into  the  HTL. 
This  can  account  for  the  lower  Vth  in  the  LEDs  with  treated 
ITO;  the  power-law  behavior  of  I(V)  above  Vlh  is  consistent 
with  the  SCL  current  mechanism.  However,  the  observed 
dependence  of  I(V)  and  77  on  surface  morphology  of  the  ITO 
is  equivalent  to  their  dependence  on  the  ITO/HTL.  This  sug¬ 
gests  that  the  SCL  currents  are  not  only  controlled  by  TCL 
electron  transport  through  Alq3,  but  also  by  hole  injection 
from  the  ITO  to  the  HTL,  as  well  as  hole  transport  through 
the  HTL  and  injection  from  the  HTL  into  the  Alq3. 

SCL  effects  alone,  however,  cannot  completely  account 
for  the  relative  behavior  of  the  LEDs.  From  Fig.  5,  it  is  clear 
that  the  surface  of  the  —30  fl/sq  treated  ITO  is  not  the  most 
porous,  yet  the  LEDs  with  3-AS  fabricated  on  such  treated 
ITO  exhibited  the  largest  /  above  Vlh  among  all  devices  with 
either  3-AS  or  TPD.  These  results,  and  the  fact  that  the  volt¬ 
age  drop  across  the  ITO  (typically  —  1  mAX50  0^50  mV) 
is  negligible,  indicate  that  effects  other  than  the  ITO/HTL 
contact  area  are  significant.  Although  their  nature  is  unclear, 
our  results  do  show  that  with  the  optimized  aquaregia  treat¬ 
ment  enhanced  hole  injection,  EL  output,  and  EL  efficiency 
are  achieved. 

In  summary,  we  have  found  that  by  using  aquaregia- 
treated  ITO  substrates,  both  current  injection  and  EL  output 
of  double-layered  organic  LEDs  are  significantly  improved. 
The  improvement  is  believed  to  be  due  to  enhanced  hole 
injection  current  density  from  the  ITO  into  the  hole  trans¬ 
porting  layer,  resulting  from  the  increased  contact  area  as 
well  as  the  contact  conditions  between  them.  The  observed 
dependence  of  the  1(  V)  characteristics,  EL  output,  and  EL 
efficiency  on  the  ITO  surface  morphology  indicate  that 
space-charge-limited  currents,  which  are  determined  by  the 
bulk  properties  of  the  organic  materials,  are  dominant  in 
double-layered  Alq3-based  LEDs. 
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Abstract 

We  recently  gave  physical  arguments  that  the  recombination  in  a  single  type  of  dangling  bond  centers  wil1  yield  a  zero 
value  for  the  lk:it  intensity  exponent  of  the  minority  carriers,  A.  We  then  concluded  that  the  observed  A  *  0  vvue  in  a-Si:H 
indicates  the  presence  of  other  states  in  its  mobility  gap.  In  this  paper,  we  substantiate  those  arguments  by  an  analytic  proof, 
by  computer  simulations  and  by  considering  our  recent  experimental  findings  on  the  temperature  dependence  of  A.  In 
particular,  we  show  that  in  an  n-type  photoconductor  the  dangling  bond  acts  as  a  single  acceptor-like  recombination  center. 
©  1998  Elsevier  Science  B.V. 
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1.  Introduction 

The  ability  to  measure  the  ambipolar  diffusion 
length  in  hydrogenated  amorphous  silicon,  a-Si:H, 
has  provided  another  source  of  information  regarding 
the  state  distribution  in  this  material  [1-5].  However, 
one  parameter  which  can  be  derived  from  the  corre¬ 
sponding  measurements  has  hardly  been  discussed 
[1,2,5].  This  is  the  light  intensity  exponent,  A,  which 
is  associated  with  the  minority  carrier  recombination 
time,  rh,  by  the  definition,  Th  oc  G\  where  G  is  the 
carrier  generation  rate.  Recently,  we  have  applied  [1] 
a  Rose-like  [6]  argument  [7]  to  find  the  above-men¬ 
tioned  information  from  the  dependence  of  A  on 
various  controllable  parameters.  This  application 
yielded  the  conclusion  that  A  is  a  unique  parameter 
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in  the  sense  that  it  can  provide  a  qualitative  test  for 
the  presence  of  recombination  centers  other  than  a 
single  type  of  dangling  bonds.  In  particular,  we 
concluded  that  a  A  ^  0  is  an  indication  that  at  least 
one  other  type  of  recombination  center  must  be 
active.  Following  this  conclusion  we,  suggested  that 
in  a-Si:H,  one  must  assume  the  presence  of  other 
types  of  states.  Such  states  are  bandtails  [8-10] 
and/or  other  types  of  dangling  bonds  [2,4].  A  priori, 
the  well-known  [8-10]  existence  of  bandtails  would 
immediately  suggest  that  we  should  have  A  #  0  and 
thus  this  test  is  trivial.  However,  it  is  expected  [9,10] 
that  above  room  temperature,  the  recombination  rate 
in  the  bandtails  is  diminished  and  thus  the  A  #  0 
phenomenon  is  informative,  indicating  other  active 
recombination  processes.  This  makes  the  A  ^  0  test  a 
conclusive  tool  for  state  distribution  spectroscopy  for 
semiconductors  in  general  and  for  amorphous  semi¬ 
conductors  in  particular. 
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Following  these  considerations  and  since  we  gave 

[1]  only  qualitative  arguments  in  support  of  this  test, 
it  appeared  important  to  substantiate  the  above  argu¬ 
ment  by  an  analytical  proof  and  by  model  simula¬ 
tions.  This  need  was  further  amplified  when  the 
relevance  of  the  test  to  a-Si:H  was  questioned.  The 
doubts  came  following  an  analytical  argument  [5] 
that  there  is  a  possible  ‘curious  situation’  for  which  a 
A  =£  0  value  can  be  found  for  a  single  type  of  dan¬ 
gling  bonds  in  a-Si:H.  In  the  present  paper,  we  also 
disprove  that  claim  by  analytical,  numerical  and 
experimental  results. 

The  structure  of  this  paper  is  as  follows.  In  Sec¬ 
tion  2,  we  derive  analytically  the  expression  for  the 
lifetime  of  both  charge  carriers  in  the  case  of  a  single 
acceptor-like  level  showing  that  it  is  the  same  as  the 
one  obtained  for  the  single  type  of  dangling  bonds. 
Then  we  determine  the  corresponding  values  of  A.  In 
Section  3,  we  apply  the  model  to  a-Si:H  and  show 
that  it  cannot  produce  the  A  #  0  result  claimed  in 
Ref.  [5].  Finally,  in  Section  4,  we  discuss  experimen¬ 
tal  results  and  present  our  main  conclusion. 

2.  Analytical  derivation  of  the  light  intensity  expo¬ 
nents 

The  four  parameters  which  can  be  derived  from 
phototransport  measurements  are  the  electron  and 
hole  mobility-lifetime  products,  /xcrc,  /^hrh,  and  their 
corresponding  exponents,  y  and  A  [1-5].  The  light 
intensity  exponents,  y  and  A,  are  defined  by  re  a 
C7”1  and  r h  a  G\  We  note  in  passing  that  in  other 
papers  the  parameter  A  was  defined  as  2S  [1,3]  or  <£h 

[2] .  The  dependence  of  these  four  parameters  on 
£c  —  £F  has  been  used  [2,3,5]  for  deriving  informa¬ 
tion  on  the  state  distribution  in  a-Si:H.  Here,  £c  is 
the  conduction  band  edge,  £F  is  the  dark  (equi¬ 
librium)  Fermi  level. 

We  turn  then  to  prove  that  in  an  n-type  semicon¬ 
ductor,  the  single  dangling  bond  acts  as  a  single 
acceptor-like  recombination  center.  For  this,  we  as¬ 
sume  that  n(/p{  >  1,  where  n{  and  p{  are  the  steady 
state  concentrations  of  free  electrons  and  holes  under 
light  excitation.  Since  the  p-type  photoconductor 
yields  a  symmetric  situation,  we  do  not  discuss  it 
here. 

Let  us  start  our  analytic  derivation  with  the  sim¬ 


plest  version  of  the  Shockley-Read  analysis  for  a 
deeply  lying  acceptor-like  state  [6,7].  For  the  elec¬ 
trons  in  this  system  we  assume  a  capture  coefficient, 
(Tn°,  and  for  the  holes  a  capture  coefficient,  cr~.  If  the 
concentration  of  these  states  is  Na,  a  fraction  of 
them,  f~Na,  is  occupied  by  electrons  and  a  fraction 
of  them,  (1  jVa,  is  occupied  by  holes.  Under 
steady  state  conditions  the  recombination  rates  for 
electrons  and  holes  must  be  equal,  so  that: 
rtf(°n0^Xl  —  Using  the  defini¬ 

tion  of  the  ‘deep  trapping  times’  [7]  Tn°  s  l/(o*n°A^) 
and  t~  ®  \/(a~N^)  we  obtain  then  that  the  steady 
state  kinetics  is  governed  by  the  relation: 

«f0  -r)A,°=/>fr/v  (0 

Eq.  (1)  is  exactly  the  expression  obtained  in  Ref.  [5] 
for  the  n-type  photoconductor  (nf  >  pf)  in  which  the 
recombination  center  is  the  multilevel,  occupation- 
correlated,  dangling  bond.  In  particular,  it  is  impor¬ 
tant  to  note  (see  below)  that  this  result  is  found  in 
Ref.  [5]  for  both  the  high  and  low  nf/pf  cases.  This 
finding  and  the  fact  that  Eq.  (1)  was  derived  here 
using  the  single  level  Shockley-Read  model  justi¬ 
fies,  a  posteriori,  our  argument  [1]  that  the  dangling 
bond  acts  as  a  single  level-like  recombination  center. 
In  turn  it  disproves  the  claim  in  Ref.  [5]  that  the  low 
n(/p{  case  (the  ‘curious  situation’)  ‘does  not  exist  in 
the  classical  Shockley-Read  recombination  scheme’. 
Hence,  to  make  Eq.  (1)  pertinent  to  the  case  of  a 
single  type  of  dangling  bonds,  all  that  has  to  be  done 
is  simply  replace  Na  by  the  concentration  of  dan¬ 
gling  bonds,  Ndb. 

Let  us  now  discuss  y  and  A  which  are  relevant  to 
the  recombination  via  acceptor-like  states.  For  this 
purpose,  we  repeat  briefly  the  argument  given  by 
Rose  [6]  (for  y)  and  by  Balberg  [1]  and  Balberg  and 
Lubianiker  [3]  (for  A),  but,  this  time,  in  relation  to 
Eq.  (1). 

Consider  then  two  distinct  situations  under  the 
condition  (both  assumed  in  Refs.  [1,5])  that  Na  »  nf, 
p(.  The  first  case  is  when  f~  *  1  —  f  =*  0.5.  In  this 
case  [1],  the  increase  of  G  and  thus  the  increase  of 
n{  and  p(  will  hardly  effect  the  electron  occupation 
of  the  recombination  centers,  and  their  recombina¬ 
tion  times,  Tn°/(1  — /“)  and  r~/J^y  will  be  indepen¬ 
dent  of  G.  We  will  have  then  that  y  =  1  and  A  =  0. 
The  other  extreme  case  is  when  f~>  1  so  that 
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(1  — /*)Wa  >/if.  In  this  case,  the  charge  neutrality 
restriction  will  not  allow  n(  to  increase  with  G 
without  a  corresponding  increase  of  (1  —  so  that 
the  condition  nf  —  (1  —  f~ ) Afa  must  be  obeyed.  If  so, 
since  the  steady  state  recombination  rate  equals  G, 
i.e.,  (see  Eq.  (1))  Ca«f(l-f  )JVa,  we  get  that 
Ca/i;  and  thus  y=  1/2.  On  the  other  hand,  f~N* 
is  essentially  a  constant  and  since  G  a  pf/~jVa  we 
get  that  A  =  0.  The  transition  from  the  first  case  to 
the  other  can  be  obtained  [1]  by  moving  £F  towards 
Ec.  Consequently,  the  value  of  y  shifts  from  1  to 
1/2  with  increase  of  £F,  while  A  =  0  is  maintained 
throughout  the  increase  of  £F. 

In  principle,  it  is  possible  (see  Section  3)  that, 
even  if  nf  >  p(y  under  an  extreme  ratio  of  the  recom¬ 
bination  capture  coefficients,  we  will  still  have  a 
situation  where  f  1  so  that  the  charge  neutrality 
condition  is  pf=f~ Na- 

This  ‘curious  situation’  is  noted  in  Ref.  [5]  in 
which  (1  -/“)Na  is  practically  a  constant,  and  G  a 
pff~ Hence,  under  the  requirements  of  charge 
neutrality  we  get  that  G  a  p]  and  that  G  a  nf,  yield¬ 
ing  A  =  -1/2  and  y=  1.  It  is  important  to  note 
then,  that  in  general  we  will  never  have  a  situation 
where  y  =£  1  and  A  0  simultaneously,  since,  math¬ 
ematically,  we  cannot  have  that  both  f“  and  (1  -  f“) 
much  smaller  than  unity,  and,  physically,  since  the 
charge  neutrality  condition  can  affect  the  concentra¬ 
tion  of  only  one  of  the  charge  carriers. 


3.  Application  of  the  analytical  derivation  to  a- 
Si:H 

The  central  argument  of  Ref.  [5]  was  that  ‘in 
reality’  the  ‘curious  situation’  described  above  may 
account  for  the  observation  of  a  A  #  0  in  undoped 
a-Si:H.  We  will  show  now  that  this  is  very  unlikely 
for  a-Si:H,  by  using  the  very  recombination  parame¬ 
ters  suggested  in  Ref.  [5],  We  note  that  these  parame¬ 
ters  are  most  favorable  (within  the  range  of  known 
[8-10]  a-Si:H  parameters)  for  the  ‘curious  situation’ 
for  a-Si:H. 

In  principle,  the  general  ‘curious  situation’  may 
occur  under  the  special  combination  of  the  deep 
trapping  times  given  by  [5]: 


Tn°Ap°  «  n(/p(  rn°/Tp-  ,  (2) 

where  rp°«  is  the  deep  trapping  time  of  holes  in  the 
neutral  dangling  bonds.  The  right-hand  inequality  is 
already  the  case  implied  by  Eq.  (l)  for  the  ‘curious’ 
case  of  f «  1  -f.  As  for  the  left-hand  side  in¬ 
equality,  it  is  commonly  agreed  that  in  a-Si:H, 
Tn°/rp°  =  1.  Hence,  the  particular  way  of  writing  [5] 
unity  as  Tn°/rp0  does  not  reflect  the  multilevel  nature 
of  the  dangling  bond  but  rather  the  a  priori  assump¬ 
tion  that  we  have  an  n-type  photoconductor  (nf>  pf). 
Hence,  the  simple  single  acceptor  level-like  nature  of 
the  problem  is  maintained ,  in  the  ‘curious  situation’. 
Turning  to  the  numerical  evaluation  of  the  right  hand 
side  inequality  we  note  that  t°/t~  derived  from  the 
experimental  data  by  Ref.  [5]  (and  others  [8-10])  is 
Tn0/ t~  =  100.  This  means  (Eq.  (2))  that  the  ‘curious 
situation’  should  correspond  to  n(/p{ »  10.  Now,  for 
the  common  G  =  1019  cm"3  s'1  the  authors  of  Ref. 
[5]  find  that  the  maximum  value  of  pf  is  smaller 
than  4  X  10n  cm"3.  Accordingly,  n(  must  be  smaller 
than  1013  cm’3.  Since  Ndb  (as  given  in  Refs.  [5,8])  is 
larger  than  1016  cm'3,  we  obtain  that  jVdb/nf, 
Ndb/p{>  103. 

Returning  to  the  above  analysis  of  y  and  A  as 
given  in  Section  2,  we  remember,  that  in  order  for 
the  charge  neutrality  condition  to  affect  the  carrier 
concentration  (yielding  y  #  1  or  A  #  0)  we  must 
have  that  either  /if  -  (1  ~/’)A^db  or  /?f  ~/~ Ndb  (the 
‘curious  situation’).  In  view  of  those  requirements, 
we  conclude  from  the  above  discussion  that  for 
a-Si:H,  f~ /(l  —  f~ )  =  10~3.  Examining  Eq.  (1),  we 
see  that  even  if  n{  =  pf  (favorable  conditions  for  the 
‘curious  case’)  the  Tn0/r~  =  100  will  yield  that 
f~/(  1  10'2  rather  than  the  above  observed 

value  of  10'3.  We  saw  in  Section  2  that  as  long  as 
n(  <  (1  -  f')Ndb  or  p(  <f~NAb9  the  carrier  concen¬ 
trations  do  not  affect  the  occupation  of  the  recom¬ 
bination  centers  and  the  lifetime.  Hence,  even  under 
the  above  very  favorable  parameters  for  the  ‘curious 
case’,  we  get  that  y=  1  and  A  =  0  which  in  turn, 
explains  why  in  the  numerical  simulations  of  the 
single  type  of  dangling  bonds  model  [2]  one  always 
finds  the  A  =  0  result  for  all  reasonable  parameters 
of  a-Si:H.  The  interpretation  [5]  of  the  A  variation  as 
a  transition  between  the  ‘normal’  case  of  A  =  0  to 
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Ec-Ef  (eV) 

Fig.  1.  The  calculated  optically  excited  free  electron^and  free  hole 
concentrations  (a),  and  their  corresponding  light  intensity  expo¬ 
nents  (b),  as  a  function  of  the  Fermi  level  position,  for  the  single 
type  of  dangling  bonds  model.  Note  that  y  and  A  are  associated 
here  with  the  electrons  and  holes,  respectively.  For  the  ‘n-type’ 
material  (£C-£F^0.8  eV)  these  are  majority  and  minority 
carriers,  respectively. 

the  ‘curious  situation’  of  A  =  1/2  is  therefore  un¬ 
founded  in  a-Si:H,  and  one  must  resort  to  a  model 
that  has  at  least  one  other  type  of  recombination 
center  in  addition  to  the  single  level  dangling  bond 
centers. 

To  further  illustrate  the  above  point,  we  show  in 
Fig.  1  the  results  of  our  own  numerical  computations 
of  the  Ec  —  £F  dependence  of  nf,  pf ,  y  and  A,  when 
the  parameters  mentioned  above  were  utilized.  The 
other  parameters,  less  critical  for  the  present  work, 
were  taken  from  Refs.  [2,8]  and  the  numerical  calcu¬ 
lations  were  done  using  the  dangling  bond  statistics 
of  Ref.  [8].  We  show  these  results  in  order  to  affirm 
the  results  of  Morgado  [2]  as  well  as  (unlike  Mor- 
gado)  to  present  them  in  terms  of  nf  and  p(.  The 
‘n-type’  cases  discussed  above  are  given  in  Fig.  1  for 
the  interval  of  Ec  —  EF  ^  0.8  eV.  For  the  larger 
Ec  —  E?  values  the  material  is  ‘p-type’  and  the  be¬ 
havior  is  analogous  to  that  of  the  ‘n-type’  material.  It 
is  clearly  seen  that  for  the  ‘n-type’  range  the  decreas¬ 


ing  n(/pf  ratio  does  not  yield  the  ‘curious  situation’ 
(A  #  0)  claimed  in  Ref.  [5]  and  that  a  rather  simple 
n-type  to  p-type  transition  [1]  takes  place.  This  A  =  0 
finding  proves  then  that  the  ‘curious  situation’  is  not 
‘relevant’  to  a-Si:H  and  that  a  single  type  of  dan¬ 
gling  bonds  cannot  account  for  the  experimental 
A  ^  0  observation  in  undoped  a-Si:H. 


4.  Discussion  and  conclusions 

Following  the  results  of  Section  3,  it  is  obvious 
that  our  test  [1]  of  A  #  0  is  enough  to  establish  the 
existence  of  states  other  than  a  single  type  of  dan¬ 
gling  bonds,  but  it  is  not  enough  for  the  determina¬ 
tion  of  which  of  the  more  complex  models  is  the 
correct  one.  On  the  other  hand  the  dependence  of  A 
on  ‘external*  parameters,  such  as  Ec  -  E?  or  the 
temperature,  can  yield  some  new  conclusions  regard¬ 
ing  the  state  distribution  in  a-Si:H.  For  example,  we 
expect  [9,10]  that  the  role  of  the  recombination  in  the 
bandtails  diminishes  with  increasing  temperature. 
According  to  our  A  0  test,  we  expect  then  that  in 
this  material  we  will  find  a  A  ->  0  behavior  with 
increasing  T.  Indeed,  we  have  experimentally  found 
[11]  this  expected  behavior.  This  finding  yields  then 
further  support  to  the  applicability  of  our  test  for 
a-Si:H.  It  also  shows  that  this  test  can  yield  informa¬ 
tion  that  is  otherwise  obtained  by  comparison  of  the 
experimental  results  with  detailed  numerical  compu¬ 
tations  [9,10]. 
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OPTICAL  PROPERTIES  AND  NANOPARTICLE  SIZE  DISTRIBUTION  IN 
LUMINESCENT  POROUS  SILICON 
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Physics  Dept,  University  of  Puerto  Rico,  P.O.  Box  23343,  San  Juan,  PR,  USA  00931 
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We  measured  the  photoluminescence,  optical  transmission  and  Raman  spectra  of  free 
standing,  luminescent  porous  Si  films  of  higji  optical  quality.  The  optical  transmission  data 
were  used  to  derive  the  size  distribution  of  the  nanoparticles  in  the  porous  Si  films.  The 
procedure  is  based  on  the  assumption  that  interband  transition  is  the  principal  mechanism  for 
the  absorption  and  that  quantum  confinement  in  the  nanoparticles  determines  the  optical  gap. 

The  distribution  obtained  is  centered  around  22  Angstroms  and  is  not  symmetric,  exhibiting  a 
tail  towards  larger  sizes.  On  the  basis  of  this  distribution  we  predicted  a  shift  of  ~20  cm'1  in 
the  521  cm'1  silicon  line  in  the  Raman  spectrum. This  shift  was  indeed  observed  in  our 
samples,  thus  confirming  the  size  distribution  obtained  from  the  optical  data. 


1  Introduction 

It  is  generally  accepted  that  quantum  confinement  plays  a  crucial  role  in  the 
photoluminescence  properties  of  porous  silicon  (PSi).  The  structure  of  the  porous 
system  depends  on  the  preparation  method  but,  in  general,  nanoscale-sized  silicon 
crystallites  are  required  to  get  visible  luminescence.  Extensive  work  has  been  done 
during  the  last  seven  years  to  characterize  the  PSi  nanostructure  and  to  correlate  it 
with  its  optical,  photoluminescence  and  electroluminescence  properties.  Raman 
analysis  is  an  important  tool  to  characterize  Si  nanostructures  giving  information 
about  size  and  shape  [1].  Raman  studies  of  electrochemically  etched  PSi  report 
spherical  nanostructures  in  most  of  the  cases  [2-4].  Many  of  the  size  analyses  using 
the  Raman  technique  yield  only  the  average  size  of  the  particles.  The  particles’  size 
distribution,  F(r),  is  addressed  only  in  an  indirect  way  through  the  Fourier 
coefficients  of  the  confined  phonon  wavefunctions  [5]. 

In  order  to  study  the  nanostructure  characteristics  of  high  porosity  PSi  samples, 
we  measured  two  independent  quantities,  the  optical  transmission  and  the  Raman 
spectrum  High  optical  quality,  free-standing  luminescent  porous  silicon  was  used 
We  derived  F(r)  from  the  optical  transmission  data  and  calculated  the  Raman 
spectrum  associated  with  this  distribution  of  particles.  The  calculated  spectrum  was 
compared  with  the  measured  Raman  data. 
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Figure  1.  Transmission  spectrum  of  the  PSi  film.  Figure  2.  Raman  spectra  of  the  Si  wafer  and  the  PSi  film. 


2  Experiments  and  Calculations 

PSi  samples  were  prepared  by  electrochemical  anodization  of  a  p-type,  0.5-1. 5 
Qcm,  (100)  silicon  wafer  using  a  HF:ethanol:water  (1:2:1)  solution.  The 
anodization  current  (200  mA/cm2)  was  applied  in  280  steps  of  0.8  sec  each  and  with 
resting  intervals  of  1.4  sec  between  anodizations,  in  order  to  ensure  a  high-optical 
quality  film. 

The  optical  transmission  spectrum  was  measured  by  a  Perkin  Elmer  Lambda  6 
spectrometer.  The  results  are  shown  in  Fig.  1.  The  Raman  spectra  of  the  crystalline 
silicon  wafer  and  of  the  free  standing  sample  were  taken  by  a  ISA  T64000 
spectrometer.  The  results  are  plotted  in  Fig.  2.  The  peak  of  the  spectrum  of  the 
crystalline  Si  is  obtained  at  521  cm"1  as  expected.  The  peak  for  the  PSi  sample 
occurs  at  501  cm*1,  a  shift  of  20  cm"1.  The  optical  transmission  spectrum  of  the 
sample  was  used  as  the  starting  point  to  calculate  the  size  distribution  of  the  silicon 
particles  in  the  PSi  sample.  This  calculation  was  done  by  assuming  that  the  porous 
silicon  sample  is  a  collection  of  spherical  crystallites  whose  sizes  are  described  by  a 
distribution  function,  F(r),  and  also  that  the  absorption  is  due  to  the  interband 
transition  of  the  electrons  in  the  individual  crystallites. 

We  assume  that  each  particle  is  transparent  for  energies  (E)  lesser  than  its 
optical  gap  (Eg)  and  its  absorbance  increases  abruptly  for  E>Eg,  the  transmission 
coefficient  of  a  film  can  be  expressed  as, 

00 

T(X)  oc  1  -  C  |  F(r)A(r,  X)dr 
0 
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Figure  3.  Calculated  F(r)  and  log-noimal  fit  Figure  4.  Comparison  between  the  observed  Raman 

spectrum  (dots)  and  the  predicted  one  (solid  line) 


whe  re  C  is  a  constant  and  A(r.),)  is  proportional  to  the  energy  absorbed  by  a  particle 
of  radius  r.  A(rA)  is  approximated  by  the  Heavyside  step  function, 
A(r,X)  cc  ®(A.gr  - X)  where  XgT  =  2nhc/EgT,  and  EgJ  is  the  energy  gap  for 

particles  of  size  r.  From  the  above  expression,  we  obtain  that  the  first  derivative  of 
the  transmission  spectrum  is  proportional  to  the  size  distribution  density, 

00 

oc  -  J  F(r)8(r  -r0)drcc  -F(ro ) 

Xq  0 

where  tq  is  the  size  which  gives  an  energy  gap  corresponding  to  >-q. 

Figure  3  shows  the  size  distribution  (solid  line)  obtained  from  the  transmission 
spectrum  (Fig.1)  under  the  above  considerations  giving  a  non-symmetric  size 
distribution  peaking  at  2.2  nm.  To  obtain  this  size  distribution,  we  assumed  a 
spherical  particle  shape  [2-4]  and  used  the  relation  between  the  magnitude  of  the 
energy  gap  and  the  size  of  the  (spherical)  particle  suggested  in  reference  [6], 
—  1  3 

Eg,  r  =1-1  +  3  .65r  *  where  Eg^r  is  in  eV  and  r  in  nm. 

Using  the  size  distribution  given  in  Fig.3?  the  Raman  spectrum  for  the  sample 
was  calculated  as 

00 

I  predicted  (®)  =  1  F(r)I©  (r)dr 
0 

where  Iffl(r)  is  obtained  from  the  phonon-confinement  model  [1,  4]. 

Figure  4  shows  a  comparison  between  the  calculated  Raman  intensity  and  the 
experimental  data.  We  notice  that  the  calculation  reproduces  accurately  the  peak 
position  but  is  somewhat  wider  than  the  experimental  curve. 


cT 

dX 
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3  Conclusions 

Two  independent  measurements,  optical  transmission  and  Raman  spectrum,  were 
used  to  determine  the  size  distribution,  F(r),  in  highly  porous  luminescent  PSi  films. 
F(r)  obtained  from  the  transmission  spectrum  peaks  at  2.2  nm  and  agrees  with  the 
Raman  observation  assuming  that  the  PSi  film  is  composed  of  spherical  particles 
rather  than  wires.  A  wire-like  nanostructures  would  yield  a  size  distribution  peaking 
at  smaller  sizes  and  thus  would  predict  a  Raman  shift  larger  than  the  observed  20 
cm'1.  This  provides  support  to  the  assumption  that  PSi  films  of  high  porosity  are 
composed  of  particles  nearly  spherical  in  shape.  The  width  of  the  Raman  peak,  as 
predicted  from  F(r),  is  larger  than  the  observed  one.  Probably  because  the  size 
distribution  was  obtained  assuming  an  abrupt  increase  of  the  optical  absorption  at 
the  gap  energy  of  each  nanoparticle.  This  simplification  can  be  improved  by  the 
addition  of  a  parabolic  dependence  on  energy  of  the  absorption  coefficient  of  the 
nanoparticles.  Such  a  dependence  would  reduce  the  width  of  F(r). 
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imnoft  Jnt  if the  no";exponentiaI  luminescence  decay  and  extraction  of  additional  physical  information  from  it  are 
important  issues  for  the  spectroscopic  study  of  porous  Si.  In  this  work  a  new  approach  to  extract  the  decay  time 
distribution  function  is  proposed,  based  on  the  fact  that  the  non-exponential  luminescence  decay  of  porous  Si  is 
a  multiple  exponential.  It  was  found  that  the  integrated  intensities  of  different  exponential  components  could  be  fitted  by 
a  Gaussian  distribution  function,  characterized  by  t0,  the  decay  time  with  maximum  probability,  and  d,  the  width  of  the 
n  funct,on-  Thfe  Parameters  i0  and  d  were  found  to  vary  with  emission  wavelengths.  t0  decreases 

potentially  with  an  increase  of  the  emission  energy  and  d  has  a  smaller  value  in  the  middle  of  the  emission 
wavelengths.  ©  1999  Elsevier  Science  B.V.  All  rights  reserved. 


Keywords:  Porous  Si;  S-band  luminescence  decay 


1*  Introduction 

Porous  Si  has  become  increasingly  attractive 
since  its  efficient  visible  luminescence  at  room  tem¬ 
perature  was  reported  in  1990  [1,22].  The  S-band 
emission  (400-800  nm,  slow  decay)  has  been 
studied  extensively  because  of  its  great  potential  for 
applications  [1,22].  Much  attention  has  been  paid 
to  the  study  of  the  temporal  behavior  of  the 
luminescence,  since  it  is  important  for  understand- 
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ing  the  nature  of  the  luminescence  mechanism.  The 
luminescence  decay  of  porous  Si  was  found  to  ex¬ 
hibit  a  complicated  non-exponential  behavior.  The 
characterization  of  the  luminescence  decay  and  ex¬ 
traction  of  physical  information  from  it  have  been 
discussed  extensively  in  the  literature.  Xie  has  de¬ 
fined  the  decay  time  as  the  time  when  the  lumines¬ 
cence  intensity  has  fallen  by  a  factor  of  e  [2]. 
However,  the  fast  components  characterized  by  this 
decay  time  may  not  dominate  the  multi-exponen¬ 
tial  decay.  Vial  has  defined  an  “average  lifetime”  as 
the  integral  of  the  intensity  with  respect  to  time, 
divided  by  the  initial  luminescence  intensity  [3]. 
But  a  single  parameter  cannot  fully  characterize 
a  non-exponential  decay  behavior.  A  number  of 
workers  have  analyzed  the  form  of  the  lumines¬ 
cence  decay  by  fitting  a  stretched  exponential  or  by 
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extracting  lifetime  distributions  from  the  multi-ex¬ 
ponential  decays  [4-10].  However,  a  deviation  of 
the  decay  curves  from  the  stretched  exponential 
was  observed  in  the  longer  time  region  at  low 
temperature  [9].  Some  of  the  proposed  lifetime 
distribution  functions  seem  difficult  to  understand 
with  clear  physical  pictures. 

In  this  work,  a  new  method  to  derive  the  decay 
time  distribution  function  is  introduced  to  analyze 
the  luminescence  data  of  our  porous  Si  samples, 
assuming  a  decay  time  distribution  originating 
from  isolated  (or  confined)  multiple  optical  centers 
in  a  random  system.  It  is  found  that  the  integrated 
intensities  of  different  exponential  components  can 
be  fitted  by  a  Gaussian  distribution  function,  which 
was  characterized  by  r0,  the  decay  time  of  the 
component  with  maximum  probability,  and  d,  the 
width  of  the  Gaussian  function.  r0  and  d  for  differ¬ 
ent  emission  wavelengths  were  obtained  through 
the  fitting  procedure.  t0  was  found  to  decrease 
exponentially  with  an  increase  of  the  emission  en¬ 
ergy  and  the  smallest  d  (the  narrowest  distribution) 
was  found  to  be  at  the  middle  wavelength  region  of 
the  emission. 


2.  Samples  and  experiments 

The  samples  were  prepared  by  the  conventional 
electrochemical  anodization  method  [11,12,23]. 
The  electrolytic  solution  for  the  anodic  etching  was 
1:2:1  (HF :  ethanol :  water).  The  cell  included 
a  standard  Platinum  electrode  as  the  cathode  and 
a  Boron-doped  crystalline  silicon  wafer  with  a  re¬ 
sistivity  in  the  range  0.5-1.5  SI  cm  as  the  anode.  The 
aluminum  ohmic  contact  at  the  bottom  of  the  wafer 
was  made  by  sputtering  deposition  and  annealing. 
Previous  to  the  anodization,  the  silicon  surface  of 
the  wafer  was  cleaned  with  20%  HF-aqueous  solu¬ 
tion  for  5  min.  Samples  were  anodized  with  a  cur¬ 
rent  density  of  100mA/cm2  for  600  s  and  finally 
dried  using  the  critical  point  method  [13]. 

A  Nd :  YAG  laser  with  a  pulse  width  of  35  ps  and 
an  optical  parametric  generator  were  employed  as 
a  tunable  pulse  excitation  source.  A  double  Spex 
1403  spectrometer  was  used  to  analyze  lumines¬ 
cence  signals.  Time-resolved  emission  spectra  were 
measured  using  a  SR250  boxcar  averager  with  time 


resolution  of  0.05  ps.  Luminescence  decay  was 
measured  by  a  SR430  multi-channel  transient 
scaler.  The  temperature  of  the  samples  was  control¬ 
led  by  a  helium-close-cycle  refrigerator  from  10  K 
to  room  temperature. 


3.  Experimental  results  and  discussion 

Luminescence  decay  curves  of  different 
wavelengths  at  11  K  and  time-resolved  emission 
spectra  at  room  temperature  were  measured  under 
excitation  with  532  nm  laser  pulses.  Typical  results 
are  shown  in  Figs.  1  and  2,  respectively,  and  are 
similar  to  those  reported  in  literature  [2-10].  The 
decay  curves  are  not  exponential.  In  Table  1,  the 
values  of  t0  and  re  are  listed  for  different  emission 
wavelengths.  xe  is  the  time  when  the  fluorescence 
intensity  falls  to  \je  of  the  initial  value.  t0  is  the 
decay  time  of  the  component  with  maximum  prob¬ 
ability,  derived  using  the  new  approach  introduced 


0  50  100  150  200  250  300  350 

Decay  time  ( ps) 

Fig.  1.  Luminescence  decay  of  different  wavelengths  under  exci¬ 
tation  with  532  nm  at  11  K.  (1)  700  nm;  (2)  640  nm;  (3)  600  nm; 
(4)  540  nm. 
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in  this  work.  The  luminescence  decay  becomes 
slower  at  longer  emission  wavelengths.  t0  increases 
from  30.8  ns  at  540  nm  to  432.6  jj.s  at  750  nm. 

It  is  generally  believed  that  the  blue-shifted  vis¬ 
ible  emission  from  porous  Si  samples  is  the  result  of 
quantum  confinement  [1-3,14,15].  Structure  stud¬ 
ies  revealed  that  porous  Si  has  a  wire-like  micro¬ 
structure  with  wire-width  randomly  distributed 


Fig,  2.  Time-resolved  emission  spectra  of  a  porous  Si  sample  at 
room  temperature,  excited  at  532  nm.  (1)  4  ps;  (2)  14  ps;  (3)  24  ps; 
(4)  44  ps;  (5)  64  ps;  (6)  84  ps. 


[14,15].  Electrons  participating  in  the  visible 
luminescence  are  those  excited  in  those  highly  con¬ 
fined  regions. 

For  the  time-resolved  emission  spectra,  the  peak , 
position  shifts  from  640  to  690  nm  with  an  increase 
of  the  delay  time.  The  band  width  was  found  to 
have  no  significant  broadening  with  increasing 
temperature  from  10  to  300  K.  All  the  results  imply 
that  the  emission  spectra  are  inhomogeneously 
broadened  and  are  an  overlap  of  the  emission 
bands  from  different  confined  regions,  each  having 
an  exponential  decay  with  its  own  characteristic 
decay  time. 

The  questions  raised  are  how  the  exponential 
components  distribute  (intensity  versus  lifetime) 
and  which  component  dominates. 

It  is  assumed  that,  due  to  strong  quantum  con¬ 
finement,  the  excited  electron-hole  pairs  are  localiz¬ 
ed  and  spatial  energy  transfer  is  negligible.  The 
emission  of  each  confined  region  has  an  exponen¬ 
tial  decay.  The  measured  luminescence  decay  at 
certain  wavelength  can  then  be  expressed  as 

/  =  £  /,  exp(  - 1 f  *()  =  £  — exP(  ~  (1) 

i  J  Zi 

where  Qt  =  /fTi,  it  is  the  integrated  intensity  of  ith 
exponential  component.  After  careful  study  of  the 
decomposition  of  decay  curves  of  different  emission 
wavelengths,  it  was  found  that  Q,  =  /jtf  can  be 
fitted  by  a  Gaussian  function  of  In  t(. 

Qi(?i)  =  A  •  exp{  -  4  In  2((ln  tf  -  In  t0)/d)2),  (2) 

where  t0  and  d  are  the  fitting  parameters,  t0  is  the 
decay  time  of  the  exponential  component  which 
has  maximum  distribution  probability,  d  is  the 
full-width  at  the  half-height  and  A  is  a  constant 


Table  1 

Fitting  parameters  t0  and  d  of  a  sample  for  different  emission  wavelengths.  Excitation  at  532  nm  at  11  K. 


X  (nm)  540  565  600  630  660  720  750 


to(Hs) 

d 


30.8 

48.4 

80.0 

117.7 

162.7 

278.0 

432.6 

±0.8 

±1.2 

±2.0 

±2.9 

±4.1 

±7.0 

±  10.8 

2.78 

2.61 

2.41 

2.15 

2.09 

2.60 

2.96 

±0.07 

±0.07 

±0.06 

±0.05 

±0.05 

±  0.07 

±  0.07 

9.4 

18.2 

33.6 

58.2 

82.2 

102.5 

116.7 

±0.1 

±0.2 

±0.3 

±0.6 

±0.8 

±1.0 

±1.2 

t.(MS) 
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determined  by  the  initial  intensity  of  the  fitted 
decay  curve. 

The  Gaussian  distribution  function  (2)  was  used 
to  fit  the  decay  curves  of  different  emission 
wavelengths  at  different  temperatures.  All  decay 
curves  were  fitted  quite  well. 

The  following  equation  was  used  for  routine  fit¬ 
ting  procedure: 

xexp(-i) 

“'4J.(i“P(-4,n2 

where  t(  =  n1  x0.  The  parameter  A  can  be  deter¬ 
mined  by  the  initial  intensity  /(0) 

A  =  /(0)ro j  £  n~' exp 4 In 2 j-  (4) 

Therefore, 


/(r)  =  7(0)  [  f  n-Mn(-41n2^p)2) 

xe*p(-41„2(ii^y)]. 


(5) 


The  number  of  terms  in  the  sum  is  2m  -I- 1.  m,  as  an 
appropriate  integer  number,  is  chosen  to  be  3  or 
larger.  The  ratio  n  will  be  automatically  optimized 
during  the  fitting  procedure.  Two  fitting  para¬ 
meters  t0  and  d  in  Eq.  (5)  can  be  easily  obtained  by 
a  simple  fitting  program.  The  fitting  parameters  are 
not  sensitive  to  the  initialization  values.  There  are 
a  unique  set  of  fitting  parameters  for  a  given  decay 
curve.  User’s  function  in  Microcal  Origin  was  used 
in  this  work  as  a  routine  fitting  procedure  to  ana¬ 
lyze  the  decay  curves  of  the  porous  Si  samples.  As 
an  example,  the  fitting  of  the  decay  curve  and  the 
distribution  Qt  versus  In  x,-  with  m  -  3  for  the 
660  nm  emission  at  290  K  are  shown  in  Fig.  3a,  b, 


(a)  Decay  time  (ns) 


i  =  0 


Fig.  3.  (a)  A  decay  curve  fitted  by  a  multi-exponential  function, 
(b)  Qi  versus  In  r  of  the  fitting  function. 
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respectively.  Only  a  small  number  of  the  experi¬ 
mental  points  in  Fig.  3a  were  used  for  clarity. 

The  suitability  of  the  fitting  function  was  checked 
using  data  published  by  other  authors.  For 
example,  Ookubo  found  that  decay  curves  of  their 
samples  fit  a  stretched  exponential  function  at  tem¬ 
peratures  above  20  K,  and  deviated  in  the  long  time 
region  at  20  K  [9].  We  reproduced  some  of  their 
data  from  the  parameters  in  the  paper  and  fit  decay 
curves  for  T  ^  20  K  satisfactorily  in  whole  time 
region  with  our  fitting  function. 

I  The  fitting  parameters  t0  and  d  of  a  sample  for 
different  emission  wavelengths  at  11  K  are  listed  in 
Table  1.  Values  of  xt  were  included  in  the  table  for 
comparison.  The  decay  time  t0  increases  with  an 
increase  of  the  emission  wavelength  and  varies 
from  30.75  ps  for  540  nm  to  432.6'  ps  for  750  nm. 

I  The  values  of  r0  and  t,  are  quite  different.  The 

1  width  d  of  the  Gaussian  distribution  is  smaller  near 
the  center  of  the  emission  band  and  larger  at  the 
extreme  wavelengths.  For  660  nm,  d  =  2.09,  and 
the  ratio  of  the  decay  time  of  the  half  height  at  the 
longer  side  to  that  of  the  shorter  side  was  found  to 
be  about  8,  while  at  750  nm  the  width  d  is  2.96  and 
the  ratio  is  about  19.  The  decay  time  distribution 
functions  Q  against  In  t  for  different  emission 
wavelengths  are  shown  in  Fig.  4a.  Redrawing  the 
distribution  functions  Q  in  the  abscissa  of  t  as 
shown  in  Fig.  4b,  it  is  clearly  seen  that  Q  has  a  long 
tail  in  the  longer  decay  side.  This  is  similar  to  what 
was  expected  for  disordered  systems  as  discussed 
by  Scher  and  others  [16]. 

It  is  interesting  to  note  that  r0  and  the  emission 
photon  energy  w  satisfy  the  following  exponential 
relation: 

To  =  tc  exp(  -  w/we),  (6) 

where  tc  and  wc  are  constants.  A  typical  result  for 
the  emission  of  a  sample  excited  at  532  nm  and 
11  K  is  shown  in  Fig.  5  with  the  fitting  parameters 
tc  =  2.77  x  10s  ps  and  wc  =  0.25  eV.  te  is  also  plot¬ 
ted  in  the  figure  for  comparison. 

The  decrease  in  the  decay  time  with  increasing 
emission  energy  has  been  ascribed  to  the  increases 
of  the  oscillator  strength  of  the  radiative  transition 
at  the  high-energy  side  resulting  from  stronger 
quantum  confinement  [17-19].  It  was  found  that 
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Fig.  4.  The  Gaussian  distribution  function  of  decay  times  for 
different  emission  wavelengths  excited  at  532  nm  at  1 1  K,  plot¬ 
ted  in  (a)  vs.  In  t  and  in  (b)  versuss  t  (1)  540  nm;  (2)  600  nm;  (3) 
660  nm;  (4)  720  nm;  (5)  750  nm. 


the  Si  conduction  band  edge  shifted  to  higher  en¬ 
ergy  and  the  oscillator  strength  of  the  radiative 
transitions  increased  with  a  decrease  of  particle 
sizes  [20,21].  The  variation  of  the  decay  time  with 
temperature  in  the  high-temperature  region  has 
been  reported  by  Vial  et  al.  [3].  They  found  that  the 
decay  rate  increased  exponentially  with  the 
luminescence  energy  up-shift.  This  observation  is 
consistent  with  our  result  described  by  Eq.  (6). 
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Emission  photon  energy  w  (ev) 

Fig.  5.  In  t0  versus  emission  enerev  w.  nnrW 


532  nm  at  1J  K. 


energy  w,  under  excitation  at 


4.  Summary 
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decay  of  fluorescence  of  porous  Si  is  proposed  It 
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xponential  components  could  be  fit  by  Gaussian 
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d,  the  width  of  the  Gaussian  function.  The  decay 
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ABSTRACT 

In  a  recent  letter  we  have  reported  the  first  observation  of  the  phenomenon  of  minority 
carrier-lifetime  sensitization  in  hydrogenated  amorphous  silicon  (a-Si:H).  We  find  now  that 
combining  the  study  of  this  phenomenon  with  the  study  of  the  well-known  phenomenon  of 
majority  carrier  lifetime  sensitization,  in  this  material,  can  provide  direct  information  on  its 
density  of  states  (DOS)  distribution.  This  finding  is  important  in  view  of  the  limitations 
associated  with  other  methods  designed  for  the  same  purpose.  We  have  carried  out  then  an 
experimental  study  of  the  effect  of  light  soaking  on  the  phototransport  in  a-Si.H.  We  found  that 
the  increase  of  the  dangling  bond  concentration  with  light  soaking  affects  the  sensitization  and 
thermal  quenching  of  the  majority  carriers  lifetime.  Using  computer  simulations,  we  further 
show  that  the  details  of  the  observations  associated  with  the  sensitization  effect  yield 
semiquantitative  information  on  the  concentration  and  character  of  the  recombination  centers  in 
a-Si:H. 


INTRODUCTION 

The  phenomenon  of  thermal  quenching  of  the  photoconductivity  in  intrinsic  a-Si:H  is  known 
for  many  years  [1].  The  interpretation  of  this  phenomenon  is  two-fold.  First,  it  requires  that  the 
capture  coefficients  of  the  electrons  is  much  smaller  for  the  bandtail  states  (BTS)  than  for  the 
dangling  bonds  (DB),  and  second,  it  requires  that  with  increasing  temperature  (T)  the  number  of 
BTS  available  for  recombination,  Pv,  decreases  until  at  some  T  it  equals  the  concentration  of  the 
dangling  bonds  ND.  While  this  basic  picture  is  in  accord  with  the  results  of  numerical  simulations 
[2,3],  we  do  not  know  of  a  direct  experimental  proof  that  establishes  this  model  and  that  ties  this 
phenomenon  with  the  phenomenon  of  sensitization.  The  latter  phenomenon  is  that  the  addition  of 
excess  recombination  centers  can  prolong  the  lifetime  of  the  charge  carriers  in  a  given 
photoconductor  [4],  In  this  paper  we  study  this  effect  by  applying  the  tool  of  light  soaking  (LS) 
in  order  to  change  the  concentration  of  desensitizing  recombination  centers  in  a-Si:H.  We  do  that 
here  in  the  “reverse  mode”  to  the  one  used  classically  [4],  i.e.  by  keeping  the  total  concentration 
of  the  sensitizing  centers  nearly  constant  and  changing  the  concentration  of  the  desensitizing 
recombination  centers. 

After  establishing  experimentally  the  above  basic  picture  of  the  sensitization,  we  have  turned 
to  numerical  simulations  that  enabled  us  to  evaluate  the  relation  between  the  manifestations  of 
the  phenomenon  of  the  sensitization  and  various  parameters  of  the  concentration  and  capture 
coefficients  of  the  recombination  centers  in  the  mobility  gap  of  a-Si:H.  As  we  show  below,  these 
features  give  further  semiquantitative  spectroscopic  information  on  the  concentration  and 
character  of  the  recombination  centers  in  the  material. 
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EXPERIMENTAL  STUDIES 


The  samples  used  in  the  present  experimental  study  were  device  quality  a-Si:H  films  which 
were  deposited  [5]  on  a  Coming  7059  glass  and  had  a  thickness  of  1.3pm.  The  temperature  of 
the  substrate  at  deposition  was  250°C.  Two  planar  contacts  of  NiCr  were  predeposited  on  the 
substrate  in  a  configuration  that  we  have  described  previously  [6].  The  photoconductivity  oPh 
was  measured  as  a  function  of  the  (HeNe  laser)  illumination  intensity  and  T,  in  the  T  range  of  80 
<  T  ^  300K.  The  latter  yielded  a  carrier  generation  rate,  G,  of  5xl018  £  G  <  5xl020cm'V l.  The 
value  of  the  mobility-lifetime  product  (pt)e  of  the  electrons  (the  majority  carriers)  was  derived 
from  the  relation  (pT)e  =  oPh/qG  where  q  is  the  electronic  charge.  It  was  thus  implicitly  assumed 
that  the  carrier  generation  quantum  efficiency  is  unity  and  the  absorption  of  the  light  is  uniform. 
To  find  the  light  intensity  exponent,  ye,  at  a  given  T,  the  best  fit  of  the  data  to  a  linear  relation 
between  log  aPh  and  log  G  has  been  used.  The  importance  of  the  value  of  yc  is,  as  will  be 
discussed  below,  that  it  reflects  semiquantitatively  some  basic  features  of  the  recombination 
centers  [4].  In  particular,  when  ye  >  1,  the  T  dependence  of  ye  yields  information  on  their  energy 
position,  their  concentration  and  their  capture  coefficient.  In  the  present  study  we  have  measured 
the  T  dependence  of  (px)e  and  ye,  before  and  after  LS  of  the  a-Si:H  film.  This  was  done  in  order 
to  find  out  what  spectroscopic  information  can  be  derived  from  these  dependences.  In  the  present 
study  the  light  soaking  was  carried  out  under  a  HeNe  illumination  of  100  mW/cm2  for  10  hours. 

The  results  of  our  measurements  are  presented  in  Fig.  1.  In  this  figure  we  show  the  T 
dependence  of  (px)e  before  and  after  LS.  It  is  clearly  seen  that  before  LS  there  is  a  clear  change 
in  the  sign  of  the  slope  of  the  aph(T)  dependence  around  150°K.  We  call  the  temperature  of  this 
change  the  T  of  thermal  quenching, 

Tu,.  After  LS  the  (px)e  drops,  as 
expected,  and  the  To,  shifts  to  a 
somewhat  lower  T.  This  is  in 
agreement  with  previous 
experimental  results  and  numerical 
simulations  [2]  of  the  effect  of  the 
increase  of  the  DB  concentration  > 
on  the  (px)c  dependence  on  T.  The  *| 
sensitization,  as  mentioned  above,  g, 
is  manifested  by  a  ye  >  1  value  [4]. 

In  the  present  case  we  find  that,  by 
comparison  of  the  above  behavior  3 
with  the  ye(T)  dependence  given  in 
the  inset  of  Fig.  1,  the  sensitization 
and  the  thermal  quenching  are  two 
aspects  of  the  same  phenomenon. 

This  is  since  we  see  that  there  is  a 
clear  downward  shift  in  the  yc(T) 
peak  as  well  as  a  decrease  in  its 
amplitude  around  Ttq.with  the  Temperature  (K) 

increase  of  LS 

The  above  features  of  the  (px)c  Figure  1.  The  measured  T  dependence  of  the  mobility- 

and  y^  dependences  on  T  are  simply  lifetime  product  and  the  corresponding  light  intensity 

explained  along  the  lines  of  the  exponent  for  atypical  sample  of  intrinsic  a-Si:H. 
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general  Rose  model  [4]  as  follows.  The  fact  that  the  (jix)e  values,  before  LS,  are  much  closer  to 
each  other  at  low  T  indicates  that  the  recombination  is  not  dominated  there  by  the  recombination 
centers  that  are  induced  by  the  LS.  This  is  in  contrast  with  the  situation  at  the  higher  T. 
Following  the  well-known  observation  that  LS  increases  the  concentration  of  the  dangling  bonds 
Np  but  does  not  cause  an  increase  in  the  concentration  of  the  valence  band  tail  (VBT)  states  [7,8] 
it  is  expected  that  the  latter  states  are  responsible  for  the  recombination  at  low  temperatuares.  It 
thus  appears  that  at  the  higher  T  the  recombination  is  controlled  by  the  dangling  bonds  (DB), 
while  at  the  lower  T  it  is  controlled  by  the  VBT  states.  We  consider  the  VBT  states  rather  than 
the  conduction  band  tail  (CBT)  states  since  the  lowering  of  T_  causes  the  shift  of  the  quasi-Fermi 
level  for  holes  towards  the  valence  band  edge  thus  increasing  the  concentration  of  centers 
available  for  electron  recombination.  The  analysis  of  the  above  data  in  view  of  the  classical 
picture  of  sensitization  suggests  then  that  as  T  is  lowered  there  is  a  sensitization  of  the  DB  states 
by  the  VBT  states,  and  that  around  T^  the  conditions  for  sensitization  [4]  are  fulfilled.  As 
(following  the  LS)  there  are  more  desensitizing  DB  the  conditions  for  sensitization  require  more 
sensitization  centers  i.e.,  in  our  case,  the  decrease  of  T.  In  all  of  this  we  must  assume  that  the 
electron  capture  coefficient  of  the  BT  states  is  considerably  smaller  than  that  of  all  relevant 
(known  as  the  D  and  D*)  DB.  Hence,  all  the  above  results  are  consistent  with  our  suggested 
explanation  of  the  Ttq  and  with  the  observed  ye  >  1  values.  The  question  arises  then  whether  the 
data  shown  in  Fig.  1  can  yield  information  beyond  the  above  general  picture.  To  consider  this 
question  we  will  discuss,  in  the  next  section,  the  T  dependence  of  ye.  This  dependence  is  used 
since,  as  can  be  appreciated  from  Fig.  1,  ye(T)  has  more  conspicuous  features  than  the  (px)e 
dependence.  For  example,  as  far  as  we  know  the  shape  of  yc  has  not  been  considered  before.  We 
will  find  out  then  the  significance  of  the  features  in  the  of  ye(T)  dependence  by  using  numerical 
simulations. 

NUMERICAL  SIMULATIONS 

Let  us  start  by  establishing  the  basic  sensitization  and  thermal  quenching  picture  described 
above.  In  Fig.  2(a)  we  show  the  T  dependence  of  (px)e  for  the  standard  model  of  a-Si:H  with 
parameters  similar  [9,10]  to  those  used  previously  by  Tran  [2],  in  his  “Bl”  model,  but  with  four 
values  of  No.  For  the  highest  No,  we  have,  below  room  temperature  (RT),  almost  a  T 
independent  (px)e,  while  for  the  lowest  No  we  have  a  strongly  dependent  (nx)c.  (The  behavior 
above  RT  has  to  do  with  the  thermal  excitation  of  the  carriers  and  will  be  discussed  elsewhere). 
The  above  results  can  be  simply  interpreted  by  the  fact  that  the  quasi-Fermi  levels  for  electrons 
and  holes  Epn  and  Epp  embrace  all  the  DB  levels  in  the  entire  T  range  so  that  the  concentration  of 
DB  centers  available  for  electrons  recombination.  No’,  does  not  change  in  the  T  range  under 
study.  In  contrast,  the  concentration  of  the  VBT  states  available  for  electron  recombination, 
Pv(T),  is  sensitive  to  the  shift  of  Epp  because  of  its  “movement”  through  the  exponential 
distribution  of  the  density  of  states  (DOS)  in  this  tail.  Hence,  when  we  eliminate  the  DB 
(smallest  ND)  we  get  (considering  the  power-law  T-dependence  of  Py(T)  in  an  exponential  tail 
[4])  that  the  recombination  is  dominated  by  PV(T).  This  interpretation  leads  naturally  to  the 
expected  behavior  of  (px)e  for  the  intermediate  values  of  ND.  For  these  values  there  is  a  T  («Tuj) 
below  which  the  recombination  via  the  available  DB  states  Nd’Cch,  equals  the  recombination  via 
the  VBT  states  Pv’Cvm.  Here  ND’  is  the  concentration  of  “hole”  occupied  DB,  Con  is  their  capture 
coefficient,  PV’(T)  is  the  available  concentration  of  holes  in  the  VBT  states  and  Cvn  is  their 
capture  coefficient.  When  T  is  lowered  so  that  Pv’Cvn  >  No’Con  the  sensitization  is  completed. 
The  latter  is  of  course  true  only  if  Cv„  «  Con  as  discussed  by  Rose  [4].  Correspondingly,  we 
expect  the  transition  from  a  recombination  via  the  DB  to  a  recombination  via  the  VBT  yielding 
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the  effect  of  the  thermal  quenching,  as  is 
clearly  confirmed  by  the  numerical  results 
seen  for  the  intermediate  ND  values.  We 
further  see  that  Tu,  shifts  towards  low  T  with 
increasing  No.  Again,  this  is  in  complete 
agreement  with  the  above  picture  since  the 
above  condition  will  be  fulfilled  then  for 
larger  Pv(T)  values,  i.e.  at  lower  T.  Hence, 
the  above  picture  is  consistent  with  both  the 
experimental  results  and  the  numerical 
simulations. 

All  the  above  becomes  even  much  more 
conspicuous  when  we  consider  the  T 
dependence  of  ye.  This  dependence  is  shown 
in  Fig.  2(b).  When  there  is  no  thermal 
quenching  the  ye  values  are  smaller  than 
unity  while  when  thermal  quenching  is 
present  a  very  distinct  peak  with  a  yc  >  1 
value  is  observed  around  Tu,.  This  is  the 
well-known  manifestation  of  the 
sensitization  effect  [4].  The  clear  shift  of 
the  position  of  the  ye  peak  with  increasing 
Nd  to  a  lower  T,  as  in  Fig.  1  establishes  the 
association  of  this  effect  with  the  thermal 
quenching.  We  further  see  that  with  the 
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Figure  2.  The  computed  T-dependence  of  the 
mobility-lifetime  product  (a)  and  the  light 
intensity  exponent  (b)  for  four  values  of  the 
dangling  bond  concentration  (in  cm*3) 


increase  of  No  the  value  of  yc  decreases  and 

|ts  width,  at  its  peak,  increases.  Again,  this  is  consistent  with  the  experimental  observation  shown 
in  Fig  1,  and  the  above  interpretation,  since  a  “wider”  section  of  the  VBT  (and  thus  a  wider  T- 
range)  is  required  in  order  to  dominate  the  recombination  process.  Hence,  the  basic  picture 
described  in  the  context  of  Fig.  1  is  confirmed  by  the  simulations.  On  the  other  hand,  the  details 
of  the  T  dependences  of  (|ix)e  and  yc  can  convey  many  more  details  regarding  the  recombination 
centers  in  a-Si:H.  In  the  present  case  the  lower  and  wider  the  ye(T)  peak  the  larger  the  ratio 
between  NDand  the  concentration  of  the  VBT  states.  In  the  rest  of  this  section  we  will  derive 
such  semi  quantitative  information  from  the  yc(T)  dependences  around  Tu,  showing  that  one  can 
to  I  low  the  changes  in  the  electronic  structure  in  a  series  of  samples  by  just  measuring  the  above 
easily  measured  phototransport  properties. 

Our  first  check  was  then  to  find  out  the  behavior,  when  the  bandtails  are  eliminated,  in  order 
to  evaluate  the  role  of  the  bantails.  We  see  in  Fig.  3  that  this  elimination  causes  the  absence  of 
sensitization  and  that  above  about  RT  ye(T)  is  determined  by  the  thermal  excitation  of  carriers 
from  the  DB.  The  fact  that  we  have  found  a  similar  result  in  the  presence  of  the  CBT  without  the 
presence  of  the  VBT  establishes  the  conclusion  discussed  above  that  the  VBT  states  are  the 
sensitizing  centers  of  the  DB.  Hence,  for  reasonable  parameters  which  are  consistent  with  all  we 
know  about  the  CBT  states,  these  states  are  not  the  sensitizing  centers  of  the  DB.  We  also 
conclude  from  the  simulations  that  an  increase  in  the  peak  value  (to  ye  >  1)  without  a  significant 
shift  m  Tu,  indicates  a  change  in  P,(T)  without  a  change  in  ND. 
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To  check  the  effect  of  the  first  condition  for 
sensitization,  i.e  that  Cd*/Cv „  will  be  large,  we 
carried  out  a  simulation  with  two  values  of  this 
ratio.  Indeed,  as  seen  in  Fig.  4,  if  this  ratio  is  3  the 
value  of  ye  does  not  exceed  1  while  if  this  ratio  is 
30  it  does  exceed  the  value  of  1 .  We  note  then  that 
since  the  capture  coefficients  are  not  expected  to 
change  significantly  with  preparation  conditions, 
the  above  interpretation  of  the  variation  in  ye 
values  when  a  series  of  samples  is  studied,  is 
valid. 


The  value  of  Pv(T)  is  also  determined  by  the 
energy  separation  of  the  dark  Fermi  level  EF  and 
the  valence  band  edge  Ev.  As  to  be  expected,  the 
smaller  the  Ep-Ev  the  deeper  the  EFp  for  a  given 
illumination  intensity  and  T,  and  thus  the  larger 
the  Pv(T).  Indeed,  the  comparison  made  in  Fig.  5 
of  the  ye  values  for  Ep-Ey  -  0.9  eV  and  EF-Ey  =  1 .3 


o  2  4  6  8  10  12 

iooorr(®K) 

Figure.  3.  The  computed  T-dependence 
of  the  Ye  for  dangling  bonds  in  the  case 
of  negligible  concentration  ofbandtail 
states. 


eV  reveals  this  point.  In  the  former  we  can  get 

sensitization  with  yt  >  1  while  in  the  latter  ye  <  1.  This  is  important  since  it  indicates,  as  was 
found  in  many  experimental  studies  [11-13],  that  in  doped  (n-type)  samples  no  sensitization  or 
t  ermal  quenching  will  be  observed.  Since,  most  of  the  recombination  center  parameters  are 
ixed  in  a-Si:H,  we  can  conclude  that  the  absence  of  sensitization  and  thermal  quenching  is 
usually  a  signature  of  a  large  EF-Ev,  i.e.  of  a  doped  or  contaminated  sample.  Whether  this  is  the 
reason  for  the  absence  of  sensitization  can  be  distinguished  from  the  reasons  mentioned  above, 
by  finding  (from  the  T  dependence  of  the  dark  conductivity  in  a  corresponding  set  of  samples) 
whether  a  change  in  Ep-Ey  took  place  v  ’ 


From  all  the  above  we  can  further  conclude  that  the  appearance  of  a  Yc  >  1  limits  the 
parameter  space  to  an  EF  located  just  above  the  midgap,  to  an  ND  concentration  of  about  10l6cm‘3 
an  to  a  VBT,  of  a  width  of  about  0.05  eV  and  a  DOS  at  the  valence  band  edge  of  about  1021eV 


1000rt*(°K) 

Figure  4.  The  computed  temperature 
dependence  of  the  light  intensity  exponent 
for  two  values  of  the  Con/Cvn. 


1000/T(°K) 

Figure  5.  The  computed  T-dependence  of 
the  light  intensity  exponent  for  two 
positions  of  the  Fermi  level. 


!cin 3.  These  values  are  consistent  with  results  derived  from  very  different  measurements  [2]  and 
thus  shows  that  combining  the  computer  simulation  with  the  experimental  features  associated 
with  the  sensitization  can  indeed  narrow  the  parameter  space  of  the  recombination  centers,  thus 
yielding  quantitative  spectroscopic  information  on  these  parameters. 

CONCLUSIONS 

We  have  shown  that  the  thermal  quenching  and  sensitization  of  the  electrons’  lifetime  in  a- 
Si:H  follows  the  Rose  model  of  two  types  of  recombination  centers.  The  VBT  states  are  the 
sensitizing  centers  of  the  DB.  The  examination  of  the  T  dependence  of  the  light  intensity 
exponent  in  the  vicinity  of  the  thermal  quenching  T  reveals  then  non-trivial  information  on  the 
recombination  centers.  This  is  useful  in  particular  when  a  series  of  materials  is  considered  as 
follows;  A  shift  of  the  ye  peak  towards  low  T  indicates  an  increase  of  the  concentration  of  the 
DB.  A  decrease  and  broadening  of  the  ye  peak  without  a  T  shift  indicates  a  reduction  in  the 
concentration  of  the  VBT  states,  provided  there  is  no  significant  change  in  the  (easily  measured) 
activation  energy  of  the  dark  conductivity.  Hence,  given  a  series  of  samples,  the  simple 
measurements  of  the  T  dependence  of  the  conductivity  and  photoconductivity  can  yield 
semiquantitative  information  on  the  recombination  centers  in  a-Si:H.  It  will  be  also  shown  that 
combining  such  data  with  that  of  the  ambipolar  diffusion  length  can  yield  a  more  quantitative 
information. 
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PLASMONS  ON  LUMINESCENT  POROUS  SILICON  PREPARED  WITH  ETHANOL 


bonds  are  involved  in  the  luminescence  [3,4],  The  study  of  the  plasmon  energies  on  luminescent 
PSi  is  of  interest  for  two  main  reasons.  The  nanostructure  of  PSi  was  shown  [4,5]  to  have  a 
coral-like  structure  consisting  of  a  continuous  hierarchy  of  columns  and  pores.  Typical  lengths 
of  the  columns  are  a  few  tens  of  nanometers,  while  their  average  radius  is  a  few  nanometers.  As 
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ABSTRACT 

The  luminescent  properties  of  porous  silicon  (PSi)  films  in  the  visible  region  were  used  to 

improve  the  photoresponse  of  PSi/Si-wafer  and  PSi/Si  JunCJ^St^0^fsf Vn 
the  reeion  below  500nm.  A  luminescent  PSi  overlayer  was  formed  on  top  of  the  Si  waters  ana  p-n 

junctions*! ty  electrochemical  anod.zat.on.  These  overlayers  have 

t  tL  pci/Qi  wafer  the  PSi  film  was  produced  with  a  high  optical 

prototypes  without  PSi.  Details  about  the  enhancement  process  are  discussed. 

INTRODUCTION 

The  use  of  silicon  for  UV-photodetcction  applications  faces  the  problem  of  the  silicon 
order^of  5  times  smaller  than  at  800nm  around  which  the  detector  has  its  maximum  sensiti  ty. 

IU,  con.ito.bte  e|T« i  Ki'ctoges wiO. 

detectors.  Sonteof.he  studied  photodiode  structures  are.  melnl/PSt/c-St.  ^PS^PS.tc 
71  The  performance  of  those  structures  is  based  on  the  UV  response  for  electnc  g  J? 
along  th^  direction  perpendicular  to  the  layers.  ™amctenstic h»d 
photosensitivity  with  strong  dependence  on  the  reverse  bias  voitag 

interfacefftf  djed  ^  photoresponse  of  lw0  lypes  of  structures  under  UV  illuminatmm  P^Si/AI 
.  pc;/n  c;/n  si/Al  The  first  structure  was  studied  by  measuring  the  change  in  th  8  P’ 
i  on  electric  cnnenl  passing  through  the  p-Si  layer.  I"  «« <•» 
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PSi  overlayer  was  used  as  a  UV-to- visible  conversion  interface.  The  second  heterostructure  was 
studied  by  measuring  the  photovoltage  across  thep-n  junction/PSi  structure 
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,  the  presented  results,  one  can  conclude  that  the  re-emission  efficiency  is  large  enough  to  overpass 

the  small  fraction  of  UV  photons  that  causes  direct  free  carriers  photogeneration  in  c-Si.  As  a 
Figure  4  .  Photoresponse  of  the  first  prototype  (circles)  fma!  remark>  due  t0  the  geometry  of  our  prototypes,  there  is  a  possibility  that  free  carriers 

as  compared  with  the  control  (diamonds)  generated  at  the  PSi  layer  can  diffuse  into  the  active  region.  This  other  mechanism  will  also 

contribute  to  the  enhanced  photoresponse  of  the  detector. 
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Previous  Raman  scattering  studies  of  the  effect  of  hydrogen  on  the  atomic  network  disorder  in 
various  hydrogenated  amorphous  silicon  (a-Si:H)  materials  resulted  in  contradicting  conclusions. 
We  resolve  these  contradictions  by  showing  that  the  surface  and  the  bulk  of  a-Si:H  films  can  behave 
differently  due  to  their  different  hydrogen  contents.  In  particular,  we  establish  that  hydrogen  has  a 
relatively  moderate  effect  in  improving  the  short-range  order  but  a  profound  effect  in  improving  the 
intermediate-range  order  of  the  atomic  network.  ©  1999  American  Institute  of  Physics. 
[S0003-695 1  (99)05244-4] 


The  effect  of  the  disorder  in  the  atomic  network  on  the 
electronic  properties  of  amorphous  semiconductors  has  been 
discussed  extensively  in  the  literature.*’2  In  particular,  for 
hydrogenated  amorphous  silicon  (a-Si:H)  an  obvious  issue 
is  the  role  of  hydrogen  in  the  network  disorder.  While  this 
issue  has  been  considered  already  a  decade  ago,3  it  has  re¬ 
gained  recent  interest4  in  relation  to  the  Staebler— Wronski 
(SW)  effect.5  This  relation  is  a  priori  not  expected  since  the 
SW  effect  has  been  attributed5"7  previously  to  the  local 
breaking  of  “weak”  Si-Si  bonds.  Following  recent  experi¬ 
mental  findings,  suggestions  were  made  that  the  SW  effect 
yields  changes  beyond  the  short-range  (nearest-neighbor)  or¬ 
der  (SRO)  of  the  atomic  network.7'9  While  experiments9 
have  not  shown  that  such  a  change  takes  place  in  a- Si,  they 
have  shown  that  quite  a  significant  change  takes  place  in 
a-Si:H.  Thus,  the  question  of  the  importance  of  hydrogen  in 
the  determination  of  network  order,  in  particular, 
intermediate-range  (beyond  nearest-neighbor)  order  (IRO)  in 
o-Si:H  became  of  wide  interest. 

The  main  tool  for  the  characterization  of  the  degree  of 
order  of  the  atomic  structure  in  a-Si:H  has  been  Raman  scat¬ 
tering  spectroscopy.  It  is  quite  generally  accepted1’3'10  14  that 
in  a-Si  and  a-Si:H  the  features  of  the  Raman  spectrum  that 
are  associated  with  the  crystalline  silicon  transverse-optical 
(TO)  vibration  mode  reflect  the  degree  of  the  SRO.10'"  The 
lower  the  peak  position  of  this  spectral  band,  o>TO ,  and  the 
larger  its  width,  rTO,  the  lower  the  SRO.  On  the  other  hand, 
the  transverse  acoustic  (TA)  mode  that  involves  triads  of 
atoms  is  associated  with  bond  bending.3’14  Correspondingly, 
the  area  under  the  TA  spectral  band  /Ta.  is  considered  to  be 
proportional  to  the  density  of  dihedral  angle 

fluctuations,3’*2'*3  i.e.,  to  the  IRO  in  the  silicon 

network.  "’13’* 4  For  the  experimental  normalization  of  the  lat¬ 
ter  quantity,  one  uses3’12  *3  the  area  of  the  TO  spectral  band 
/T0,  which  represents  the  density  of  the  nearest-neighbor 
bond-angle  fluctuations.'0,11’14 

Turning  to  the  effect  of  hydrogen  on  the  above  three 
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features,  one  would  expect  that  the  higher  the  deposition 
temperature  Ts ,  the  higher  the  degree  of  order  in  the  net¬ 
work.  This  is  because  the  system  can  relax  from  its  meta¬ 
stable  state  to  a  lower-energy  state  that  is  closer  to  the  or¬ 
dered  crystalline-like  state.  Indeed,  for  a -Si  films,  it  has  been 
found'2-'5  that  both  Tjo  and  IjA/Ijo  decrease  with  increas¬ 
ing  Ts .  In  a-Si:H  we  have,  however,  in  addition,  the  effect 
of  the  hydrogen  content  CH ,  which  decreases  with  the  in¬ 
crease  of  Ts.  Maley  and  Lannin3  found  that  both  the  rT0 
and  /TA//T0  features  decrease  with  CH .  They  suggested  that 
the  decrease  in  atomic  network  order  with  decreasing  CH 
(which  follows  the  increase  of  Ts,  hereafter,  the  “H  role”) 
offsets  the  increase  of  this  order  due  to  the  above  thermal 
relaxation  effect  (hereafter  the  “7"  role”).  While  the  same 
behavior  was  found  in  some  studies,9’*2  the  reverse  behavior 
has  been  found  in  others.'516  Our  own  comprehensive 
study17  on  samples  deposited  by  glow-discharge  and  sputter¬ 
ing  techniques  has  shown  a  very  small  effect  of  Ts  on  both 
the  SRO  and  the  IRO.  These  results  could  be  interpreted  as 
caused  by  a  delicate  balance  between  the  above  two  roles. 

The  above  interpretation  appears,  however,  very  unsatis¬ 
factory  when  one  considers  the  results  reported  for  materi¬ 
als  prepared  by  the  “hot-wire”  technique.'  ’  In  these  ma¬ 
terials,  CH  varies  by  about  two  orders  of  magnitude  (20-0.2 
at.  %  of  hydrogen)  when  Ts  is  raised  from  50  to  550  °C.  This 
makes  these  materials  excellent  candidates  for  the  study  of 
the  H  role  in  the  atomic  network  order.  In  sharp  contrast  with 
indications  from  transport  measurements,21  the  earlier 
conclusion3  and  the  theoretical  expectations,22  the  Raman 
spectra18  reported  for  these  materials  have  hardly  shown  any 
change  in  rjo>  as  a  function  of  Ts,  in  spite  of  the  large 
change  in  .  (The  a»xo  an£I  ^ta^to  features  have  not 
been  reported.)  Moreover,  results  of  x-ray  diffraction  mea¬ 
surements  have  been  interpreted'8  to  indicate  the  increase  of 
IRO  with  Ts  in  the  “hot-wire”  materials  (see  below). 

Following  this  rather  puzzeling  behavior,  we  have  car¬ 
ried  out  a  comprehensive  Raman  scattering  study  of  the 
“hot-wire”  materials.  These  materials  that  are  very  interest¬ 
ing  in  their  own  right,'9’2'’23  exhibit  a  “weaker”  SW  degra¬ 
dation  than  the  conventional  “glow-discharge  materials.” 
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FIG.  1 .  Raman  spectra  of  films  prepared  by  the  hot-wire  technique  as  de¬ 
termined  by  using  an  excitation  source  of  light  wavelength  of  514.5  nm.  The 
samples  differ  by  the  substrate  temperature  Ts  at  which  they  were  deposited. 
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FIG.  3.  Raman  spectra  of  films  prepared  by  the  hot-wire  technique  as  de¬ 
termined  by  using  an  excitation  source  of  light  wavelength  of  583.2  nm.  The 
samples  differ  by  the  substrate  temperature  Ts  at  which  they  were  deposited. 


The  series  of  “hot-wire”  films  used  here  was  deposited  and 
characterized  (especially  for  their  hydrogen  content19,20)  in 
the  same  laboratory  at  which  the  previous  Raman  scattering 
studies  were  carried  out.18  Applying  the  micro-Raman  tech¬ 
nique  and  data  analysis  of  our  previous  studies17,24  of 
a-Si:H,  we  repeated  first  the  previous  measurements18  using 
the  standard  5 14.5-nm- wavelength  illumination  of  an 
Ar+-ion  laser  as  the  excitation  source. 

The  spectra  obtained  are  shown  in  Fig.  1.  In  these  spec¬ 
tra  the  “TO  band”  (around  480  cm-1)  and  the  “TA  band” 
(around  160  cm-1)  are  the  conspicuous  features  of  the  data. 
The  results  of  the  analysis  of  the  spectra  shown  in  Fig.  1  are 
presented  in  Fig.  2.  It  is  seen  that  there  is  a  slight  decrease  of 
Fto  [Pig*  2(a)],  an  increase  in  o>TO  [Fig.  2(b)],  and  a  de¬ 
crease  of  /TA//To  [Pig-  2(c)]  with  Ts ,  In  order  to  emphasize 
the  relation  between  Tto  and  the  SRO,  we  have  added  in 
Fig.  2(a)  the  rms  bond-angle  deviations  as  suggested  in  Refs. 
10  and  11.  All  these  data  indicate  the  dominance  of  the  “7 
role”  over  the  “H  role,”  confirming  the  previous  data18  and 


FIG.  2.  Dependence  of  the  main  features  of  the  Raman  spectra  on  the 
substrate  temperature  at  deposition  Ts  as  derived  from  the  data  shown  in 
Fig.  1. 


leaving  us  with  the  above  “puzzle”  of  the  minute  “H  role” 
in  spite  of  the  large  change  in  CH . 

We  found  a  plausible  clue  to  this  puzzle  noting  that  CH 
near  the  surface  of  “hot-wire  films”  is  higher  than  in  the 
bulk,  and  that  the  surface  is  weakly  effected  by  the  out- 
diffusion  of  the  hydrogen  during  the  deposition  process.20 
Since  the  514.5  nm  radiation  is  absorbed  within  50  nm,  it 
seemed  to  us  reasonable  to  explain  all  the  above  results  as 
reflecting  the  network  disorder  at  the  films’  surface,  where 
the  Ch  is  varied  only  by  a  factor  of  2  (Ref.  20)  throughout 
the  wide  Ts  range  employed  in  this  study.  Hence,  we  con¬ 
jectured  that  the  large  CH  variation  with  Ts  is  a  bulk  effect 
that  could  not  have  been  probed  by  the  present  and 
previous18  Raman  studies,  which  utilized  the  514.5  nm  laser 
light. 

To  test  this  conjecture  we  have  repeated  our  measure¬ 
ments  but  with  a  laser  light  of  a  wavelength  of  583.2  nm.  At 
this  wavelength  the  absorption  depth  is  300  nm.  Noting  that 
the  samples’  thickness  is  of  the  order  of  1  /xm,  the  observed 
behavior  can  be  considered  to  be  associated  with  the  bulk  of 
the  film.  Furthermore,  we  note  that  this  depth  is  much  be¬ 
yond  the  range  where  the  CH  gradient  takes  place.20 

For  obtaining  the  583.2  nm  radiation,  we  have  employed 
a  dye  laser  (CR-599,  rhodamine  6G)  pumped  by  an  Ar+-ion 
laser.17  The  spectra  obtained  with  this  “bulk  probe”  are 
shown  in  Fig.  3,  and  the  analyzed  7 s  dependence  of  the  main 
spectral  features  is  presented  in  Fig.  4.  Considering  the 
above  discussion,  both  Figs.  4(a)  and  4(b)  indicate  an  overall 
decrease  in  the  SRO  with  increasing  Ts .  Similarly,  Fig.  4(c) 
indicates  a  decrease  of  the  IRO  with  increasing  Ts  (decreas¬ 
ing  C„).  The  latter  effect  is  significantly  stronger  in  compari¬ 
son  with  the  one  found  in  a-Si:H  materials,  which  were  de¬ 
posited  by  other  techniques.17  Furthermore,  the  /TA//T0 
ratios  observed  for  the  higher  Ts  (smaller  CH)  values  are 
larger  than  any  ratio  ever  encountered  before  in  a -Si  and 
a-Si:H.  All  the  above  results  show  then  a  clear  domination 
of  the  “H  role”  over  the  “7  role”  when  the  variation  in 
is  large.  We  can  reconcile  our  conclusion  with  the  x-ray 
diffraction  results18  (see  above),  that  also  sample  the  bulk  of 
the  material,  by  noting  that  such  results  have  been  attributed 
previously25  to  the  SRO  parameter  that  is  manifested  in  the 
Raman  spectrum  by  Tto.  If  we  adopt  the  latter  interpreta- 
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FIG.  4.  Dependence  of  the  main  features  of  the  Raman  spectra  on  the 
substrate  temperature  at  deposition  Ts  as  derived  from  the  data  shown  in 
Fig.  3. 

tion,  we  get  a  consistent  picture  of  the  x-ray  data  and  our 
data  as  shown  in  Fig.  4(a),  i.e.,  that  for  the  same  Ts  range 
(200-400  C)  the  SRO  is  slightly  increasing  with  T s . 

The  above  solution  of  the  “hot-wire”  puzzle  justifies 
then  the  physical  interpretation  of  Maley  and  Lannin3  show¬ 
ing  that  the  effect  of  hydrogen  is  much  stronger  on  the  num¬ 
ber  of  bond-bending  mode  fluctuations  than  on  the  number 
of  bond-angle  mode  fluctuations.  In  addition,  the  present 
study  solves  the  problem  of  the  contradicting  reports  in  the 
literature  since  it  shows  that  all  the  previously  reported  re¬ 
sults  obtained  by  using  the  common  514.5  nm  laser  line  are 
associated  with  the  “surface”  of  a-Si:H,  and  that  the  at¬ 
tempts  to  correlate  them  with  defect  concentration  and  trans¬ 
port  in  the  bulk  should  be  examined  with  care. 

Our  results  are  also  a  further  confirmation  of  our  general 
conclusion26  that  in  a-Si:H  there  is  always  a  “surface”  (a 
few  tens  of  nm  thick)  layer,  which  has  a  lower  SRO  and  IRO 
than  the  “bulk.”  However,  we  can  add  now  that  once  CH  is 
significantly  different  in  them,  as  in  the  “hot-wire”  materi¬ 
als,  the  reverse  may  occur . 

In  summary,  we  have  established  experimentally  that  hy¬ 
drogen  causes  a  relaxation  of  the  atomic  network  in  a-Si:H 
and  that  this  effect  is  much  more  apparent  for  the 
intermediate-range  order  than  for  the  short-range  order.  This 
reconfirmation  resolves  the  seemingly  contradicting  avail¬ 
able  data  as  follows.  It  appears  that,  quite  generally,  the  sur¬ 
face  has  a  near-constant  CH ,  and  thus,  the  degree  of  its  order 
is  determined  by  the  temperature-induced  network  relax¬ 
ation.  Since  in  the  bulk  the  value  of  CH  depends  strongly  on 
the  deposition  conditions,  the  role  of  hydrogen  on  the  net¬ 


work  disorder  is  much  more  pronounced.  It  is  expected  then 
that  intermediate-range-order  effects  in  the  network  (as  has 
been  suggested  for  the  SW  effect)  are  associated  with  the 
presence  of  hydrogen. 
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Abstract 


HT  fng  rep°rtf  c0ncermng  the  effect  oflight  soaking  (LS)  on  the  structure  of  hydrogenated  amorphous 
s.l.con  (a-S.:H)  we  have  earned  out  a  comprehensive  Raman  scattering  study  of  this  effect  on  a-Si:H  films  prepared  by 
Atoent  (tepositum  techniques.  We  have  found  that  first,  the  LS  causes  a  major  rearrangement  of  atoms  as  suggested  by 

hnnH  iu  T '“I!  m0de  S'  SeCOnd'  the  short'range  order  decreases  upon  LS  as  we  expect  from  the  breaking  of  Si-Si 
a  d  Thlrd  and  counterintuitive,  the  intermediate  range  order  may  increase  upon  LS.  ©  2000  Elsevier  Science  B  V 
All  rights  reserved. 


1.  Introduction 

The  effect  of  illumination  induced  degradation 
on  the  electronic  properties  of  amorphous  semi¬ 
conductors,  known  as  the  Staebler-Wronski  (SW) 
effect,  has  been  a  subject  of  study  in  the  last  20 
years  [1,2].  In  spite  of  its  importance  for  under¬ 
standing  the  covalent-amorphous  system  and  in 
spite  of  its  relevance  to  opto-electronic  applica¬ 
tions,  it  is  still  poorly  understood  and  hardly 
controlled.  While  the  recovery  of  the  SW  effect  by 
annealing  indicates  that  the  effect  is  associated 
with  a  modification  of  the  covalent  amorphous 
network,  the  consideration  of  the  structural 
mechanism  underlying  the  SW  effect  has  become  a 
subject  of  interest  only  in  the  last  few  years  [3],  In 
contrast  with  the  traditional  models  of  the  SW 
effect  that  were  concerned  with  the  single  dangling 
bond  and  its  nearest  neighbors,  there  is  now  ex- 
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perimental  evidence  [4-12],  supported  by  theoret¬ 
ical  calculations  [13,14],  that  indicates  that  the  SW 
effect  is  a  global  effect.  This  effect  is  a  result  of  a 
reorganization  of  the  network,  which  follows  the 
carrier  excitation  and/or  recombination. 

While  the  very  fact  that  this  reorganization  of 
the  atomic  network  appears  now  to  be  a  funda¬ 
mental  property  of  the  SW  effect,  its  description  is 
still  at  a  rudimentary  level  in  our  opinion.  In  fact 
while  the  theoretical  predications  for  both  unhy¬ 
drogenated  amorphous  silicon  (a-Si)  and  hydro¬ 
genated  amorphous  silicon  (a-Si:H)  have  dealt 
with  this  reorganization,  they  did  not  address  ex¬ 
plicitly  the  question  of  the  order.  We  note  however 
that  recent  theoretical  calculations  [4,5,14]  of  the 
SW  effect  have  implied  the  (intuitively  expected) 
increase  of  the  root  mean  square  (rms)  of  the  bond 
angle  distribution,  i.e.,  a  decrease  of  the  short- 
range  order  (SRO)  upon  light  soaking  (LS)  in 
a-Si.H.  Experimental  data  and  variable  coherence 
transmission  electron  microscopy  (VC-TEM)  re¬ 
sults  [15]  (which  do  not  monitor  the  SRO)  sug¬ 
gested  that  the  LS  causes  a  decrease  in  the 
intermediate  range  order  (IRO),  i.e.,  the  order 
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beyond  the  nearest  neighbor,  upon  the  LS  of 
a-Si:H.  On  the  other  hand,  the  results  of  the  most 
common  method  to  follow  the  degree  of  both  SRO 
and  IRO,  i.e.,  the  Raman  scattering  technique,  are 
not  conclusive.  Considering  that  the  effect  of  LS  in 
one  work  [16]  did  not  show  a  change  in  the  fea¬ 
tures  (see  below)  that  account  for  the  SRO  and 
IRO,  in  another  work  [17],  a  decrease  in  the  SRO 
(in  agreement  with  the  theoretical  expectations 
[4,5,14])  and  an  increase  in  the  IRO  (in  contrast 
with  the  VC-TEM  data  [15]  and  these  expectations 
[4,5])  was  found.  In  particular  we  note  that  dif¬ 
ferent  samples  and  different  techniques  have  been 
used  for  the  evaluation  of  this  effect  so  that  all 
these  findings  do  not  sum  to  consistent  experi¬ 
mental  results  as  to  the  effect  of  light  soaking  on 
the  network  disorder  in  a-Si:H.  In  the  present 
work  we  intended  to  resolve  this  problem  and  to 
provide  a  consistent  experimental  analysis  of  the 
evolution  of  the  disorder  upon  LS,  by  using  the 
same  technique  on  different  a-Si:H  materials. 


2.  Experimental 

We  have  carried  out  our  study  on  samples 
prepared  by  different  deposition  techniques  and 
different  laboratories.  The  first  set  of  samples  used 
in  our  study  were  1  pm  thick,  device  quality,  a- 
Si:H  films,  whose  deposition  has  been  described 
previously  [18,19].  The  substrate  used  for  the 
samples  was  a  7059  Corning  glass  and  the  depo¬ 
sition  technique  was  the  dc  glow  discharge  (GD) 
decomposition  of  silane  with  a  substrate  temper¬ 
ature,  rs,  of  260°C.  We  have  also  reported  previ¬ 
ously  [18,19]  the  effects  of  light  soaking  on  the 
electronic  properties  of  these  films.  These  effects 
are  similar  to  those  well  documented  in  the  liter¬ 
ature  [1,2].  The  second  set  of  samples  were  pre¬ 
pared  by  the  hot  wire  (HW)  technique  with 
thickness  and  substrate  similar  to  the  GD  samples. 
For  this  set  we  have  examined  samples  prepared 
with  varying  substrate  temperatures,  Ts.  Again,  we 
have  reported  the  deposition,  the  electronic  prop¬ 
erties  and  the  effect  of  LS  on  these  properties 
previously  [20,21]. 

The  micro-Raman  set-up  and  the  measurement 
procedure  and  data  analysis  used  in  the  present 


study  have  been  described  [22,23].  The  514.5  nm 
laser  illumination  was  obtained  from  an  Ar-ion 
laser.  For  the  present  measurements,  however,  the 
power  of  the  incident  radiation  was  7.5  mW,  the 
beam  diameter  was  5  pm  and  the  measurement 
was  carried  out  for  300  s  to  avoid,  as  much  as 
possible,  the  light  soaking  by  the  ‘Raman-probe’ 
beam.  The  controlled  light  soaking  was  carried  out 
using  a  (100  W  Oriel  Q)  lamp  with  a  broad  band 
quartz  fiber  optics  under  a  power  of  400  mW/cm2. 
This  illumination  was  applied  for  2  h.  The  mea¬ 
surements  and  light  soaking  were  carried  out  at 
room  temperature  (300  K)  and  at  liquid  nitrogen 
temperature  (70  K).  The  properties  found  were 
much  the  same  at  both  temperatures. 


3.  Experimental  results 

Turning  to  the  results  of  our  measurements  we 
show  in  Fig.  1  the  Raman  spectra  obtained  before 
and  after  light  soaking  at  70  K  for  a  GD  sample. 
Similar  data  were  obtained  at  300  K  and  for  the 
HW  samples.  We  note,  as  is  common  in  the 


Fig.  1 .  The  Raman  spectra  of  a  sample  of  the  GD  set  before  (a) 
and  after  (b)  light  soaking. 
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upon  LS  in  accordance  with  the  above  and  the 
previous  results  [17].  This  conclusion  applies  to  all 
the  samples  (50  ^  Ts  <  520°C)  used.  On  the  other 
hand,  the  results  associated  with  the  IRO  were  the 
opposite  of  those  given  above  for  the  GD  material. 
This  difference  is  shown  by  the  summary  of  the 
corresponding  results  which  are  presented  in  Fig. 
2.  In  this  figure  we  see  a  decrease  of  the 
feature  upon  LS.  This  change  as  explained  above, 
indicates  the  increase  of  the  IRO  with  LS.  We  see 
then  that  in  agreement  with  Ref.  [17]  the  LS  causes 
an  increase  in  the  IRO  in  spite  of  the  fact  that  there 
is  a  decrease  in  the  SRO. 


4.  Discussion  and  conclusions 

From  our  results  we  can  conclude  that  the  re¬ 
organization  of  the  atomic  network  in  a-Si:H,  due 
to  the  SW  effect,  is  accompanied  by  the  decrease  in 
the  SRO  and  a  change  in  the  IRO.  The  observed 
opposite  directions  of  this  change  have  to  do  with 
the  more  subtle  properties  of  the  material  under 
study.  This  conclusion  explains  then  that  the  dif¬ 
ferences  between  previous  results  are  due  to  the 
differences  in  the  materials  and  not  to  the  differ¬ 
ence  in  the  experimental  techniques  and  their  in¬ 
terpretation. 

The  above  conclusion  on  the  two  possible 
properties  of  the  IRO  in  a-Si:H  has  to  be  consid¬ 
ered  when  models  of  the  SW  effect  are  proposed. 
In  particular  the  simultaneous  decrease  of  the 
SRO,  with  the  possible  increase  or  decrease  of  the 
IRO,  has  to  be  accounted  for  in  any  future  model 
of  the  SW  effect.  The  changes  in  SRO  are  expected 
due  to  the  fact  that  the  weak  Si-Si  bond,  or  the  Si¬ 
ll  bond,  retain  the  similarity  to  the  c-Si  tetrahedral 
configuration  better  than  a  Si  atom  with  a  dan¬ 
gling  band.  On  the  other  hand,  the  fact  that  the 
IRO  change  does  not  follow  the  SRO  change  in¬ 
dicates  that  the  atomic  network  rearrangement 
that  follows  is  not  trivial.  Tentatively,  we  explain 
these  changes  by  a  competition  of  two  effects.  The 
first  is  that  the  breaking  of  bonds  is  also  the 
breaking  of  dihedral  bonds.  This  effect  will  de¬ 
crease  the  correlation  between  adjacent  tetrahedra 
and  thus  produce  a  decrease  of  the  IRO.  On  the 
other  hand,  the  breaking  of  the  weaker  bonds  can 


be  followed  by  a  relaxation  of  the  network  in  the 
environment  of  the  broken  bond.  If  the  network  is 
not  relaxed  the  breaking  of  bonds  by  the  LS, 
which  is  manifested  by  the  decrease  of  the  SRO, 
may  lead  to  a  more  relaxed  network.  If  this  effect 
overcomes  the  decrease  in  the  correlation,  the  end 
result  will  be  an  increase  in  the  IRO.  On  the  other 
hand,  if  the  network  is  a  priori  more  relaxed,  the 
decrease  in  the  correlation  will  dominate  and  an 
overall  decrease  of  the  IRO  will  follow.  An  inde¬ 
pendent  check  of  this  suggestion  can  be  made  by 
various  studies.  For  example  a  study  of  the  va¬ 
lence-band  tail  width  in  materials  that  show  dif¬ 
ferent  effects  of  LS  on  the  IRO  is  called  for  since  it 
is  expected  that  a  wider  tail  is  associated  with  a  less 
relaxed  network. 
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Abstract 

We  have  synthesized  green  and  red  luminescent  silicon  nanocrystals  in  a  Si02  matrix  by  RF  co-sputtering  on  a  quartz 
substrate.  The  transmission  coefficient  measurements  were  used  to  estimate  the  nanocrystal  size  distribution.  The  size 
distribution  reveals  peaks  in  the  range  1.1 -2.6  nm  with  a  long  tail  towards  the  larger  size.  As  the  nanocrystal  size  reduces 
photolum.nescence  spectrum  shifts  from  red  to  green  wavelengths.  The  measured  PL  emission  energy  is  in  agreement 
with  the  corrected  LDA  calculations.  With  decreasing  nanocrystal  size,  the  phonon  Raman  spectra  exhibit  softening 
accompanied  with  ncreasing  asymmetrical  broadening.  The  observed  line  shape  is  explained  by  considering  phonon 
confinement  in  a  spherical  nanocrystal.  The  major  contribution  to  the  phonon  line  shape  comes  from  those  nanocrystals 
that  favor  resonance  interaction  with  either  incoming  or  outgoing  photon.  ©  1999  Elsevier  Science  B.V.  All  rights 

reservpH  & 


Keywords:  Silicon  nanocrystals;  Si-Si02  film;  Photoluminescence;  Raman  scattering;  Quantum  confinement 


1 1.  Introduction 

Nanometer  size  silicon  crystals  show  marked  im¬ 
provement  in  photoluminescence  (PL)  quantum  ef¬ 
ficiency  compared  to  bulk  silicon  and  are,  therefore, 
attractive  for  silicon  based  optoelectronic  applica¬ 
tions  [1],  The  PL  observed  from  Si  nanocrystals  is 
attributed  to  quantum  size  effects.  The  effect  of 
carrier  confinement  and  band  gap  upshift  in  silicon 
nanocrystals  has  been  discussed  intensively  in  the 
i  literature  [1],  The  photoluminescence  energy  for 
silicon  nanocrystals  have  been  calculated  by  using 
the  empirical  pseudopotentials  (EPS)  [2],  the 
I  third-nearest-neighbor  tight  binding  [3]  and  time- 
dependent  tight  binding  techniques  [4],  Despite  the 
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apparent  disagreement  in  predicting  the  PL  emis¬ 
sion  energies,  these  calculations  demonstrate  that 
the  PL  energy  increases  as  the  nanocrystal  size 
decreases  and  the  dominant  contribution  to  the 
visible  light  emission  comes  from  nanocrystals 
smaller  than  2  nm.  The  resonantly  excite  PL  spec¬ 
trum  has  shown  long  radiative  lifetime  in  the  Si 
nanocrystals  and  provides  an  important  evidence 
that  the  phonons  are  directly  involved  in  the 
radiative  recombination  process  [5].  Clearly,  this 
suggests  that  the  Si  nanocrystal  electronic  band 
structure  remains  indirect  type.  Theoretical  calcu¬ 
lations,  however,  indicate  that  quasi-direct  gap  in 
nanocrystal  silicon  is  possible  for  sizes  1.0-1. 5  nm 
in  diameter  [3], 

Raman  scattering  technique  is  commonly  used 
for  semiquantitative  determination  of  size  effects  on 
vibrational  modes,  such  as  confined  optical 
phonon,  surface  phonon  and  confined  acoustic 
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phonon,  in  nanocrystals  [6-9].  An  adequate  know¬ 
ledge  of  the  bulk  phonon  dispersion  is  a  prerequi¬ 
site  for  understanding  lattice  dynamical  properties 
and  electron-phonon  interactions  in  nanocrystals. 
The  effect  of  carrier  confinement  on  optical 
phonons  in  silicon  nanoparticles  has  been  widely 
studied  by  Raman  scattering  [6-8].  The  size  effect 
is  reflected  in  a  finite  wave  vector  of  nanocrystal 
and  corresponding  lowering  of  phonon  frequency 
from  those  of  the  bulk.  A  quantum  confinement 
model  that  takes  into  account  the  finite  phonon 
wave  vector  and  Gaussian  size  distribution  for 
spherical  nanocrystals  is  generally  invoked  to  ex¬ 
plain  the  observed  Raman  line  shape.  Here,  we 
investigate  the  size-dependent  photoluminescence 
and  Raman  scattering  from  silicon  nanocrystals 
dispersed  in  a  Si02  matrix  grown  on  a  quartz 
substrate.  We  emphasize  the  role  of  nanocrystals 
that  favor  the  resonant  interaction  with  incoming 
or  outgoing  photon. 


2.  Experiment 

The  Si-Si02  film  was  synthesized  on  the  quartz 
substrate  by  RF  co-sputtering.  The  substrate  was 
a  6-in  long  quartz  strip.  The  chamber  was  first 
evacuated  down  to  1.0  x  10' 7  Torr  and  the  depos¬ 
ition  was  performed  under  an  argon  pressure  of 
2.0x10  2  Torr  at  deposition  temperature  of 
115  C.  On  a  Si02  target  (6  in  in  diameter)  small 
pieces  of  Si  tips  (area  ~  0.1  in2)  were  placed  and 
they  were  co-sputtered  for  3  h.  The  number  of  Si 
tips  and  their  position  on  the  Si02  were  used  to 
control  the  size  of  silicon  nanoparticles.  The  film 
was  subsequently  annealed  in  nitrogen  ambient  at 
1100°C  for  30  min.  The  deposited  film  has  an  area 
of  5  x  0.2  in  and  thickness  of  the  order  of  1.7  pm. 
Both  the  size  and  density  of  silicon  nanocrystals 
vary  continuously  along  the  length  of  the  film.  The 
transmission  measurements  were  made  in  the  range 
190-820  nm  using  a  HP-8452A  spectrophotometer. 
The  photoluminescence  and  Raman  measurements 
were  performed  in  a  quasi  back-scattering  geo¬ 
metry  at  room  temperature  using  the  514.5  nm  line 
of  an  Ar  laser.  The  spectra  were  dispersed  by 
a  TRIAX  320  spectrometer  and  detected  by 
a  cooled  photomultiplier  with  GaAs  photocathode. 


3.  Results 

Fig.  1  shows  typical  transmission  spectra  taken 
from  different  positions  on  the  Si-Si02  film.  The 
deposited  film  is  divided  in  50  equal  segments  and 
designated  a  number,  number  1  is  assigned  to  Si02 
rich  end  and  50  to  the  Si  rich  end  of  the  film.  Each 
of  these  segments  has  different  nanocrystal  density 
and  size.  In  order  to  study  the  size  dependence  of 
optical  properties,  we  measured  the  transmission 
spectrum  from  each  one  of  these  segments.  The 
optical  absorption  edge  shifts  to  higher  energy  as 
one  moves  from  silicon  rich  side  (position  50)  to 
Si02  rich  side  (position  1)  of  the  sample  due  to 
nanocrystal  size  reduction.  The  transmission  coef¬ 
ficient  of  the  film  is  related  to  the  energy  absorbed 
by  the  Si  nanocrystal  of  size  R  imbedded  in  the 
Si02  matrix.  Assuming  that  the  energy  upshift  is 
given  by  A/R13  ( A  is  a  constant),  and  the  absorp¬ 
tion  is  primarily  due  to  interband  transitions 
between  the  quantum  confined  states  of  spherical 
Si  nanocrystals,  the  size  distribution  can  be 


Fig.  1.  Transmission  spectra  from  different  positions  of  a 
Si-SiO,  film  annealed  at  1 100°C.  The  film  is  divided  in  50  equal 
segments;  the  number  on  each  spectrum  designates  the  segment 
number.  The  number  1  and  50  are  assigned  to  the  SiO,  rich  and 
Si  rich  end  of  the  film,  respectively.  The  absorption  edge  shifts  to 
lower  wavelengths  with  decreasing  silicon  nanocrystai  size  and 
density. 
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Fig.  2.  Calculated  nanocrystal  size  distribution  in  Si-Si02  films 
annealed  at  1 100°C  at  the  positions  where  transmission  spectra 
were  taken.  The  size  distribution  is  narrower  for  smaller  sizes, 

approximated  with  the  first  derivative  of  the  trans¬ 
mission  spectrum  [10].  We  have  used  the  transmis¬ 
sion  measurements  to  estimate  the  nanocrystal  size 
distribution  in  the  film,  the  estimated  size  distribu¬ 
tion  is  asymmetric  with  a  long  tail  towards  the 
larger  size  and  centered  around  1. 1-2.6  nm  as 
shown  in  Fig.  2. 

Fig.  3  shows  the  PL  emission  spectra  from 
Si-Si02  film  excited  with  the  green  514.5  nm 
(2.412  eV)  line  of  an  argon  ion  laser.  The  PL  spec¬ 
trum  is  broad  and  shifts  to  higher  energy  with 
decreasing  size.  In  addition,  there  is  a  large  reduc¬ 
tion  in  integrated  PL  intensity  for  smaller  sizes  due 
to  lower  density.  The  peaks  are  observed  at  782  nm 
(1.59  eV)  and  601  nm  (2.06  eV)  for  dot  sizes  of  2.6 
and  1.4  nm,  respectively.  According  to  Delerue 
et  al.  [11],  EPS  or  tight  binding  as  well  as  corrected 
ab  initio  local  density  approximation  (LDA)  tech¬ 
niques  give  reliable  predictions  for  Si  nanocrystal 
band  gap.  It  has  also  been  pointed  out  that  the 
quasi  particle  gaps  and  exciton  coulomb  energies 
of  Si  nanocrystals  require  a  large  correction  to 
the  band  gap  which  are  not  included  in  the  one 
electron  theory.  The  exciton  gap  is  defined  as 


500  600  700  800  900 
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Fig.  3.  Size-dependent  PL  spectra  from  silicon  nanocrystals 
embedded  in  SiO,.  The  average  size  of  the  nanocrystal  is  cal- 
culated  from  the  transmission  data. 

£“c  =  £^p  —  EC0U|,  where  the  quasi-particle  gap 
p  is  the  difference  in  energies  resulting  from  the 
separate  addition  of  an  electron  and  hole  to  the 
system  while  the  correction  £coul  results  from  the 
corresponding  attraction  between  these  two  quasi¬ 
particles.  Our  PL  results  are  in  agreement  with 
corrected  LDA  calculation  [11],  We  attribute  the 
featureless  broad  PL  spectrum  to  band-to-band 
recombination  in  Si  nanocrystal  and  fluctuation  in 
their  sizes.  For  the  1.4  nm  nanocrystal,  we  also  see 
structure  towards  the  high  energy  (2.2  eV)  due  to 
resonant  excitation  with  laser  line.  It  should  be 
noted  that  PL  energy  is  lower  than  the  absorption 
edge  at  each  segment  of  the  film,  and  this  difference 
is  appreciably  large  towards  the  Si02  side  of  the 
film  where  the  nanoparticle  sizes  are  smaller.  Self- 
trapped  exciton  at  the  interface  between  Si  and 
Si02  or  a  defect  state  is  considered  to  be  respon¬ 
sible  for  the  observed  Stokes  shift  of  PL  emission 
energy  [11], 

The  Raman  spectra  from  Si-Si02  film  were  ex¬ 
cited  with  514.5  nm  (2.41  eV)  laser  line  of  an  argon 
ion  laser  from  the  region  that  emits  green-orange 
luminescence.  Fig.  4  shows  representative  spectra 
from  three  different  positions  of  the  sample  charac¬ 
terized  by  the  average  nanocrystal  diameter.  The 
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observed  spectra  show  phonon  softening  and  line 
shape  broadening  with  decreasing  size.  As  the  size 
decreases  from  2.2  to  1.4  nm,  the  Raman  peak  shifts 
from  511  to  502  cm-1.  The  Raman  lines  are  asym¬ 
metrically  broadened  with  a  tail  towards  the  low- 
energy  side.  In  the  uncorrected  data,  the  integrated 
intensity  remains  fairly  constant.  The  broad  line 
shape  arises  from  size  distribution  in  an  ensemble 
of  nanocrystals,  which  leads  to  a  large  fluctuation 
in  wave  vectors,  particularly  for  size  smaller  than 
2  nm. 

In  order  to  explain  the  observed  resonance  line 
shape,  we  assume  that  the  optical  vibrations  are 
confined  to  the  spherical  nanocrystal  with  effective 
wave  vector, 


where  n„  is  nth  node  of  the  spherical  Bessel  function 
and  R  is  the  nanocrystal  radius.  A  quadratic 
bulk  Si  dispersion  gives  the  eigen  frequencies, 

col(R)  =  col-pql{RX  (2) 

where  co0  is  the  bulk  silicon  optical  phonon 
frequency  (52° 5 cm"1)  and  p  is  a  parameter 
(6.977  x  10  s  *)  describing  the  dispersion  of  the 
optical  phonon  in  the  bulk  silicon.  The  observed 
Raman  frequencies  are  in  agreement  with  the  cal¬ 
culated  eigenfrequencies  of  optical  vibrations  of 
zero  angular  momentum  (/  =  0)  mode  as  a  function 
of  R  in  a  spherical  Si  nanocrystal  embedded  in 
the  Si02.  It  should  be  pointed  out  that  a  precise 
determination  of  nanocrystal  size,  particularly  for 
R  <  1  nm,  is  necessary  for  the  estimation  of  equiva¬ 
lent  wave  vectors  and  corresponding  eigen- 
frequencies.  The  Raman  line  shape  can  be  written 
as  /(«,  R)  [9], 

7(cos,  R)  =  /0 


Raman  Shift  (cm*1) 


Fig.  4.  Size-dependent  Raman  spectra  from  silicon  nanocrys¬ 
tals.  The  average  size  of  the  nanocrystal  is  calculated  from  the 
transmission  data. 

elements  over  all  the  intermediate  states,  R  is  the 
mean  nanocrystal  radius,  E„,(R)  and  E„2(R)  are 
the  eigenenergy  for  R  in  the  /q  and  fi2  intermediate 
states,  respectively,  r ,,(/’)  is  the  intermediate  state 
(phonon)  lifetime  broadening,  co,  and  cos  are  the 
excited  laser  energy  and  scattered  photon  energy, 
respectively.  In  the  incoming  resonance  process, 
co/  coincides  with  one  of  the  electronic  states  £/q(R) 
while  outgoing  resonance  occurs  when  co,  ap¬ 
proaches  E/i2.  The  resonance  conditions  is  de¬ 
scribed  as  hco,  =  Efii(R)  for  incoming  resonance, 


[.fico/  hco i,  hco„(R)]2  +  r2  ’  ^  and  ~  hco„(R)  =  Efi,(R)  for  outgoing  reso- 

where  /„  is  the  product  of  scattering  a  nance;  We  assume  that  the  intermediate  electronic 

eigenvalues  of  f*  “  an  e,ectro"-hole  Pa*  «ate  and  the  elec 

appropriate  scattering  matrix  tron-phonon  couple  via  deformation  potential 


R.K.  Soni  et  al.  /  Journal  of  Luminescence  83-84  (1999)  187-191 


191 


Raman  Shift  (cm'1) 


Fig.  5.  Calculated  resonance  contribution  to  the  experimentally 
observed  Raman  line  shape  of  1.4  nm  silicon  nanocrystal. 


interaction.  The  resonant  nanocrystals  at  the  laser 
photons  at  2.412  eV  have  diameters  in  the  range 
1.3-1.4nm.  Fig.  5  shows  the  calculated  line  shape 
for  1.4  nm  nanocrystal.  It  is  evident  form  the 
figure  that  the  main  contribution  to  the  Raman  line 
shape  comes  from  the  resonant  radii,  though  non¬ 
resonant  nanocrystals  provide  broad  constant 
background  to  the  Raman  spectrum. 


4.  Conclusions 

We  have  synthesized  Si  nanocrystals  in  a  Si02 
matrix  on  quartz  strip  by  RF  co-sputtering.  The 
silicon  nanocrystals  are  dispersed  in  size  along  the 
length  of  the  strip.  Optical  transmission  measure¬ 
ments  were  employed  to  estimate  the  nanocrystal 
size  distribution  which  shows  peaks  in  the  range 
1.1-2.6  nm  with  an  asymmetrical  tail  towards  the 
large  particle  size.  The  PL  emission  from  the 
Si-Si02  film  exhibits  broadband  in  the  red  and 
green  region,  and  their  energies  are  in  agreement 
with  corrected  LDA  calculation  of  Si  nanocrystals. 


Size  dependence  of  the  Raman  scattering  was 
performed  on  Si  nanocrystals,  the  results  show 
phonon  softening  and  large  asymmetrical  broaden¬ 
ing.  With  decreasing  nanocrystals  size,  broadening 
increases,  further  due  to  fluctuation  in  effective 
wave  vector  associated  with  size  distribution.  The 
observed  Raman  frequency  and  asymmetry  was 
calculated  by  considering  phonon  confinement  in 
a  spherical  nanocrystal.  For  small  size,  the  line 
shape  is  dominated  by  the  resonant  nanocrystals 
that  favor  incoming  or  outgoing  resonance.  The 
electronic  state  responsible  for  the  resonant  inter¬ 
action  with  nanocrystal  has  a  large  lifetime 
broadening  of  the  order  of  25  meV.  The  weak  res¬ 
onance  enhancement  seen  under  the  resonant  con¬ 
dition  is  indicative  of  band-to-band  transitions  in 
the  Raman  process. 
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Abstract 

Optical  quality,  free  standing  and  highly  photoluminescent  PSi  was  prepared  from  both,  p-  and  n-type  Si  wafers  using 
electrochemical  anodization  of  c-Si.  A  broad  and  intense  absorption  band  is  in  the  UV-Visible  region  with  a  cut-off  edge 
at  ~  400  nm.  A  strong  emission  was  found  in  the  red  region.  At  8  K  the  broadband  luminescence  is  peaked  at  ^  670  nm. 
The  band  maximum  shifts  toward  low-energy  side  while  excitation  wavelength  increases.  The  size  distribution  of  the 
nanoparticles  was  obtained  from  the  optical  transmission  data  and  the  theoretical  relation  between  the  size  of  the 
particles  and  the  energy  gap.  Nonlinear  optical  response  measurement  was  also  performed  using  degenerate  four-wave¬ 
mixing  (DFWM)  in  a  backward  configuration  with  all  waves  in  s-polarization.  The  response  signal  consists  of  an 
instantaneous  component  followed  by  a  long  lived,  slowly  decaying  component.  The  former  is  associated  with  the 
third-order  susceptibility  of  the  material  whereas  the  latter  originates  from  the  contribution  of  surface  states  created  by 
laser  excitation.  The  response  due  to  the  surface  state  is  significant,  which  has  potential  in  practical  applica¬ 
tions.  ©  1999  Elsevier  Science  B.V.  All  rights  reserved. 

PACS:  42.65.Re;  78.55.Mb;  78.55.  -  m 

Keywords:  Porous  silicon;  Photoluminescence;  Wave  mixing;  Optical  response 


1.  Introduction 

It  is  well  known  [1]  that  the  small  and  indirect 
band  gap  of  crystalline  silicon  (c-Si)  is  characterized 
by  an  extremely  low  quantum  efficiency  of  lumin¬ 
escence  while  porous  silicon  (PSi),  obtained  by 
electrochemical  etching  of  c-Si  surfaces,  yields  ex¬ 
tremely  high  photoluminescence  (PL)  efficiencies 
up  to  10" 1  in  the  visible  range  [2-5].  The  origin  of 
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the  PL  has  been  the  subject  of  many  studies.  We 
have  studied  the  photodynamical  process  and 
Raman  spectra  of  PSi  [4,6].  It  was  found  that 
S-band  luminescence  comes  from  highly  confined 
regions  with  a  distribution  of  size.  Quantum  con¬ 
finement  increases  the  emission  energy  and  the 
oscillator  strength  of  the  radiative  transition,  de¬ 
creasing  the  radiative  lifetime.  Therefore,  the  most 
accepted  luminescence  mechanism  at  this  time  is 
associated  with  quantum  confinement  in  the  nc-Si 
formed  in  the  PSi  by  the  anodization  process  [5], 
The  presence  of  ~  3  nm  size  nanocrystallites  is 
the  size  required  by  the  confinement  mechanisms 
to  enlarge  the  band  gap.  On  the  other  hand, 
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in  nonlinear  optical  experiments  a  key  parameter 
to  be  measured  is  the  nonlinear  refractive  index 
which  is  related  to  a  variety  of  other  NLO  effects 
through  the  formalism  of  NLO  susceptibilities.  To 
study  ultrafast  nonlinear  optical  response  on  the 
picosecond  time  scale  plays  an  important  role  for 
practical  NLO  applications  such  as  switching  and 
modulation.  However  no  study  on  PSi  has  been 
reported  to  our  knowledge.  In  this  paper  we  dem¬ 
onstrate  the  size-related  emission  property  and 
show  that  the  excitation  is  mainly  localized.  The 
localized  surface  state  is  responsible  for  the  post- 
instantaneous  signal. 

2.  Sample  preparation 

The  optical  quality,  free  standing  and  highly 
photoluminescent  PSi  was  obtained  from  both,  p- 
and  n-type  Si  wafers  using  electrochemical  anodiz¬ 
ation  of  c-Si.  The  key  role  of  the  holes  in  the  etching 
process  requires  that  n-type  samples  be  illuminated 
during  the  anodization.  The  porosity  and  the  final 
structural  and  optical  properties  of  the  samples  can 
be  varied  by  changing  the  acid  concentration,  cur¬ 
rent  density,  anodization  time  as  well  as  the  charac¬ 
teristics  of  the  Si-wafer  as  resistivity  and  crystalline 
orientation.  We  prepared  a  Si  wafer  with  (10  0) 
crystalline  orientation  and  0. 5— 1 . 5  £2  cm  resistivity. 
Using  a  HF :  ethanol :  water  (1:2:1)  solution,  the 
anodization  current  (200  mA/cm2)  was  applied  in 
280  steps  of  0.8  s  each  and  with  resting  intervals  of 
1.4  s  between  anodizations,  such  that  a  high-optical 
quality,  free-standing  porous  films  were  obtained. 
The  preliminary  studies  we  have  made  including 
the  determination  of  electronic  surface  states  dens¬ 
ities,  electrical  transport  properties,  optical  and 
morphological  properties  as  well  as  surface  proper¬ 
ties  show  that  the  material  prepared  is  homogene¬ 
ous  with  high  porosity.  The  passivation  of  the 
nanostructure  of  PSi  for  increasing  stability,  such 
as  post-oxidation  of  the  samples  was  also  proceed¬ 
ed  and  tested. 

The  obtained  free-standing  sample  is  transparent, 
brownish.  The  optical  absorption  measurement 
was  taken  using  a  Hitach  U-2001  spectro¬ 
photometer  at  room  temperature.  A  sharp  cut-off 
absorption  edge  is  located  at  410  nm.  The  absorp¬ 
tion  wing  extends  to  600  nm. 


3.  Crystallite  size  variation 

In  emission  measurement,  a  broad-band  red 
luminescence  can  generally  be  observed  by  illumina¬ 
tion  of  the  sample  at  wavelength  shorter  than 
600  nm.  Fig.  1  stands  for  typical  luminescence  spec¬ 
trum  of  PSi  at  room  and  low  temperature  at  8  K 
using  Ar+  laser  as  a  light  source.  At  room  temper¬ 
ature,  the  position  of  emission  band  maximum  is 
insensitive  to  the  excitation  wavelength  whereas  at 
low  temperature  of  8  K  the  peak  shifts  toward  low 
energy  in  parallel  to  the  excitation  photon  energy. 
Comparing  the  spectra  obtained  at  different  tem¬ 
peratures,  no  distinct  correlation  between  peak 
position  and  temperature,  or  between  emission 
bandwidth  and  temperature  was  observed.  It  sug¬ 
gests  that  all  excitations  are  highly  localized.  For 
excitation  at  458  nm,  the  electrons  in  conduction 
band  with  energy  well  above  the  band  gap  could 
undergo  a  direct  radiative  transition  in  the  area 
where  they  were  excited  while  the  others  could 
migrate  away.  The  excitation  migration  would  then 


Wavelength  ( nm ) 

Fig.  1.  Emission  spectra  of  PSi  by  excitation  of  an  Ar+  laser  at 
three  different  wave  lengths.  The  upper  graph  stands  for  low 
temperature  at  8  K,  while  the  lower  one  for  room  temperature. 
The  upper  inset  is  the  absorption  spectrum,  the  lower  inset 
represents  the  size  distribution  of  particles. 
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be  trapped  in  some  different  confined  region,  and 
emit  from  that  region.  If  the  latter  were  the  case 
the  entire  emission  would  be  an  overlap  of  entire 
emission  from  different  regions,  resulting  in  a  much 
broader  and  diffusive  feature  with  an  obvious 
wavelength  dependence.  In  fact,  in  our  sample  the 
633  nm  excitation  which  provides  energy  well  be¬ 
low  the  band  gap  also  exhibits  similar  emission  and 
bandwidth  as  well. 

Based  on  optical  absorbance  data  obtained  the 
size  distribution  of  the  silicon  particles  in  the  PSi 
sample  was  calculated.  The  method  is  as  follows. 
We  assume  that  the  sample  is  a  collection  of  spheri¬ 
cal  crystallites  whose  sizes  are  described  by  a  distri¬ 
bution  function,  F(r),  and  also  that  the  absorption 
is  due  to  the  interband  transition  of  the  electrons  in 
the  individual  crystallites. 

Due  to  the  quantum  confinement  of  the  free 
carriers  in  the  nanocrystallites  it  is  expected  that 
each  particle  is  transparent  for  energies  (£)  less  than 
its  optical  gap  (£g)  and  its  absorbance  increases 
abruptly  for  E  >  £g.  The  transmission  coefficient  of 
the  PSi  film  can  be  expressed  as 

rWocl-cj*  F(r)A(r,  A)  dr. 


where  C  is  a  constant  and  A(r,  A)  is  proportional  to 
the  energy  absorbed  by  a  particle  of  radius  r.  Air,  A) 
is  approximated  by  the  Heavy-side  step  function, 
A(r,  A)  oc  0(Ag,r  -  A)  where  Ag,r  =  2nhc/E g-r,  and  Eg>r 
is  the  energy  gap  for  particles  of  size  r.  From  the 
above  expression,  we  obtain  that  the  first  derivative 
of  the  transmission  spectrum  is  proportional  to 
the  size  distribution  density, 


F(r)S(r  —  r0)  dr  oc  -  F(r0), 


where  r0  is  the  size  which  gives  an  energy  gap 
corresponding  to  A0.  To  obtain  this  size  distribu¬ 
tion  we  used  the  relation  between  the  magnitude  of 
the  energy  gap  and  the  size  of  the  (spherical)  par¬ 
ticle  suggested  in  Ref.  [7],  £g-r  =  1.1  +  3.65  r~ 13 
where  £g>r  is  in  eV  and  r  in  nm.  As  shown  in  Fig.  1 
(lower  inset)  this  calculation  shows  a  nonsymmetric 
size  distribution  peaking  at  ~  1.65  nm  in  our  PSi 
film. 


4.  Picosecond  nonlinear  optical  response 

The  nonlinear  optical  response  measurement  in 
PSi  is  difficult  in  visible  to  UV  region  because  of 
higher  absorption  that  the  free-standing  sample  is 
used  to  be  burned  out  before  getting  signal.  In  this 
experiment  we  use  a  frequency-tripled  Y3A15012 : 
Nd  laser  operated  at  355  nm  with  pulse  width  of 
45  ps  as  a  pump  source  for  ps  optical  parametric 
generation  (OPG).  The  OPG  output  wavelength 
can  be  continuously  tuned  from  400  to  700  nm.  It 
was  found  that  significant  absorption  by  PSi  free 
standing  film  in  the  region  of  550-590  nm  would 
cause  a  damage  of  sample.  In  the  region  of 
590-600  nm  the  obtained  response  signal  is  domi¬ 
nant  by  the  slow-response  component  which  orig¬ 
inates  from  the  excited  state  population  grating 
formed  by  the  two  pump-pulses  across  inside  the 
sample.  Tuning  the  wavelength  beyond  610  nm 
caused  the  DFWM  signal  too  weak  to  record  in 
transverse  geometry.  The  final  measurement  was 
therefore,  conducted  at  604  nm.  In  the  experiment 
the  laser  beam  was  split  into  two  pump  beams  and 
a  probe  beam  as  shown  in  Fig.  2.  The  latter  carried 


Fig.  2.  Experimental  setup  of  picosecond  DFWM. 
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10%  of  the  total  laser  power.  The  probe  went 
through  an  optical  delay  line,  and  was  then  incident 
on  the  sample  from  the  substrate  side.  The  two 
pump  beams  were  crossed  on  the  film  at  an  angle, 
interfering  to  form  a  transient  grating  in  the  film. 
The  delay  line  in  one  of  the  pump  beam  paths  was 
used  to  ensure  simultaneity  of  the  pulses  at  the 
sample.  The  beam  size  of  the  probe  was  about  50% 
that  of  the  pump.  Thus,  the  probe  could  uniformly 
sense  the  susceptibility  within  the  interference  re¬ 
gion  of  the  two  pump  pulses.  DFWM  methods  for 
thin  films  have  been  amply  discussed  in  the  litera¬ 
ture  [8,9].  All  the  beams  in  our  experiment  were  in 
backward  phase  conjugation  configuration,  which 
can  allow  higher  signal-to-noise  ratio.  The  pump- 
pulse  intensity  was  typically  10-20  mJ/cm2.  The 
polarization  for  all  measurements  was  perpendi¬ 
cular  to  the  incidence  plane  for  all  four  beams, 
expressed  as  <ssss).  In  this  configuration,  the  time- 
resolved  NLO  spectrum  can  adequately  be  used  to 
characterize  the  amorphous  PSi’s  NLO  properties. 

Fig.  3  shows  the  DFWM  signal  of  free  standing 
PSi  wafer  compared  with  CS2  reference.  For  CS2 
reference,  the  intense  coherent  response  signal  with 
a  half-width  of  about  55  ps  located  at  zero  delay  of 
the  probe  pulse  is  associated  with  the  third-order 
nonlinearity  of  the  CS2  material.  It  is  an  instan¬ 
taneous  response.  It  comes  from  the  change  in 
polarizability  of  grating  peak  versus  grating  valley 
when  two  pump  pulses  are  cross  inside  the  sample 
writing  a  periodic  bright-dark  interference  pattern. 
For  CS2  reference,  the  fast  response  is  well  known 
as  due  to  the  vibrational  character  of  the  inter- 
molecular  motion  in  liquid.  The  signal  width  is 
consistent  with  the  intensity  autocorrelation  of  the 
three  laser  pulses.  For  PSi  sample,  however,  the 
signal  consists  of  two  components.  One  is  similar 
to  CS2,  with  the  peak  located  at  zero  delay  of 
the  probe  pulse.  The  other  is  a  slow  response 
signal,  can  be  identified  as  to  start  at  ~  +  30  ps  of 
the  probe  pulse  delay.  The  signal  intensity  of 
the  slow  response  (SR)  is  almost  a  half  of  the  instan¬ 
taneous  response  signal,  but  it  lives  for  much  lon¬ 
ger.  By  estimation  it  can  enter  into  microseconds 
region. 

The  first  component  at  zero  delay  is  a  direct 
measure  of  x(3)  of  the  PSi  material.  It  has  two 
distinct  features.  Firstly,  the  ratio  of  peak  intensity 
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Fig.  3.  Picosecond  nonlinear  optical  response  of  PSi  compared 
with  CS2  reference. 

to  sample  thickness  is  extremely  large  compared  to 
most  bulk  crystals.  In  solid-state  materials  the  con¬ 
tribution  to  the  instantaneous  signal  component 
may  come  from  electronic  cloud  deformation  as 
well  as  the  nucleus  reorientation.  The  latter  contri¬ 
bution  in  solids  is  usually  smaller.  The  second  fea¬ 
ture  of  the  coherent  response  component  is  the 
squeezed  nature  in  time  domain. 

When  the  two  pump  pulses  are  polarized  in  the 
same  direction,  the  coherent  signal  component  is 
directly  related  to  the  third-order  susceptibility  by 
rj  —  exp(—  ad/cos  6)s\n2(dnAn/2  cos  0),  where  An  = 
l2nx{3)/n0(E2}.  E  is  optical  electric  field,  d  is  the 
grating  thickness,  and  26  is  the  crossing  angle  of  the 
two  pump  pulses.  The  intrinsic  nonlinear  optical 
response  signal  generated  in  this  circumstance  can 
be  analyzed  as  the  probe  beam  diffracted  by  the 
grating  formed  by  the  two  pump  pulses  A  and  B, 
which  intersect  at  an  angle  of  26  inside  the  thin  film. 
The  grating  vector  is  q  -  ±(kA-  *B)>  confined  in 
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the  film  along  the  x  direction.  The  electric  field 
amplitude  A  is  A  —  AAQ+lkx  +  ABe~lkx ,  where  k  rep¬ 
resents  the  wave  vector.  The  light  intensity  is 
then  modulated  as  /  =  /A  +  /„  +  2A7  cos(2/cx), 
here  A 7  =  l/2ne0AA  •  AB.  In  our  experiment, 
/a  =  hy  so  /  =  2/a(1  +  cos {qx))  for  identical  polar¬ 
izations  (ss)  of  Aa  and  AB.  The  diffracted  signal 
intensity  is  directly  related  to  the  third-order  sus¬ 
ceptibility  as  expressed  by  Eq.  (1).  Compare  the 
FWHM  of  Psi  sample  to  CS2,  the  former  is  about 
70%  of  the  latter  compared  to  the  electronic 
part.  For  isotropic  sample  *£«  =  =  hyVxy 

[14]-  We  obtained  =  Xmy  =  =  0.12  x 

10  1 1  cm3  erg,  which  is  significantly  large  com- 
pared  to  many  bulk  crystals. 

The  “squeezing”  feature  observed  in  the  PSi 
wafer  is  quite  interesting.  We  previously  reported 
a  compression  effect  of  the  DFWM  signal  in  bulk 
KNb03  crystal.  It  was  due  to  TP  A  interaction  [10] 
between  the  pump  and  probe  pulses  when  the  laser 
fluence  was  greater  than  280mJ/cm2.  The  nar¬ 
rowed  signal  was  shifted  towards  negative  delay. 
In  the  present  case,  however,  the  laser  fluence  was 
lower  and  the  peak  does  not  shift.  Further  increas¬ 
ing  the  laser  power  would,  unfortunately,  cause  the 
sample  to  be  damaged.  For  other  thin  amorphous 
samples  such  as  amorphous  quartz  wafer,  no 
squeezing  effect  could  be  observed.  The  reason  for 
this  is  not  completely  understood  yet.  In  four- 
wave-mixing,  squeezed  light  was  first  observed  by 
Slusher  et  al.  [11]  in  an  optical  cavity  (atomic 
vapor).  It  was  due  to  the  generation  of  “squeezed 
states”.  In  this  phenomenon  the  statistical  variance 
in  one  of  the  two  in-quadrature  components  of  the 
coherent  electric  field  is  reduced  below  the  shot 
noise  value,  while  variance  of  the  other  component 
is  increased  so  that  the  uncertainty  principle  is 
satisfied.  The  squeezing,  which  could  have  applica¬ 
tions  in  communications,  has  also  been  observed 
in  semiconductors,  optical  fibers,  as  well  as  in 
waveguides.  It  is  due  in  these  cases  to  a  nonlinear 
mechanism  through  the  third-order  nonlinear  sus¬ 
ceptibility,  including  Kerr  effect,  which  can  affect 
the  transverse  spatial  profile  of  the  generated  light 
[12,13].  In  our  case  it  possibly  relates  to  a  similar 
mechanism,  except  for  the  small  optical  interaction 
length  in  the  film. 


The  second  component  in  Fig.  3  is  a  SR  signal. 
This  signal  component  is  attributed  to  the  popula¬ 
tion  grating  produced  by  the  interference  of  the 
two  crossed  laser  pulses.  In  the  bright  region  of  the 
grating,  the  intrinsic  interband  transition  occurs, 
forming  highly  localized  surface  state  as  discussed 
in  the  above  section.  It  in  turn  alters  the  polariza¬ 
bility.  In  many  crystal  samples,  the  SR  signal,  fol¬ 
lowing -the  coherent  component,  used  to  be  built 
up  for  a  risetime  of  up  to  10  ps  [10],  which  corres¬ 
ponds  to  relaxation  process  of  the  excited  state 
dynamics.  In  PSi  sample,  however,  no  risetime  was 
observed,  which  further  proves  that  the  excitation 
at  604  nm  was  highly  localized  without  further 
migration  or  other  dynamical  processes. 
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Abtrasct  Optical  quality,  free  standing  and  highly  photoluminescent  PSi  was  prepared  from  both,  p- 
and  n-type  Si  wafers  using  electrochemical  anodization  of  c-Si.  The  porosity  and  the  final  structural 
and  optical  properties  of  the  samples  were  found  to  be  altered  drastically  by  changing  the  acid 
concentration,  current  density,  anodization  time  as  well  as  the  characteristics  of  the  Si-wafer.  In 
optical  measurement,  a  broad  and  intense  absorption  band  is  in  the  UV-Visible  region  with  a  cut-off 
edge  at  -  400  nm.  A  strong  red  emission  was  found  to  be  at  ~  670  nm. .  It  is  extremely  broad  even  at 
8  K  The  experimental  results  show  that  a  group  of  sites  can  be  simultaneously  excited  by  pump  laser 
due  to  excitation  energy  coincidence.  In  nonlinear  optical  measurement  the  optical  response  was 
characterized  using  degenerate  four-wave-mixing(DFWM)  in  a  backward  configuration.  The  response 
signal  consists  of  an  instantaneous  component  followed  by  a  long  lived,  slowly  decaying  component. 
They  are  associated  with  the  third-order  susceptibility  of  PSi  and  the  surface  states.  Tlie  response  due 
to  the  surface  state  is  significant. 


INTRODUCTION 

It  is  well  known  that  the  small  and  indirect  band  gap  of  crystalline  silicon  (c-Si)  is 
characterized  by  an  extremely  low  quantum  efficiency  of  luminescence  while  porous  silicon 
(PSi),  obtained  by  electrochemical  etching  of  c-Si  surfaces,  yields  extremely  high 
photoluminescence  (PL)  efficiencies  up  to  10'1  in  the  visible  range  (2).  The  origin  of  the 
PL  has  been  the  subject  of  many  studies.  We  have  studied  the  photodynamical  process  and 
Raman  spectra  of  PSi  (3,4).  It  was  found  that  S-band  luminescence  comes  from  highly 
confined  regions  with  a  distribution  of  wire  widths.  Quantum  confinement  increases  the 
emission  energy  and  the  oscillator  strength  of  the  radiative  transition,  decreasing  the 
radiative  lifetime.  Therefore  the  most  accepted  luminescence  mechanism  at  this  time  is 
associated  with  quantum  confinement  in  the  nc-Si  formed  in  the  PSi  by  the  anodization 
process.  The  presence  of  -3  nm  size  nanocrystallites  is  the  size  required  by  the 
confinement  mechanisms  to  enlarge  the  band  gap.  On  the  other  hand,  there  is  interest  in  the 
changes  of  the  size-sensitive  optical  properties  due  to  the  quantum  confinement.  The 
microstructure  of  the  vitreous  state  exhibits  short-range  order  within  a  c-Si  core  but 
long-range  disorder.  The  wave  vector  and  eigenfunction  associated  with  an  optical  process 
are  thus  no  longer  valid  due  to  the  lack  of  periodicity  of  the  structure.  In  PSi  the  network 
can  be  tailored  within  one  size-category  by  using  size-sensitive,  fast  selective  excitation. 
This  will  reveal  the  associated  optical  properties  of  the  material,  such  as  photorefractive 
sensitivity,  FWM  signal  polarization,  appearance  of  ultrafast  NLO  responses,  optical 
activity,  etc. 
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In  nonlinear  optical  experiments  a  key  parameter  to  be  measured  is  the  nonlinear 
refractive  index  which  is  related  to  a  variety  of  other  NLO  effects  through  the  formalism 
of  NLO  susceptibilities.  It  is  through  these  related  effects  (polarization  rotation,  optical 
Kerr  effect,  self-phase  modulation,  optical  mixing,  etc.)  that  much  of  the  available  data  on 
nonlinear  refraction  as  well  as  the  fast  third-order  nonlinear  properties  of  the  material  can 
be  obtained.  The  ultrafast  response  time  on  the  picosecond  time  scale  is  in  fact  a  key  point 
for  practical  NLO  applications  such  as  switching  and  modulation.  However  no  study  on 
the  NLO  response  from  PSi  has  been  reported  to  our  knowledge.  In  this  paper  we  present 
the  experimental  results  on  degenerate-four-wave-mixing  (DFWM)  measummet  as  well  as 
the  conventional  optical  spectral. 


EXPERIMENTAL 

The  optical  quality,  free  standing  and  highly  photoluminescent  PSi  was  obtained 
from  both,  p-  and  n-type  Si  wafers  using  electrochemical  anodization  of  c-Si.  The  acid 
solution  for  the  anodic  etching  was  1:2:1  (HF: ethanol: water).  The  electrolytic  cell  included 
a  standard  Platinum  electrode  as  a  cathode  and  the  crystalline  silicon  wafer  as  an  anode.  The 
aluminum  ohmic  contact  at  the  bottom  of  the  wafers  were  made  by  sputtering  deposition 
and  annealing.  Prior  to  the  anodization,  the  silicon  surface  of  the  wafers  were  cleaned  with 
20%  HF-aqueous  solution  for  5  minutes.  Samples  were  anodized  and  finally  dried  using 
the  critical  point  method.  Free  standing  samples  were  also  prepared  by  applying  a  high 
electropolishing  pulse  at  the  end  of  the  anodizing  process.  Such  a  pulse  disengages  the 
porous  layer  from  the  underlying  crystalline  silicon.  The  key  role  of  the  holes  in  the  etching 
process  requires  the  illumination  of  the  n-type  samples  during  the  anodization.  The  porosity 
and  the  final  structural  and  optical  properties  of  the  samples  can  be  varied  by  changing  the 
acid  concentration,  current  density,  anodization  time  as  well  as  the  characteristics  of  the 
Si- wafer  (resistivity  and  crystalline  orientation).  The  preliminary  studies  we  have  made 
including  the  determination  of  electronic  surface  states  densities,  electrical  transport 
properties,  and  optical  and  morphological  properties  and  surface  properties  show  that  the 
material  prepared  is  homogeneous  with  high  porosity.  The  passivation  of  the  nanostructure 
of  PSi  for  increasing  stability,  such  as  post-oxidation  of  the  samples  was  also  proceeded  and 
tested.  Optically,  nanoporous-Si  behaves  as  a  homogenous  material  in  the  visible  and  longer 
wavelength  regions,  as  tested  by  ellipsometry,  reflectivity  measurements  and  by 
transmitance  as  well.  The  sample  studied  here  was  prepared  from  a  boron  doped  Si  wafer 
with  (100)  crystalline  orientation  and  0.5-1. 5  Qcm  resistivity.  The  anodization  current 
(200mA/cm2)  was  applied  in  280  steps  of  0.8sec  each  and  with  resting  intervals  of  1.4  sec 
between  anodizations,  such  that  a  high-optical  quality,  free-standing  porous  film  was 
obtained. 
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Wavelength  ( nm ) 


FIGURE  1.  Emission  spectra  of  PSi  by  excitation 
of  an  Ar+  laser  at  3  different  wavelengths.  The  up¬ 
per  graph  stands  for  low  temperature  at  8  K,  while 
the  lower  one  for  room  temperature.  The  inset  rep¬ 
resen  represents  the  absorption  spectrum. 


FIGURE  2.  Time-resolved  emissionspectra  of 
PSi  excited  at  532  nm.  The  delay  time  from  up¬ 
per  to  lower:4,  24,  44,  64,  84  ps. 


Absorption  and  emission 

The  absorption  measurement  was  conducted  using  a  Hitach  spectrophotometer  at  room 
temperature.  A  sharp  cut-off  absorption  edge  at  the  high  energy  side  is  located  at  -  410nm 
and  the  low  energy  side  is  at  -  5 1 0  nm  (shown  in  the  inset  of  Figure  1)  with  a  tail  extended 
to  600  nm.  At  the  wavelength  greater  than  5 1 0  nm  the  absorption  becomes  very  weak,  and 
at  633  nmn  the  absorption  is  negligible. 

The  emission  measurements  were  conducted  using  lasers  setting  at  different  wavelengths 
covering  458  -  633  nm  region  in  purpose  to  provide  different  excitation  photon  energy. 
Nevertheless,  in  all  cases  a  broad  band  red  luminescence  can  generally  be  observed.  Fig.  1 
stands  for  a  typical  luminescence  spectra  of  PSi  at  room  and  low  temperature  8  K.  Three 
lines  from  an  Ar+  laser  at  458,  488  and  514  nm  were  used  for  excitation.  At  room 
temperature,  the  emission  band  maximum  does  not  shift  markedly  with  changing  excitation 
wavelength,  but  emission  intensity  varies  significantly  because  absorption  changes  from  458 
nm  to  5 1 4  nm.  At  low  temperature  8  K  however,  the  emission  peak  obviously  shifts  toward 
longer  wavelengths  while  the  excitation  tuned  to  longer  wavelength.  But  the  emission 
intensity  does  not  change  so  much  as  what  observed  at  room  temperature.  It  is  presumably 
due  to  a  flat  plateau  in  this  absorption  region  at  low  temperature.  On  the  other  hand,  to 
compare  the  spectra  obtained  at  different  temperature  no  distinct  correlation  between  peak 
position  and  temperature,  or  between  emission  bandwidth  and  temperature  was  observed. 
Furthermore,  emission  spectra  excited  at  633  nm.were  also  obtained,  and  used  to  compare 
with  those  excited  at  shorter  wavelengths.  The  photon  energy  of  633  nm  excitation  is 
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51  ^  ^  absorptionfed8e  51?  However,  the  luminescence  can  still  be  observed 
all  cases  the  common  feature  for  all  emission  spectra  is  the  large  bandwidth  From 

Therefore  V  *°  6  ?  T’  T  difference  bandwidth  was  found  to  be  less  than  10% 

Therefore,  it  is  suggested  that  all  excitations  are  highly  localized.  Fig.2  shows  the  time- 

resolved  enussion  spectra  at  room  temperature.  The  excitation  was  at  532  nm  with  the 
delay  time  from  4  to  84gs.  It  features  a  typically  glassy-like  profile.  As  the  delay  time 
increases  shorter-lived  emitter  dies  out  leaving  luminescence  due  to  longer-lived  e^tter 

Tk”  6  figUre  *Jhere  *S  n° inleraction  evidence  in  any  Emission  profile 
of 3d oM  STf Jf honon-assistedtransition  is  less  important  here,  compared  to  the  case 

Rom  aln/  T  h  ^  T  *S  1116  reaSOn  is  two-fold>  35  discussed  in  the  following 
JL  °™  ET-  U“'ng  T  may  assume  two  Possible  recombination  processes  For 
those  electrons  excited  from  valence  band  to  conduction  band  by  excitation  at  458  nm  they 

possess  excitation  energy  well  above  the  bandgap.  A  part  of  them  could  undergoTdirecit 
adiative  transition  m  the  area  where  they  were  excited  while  the  others  could  migrate  away 

fromThat  r^'1  T?!L°n  ,W°  thCn  be  traPped  m  some  different  “^ned  region  and  emit 
fnf^  ba  eg'°"  ,atter  case  were  true  the  entire  emission  would  be  an  overlap  of 

entire  emission  from  different  regions,  resulting  in  a  much  broader  and  diffusive  feature  £ith 
an  obvious  wavelength  dependence.  On  the  other  hand,  for  thoseSro^Sed t 
633nm  excitation  they  possess  energy  well  below  the  bandgap.  These  electrons  can  not 
foe r8°t  T 8131,0 T  PuCeSS  They  must  '"""ediately  be  localized  and/or  emit  hgTfrom 

L  SI™  tE18  T  emiSSi°n  Pr0file  is  suPP°sed  to  be  ^t  as  broaden*  foat 
...  *ii,  Jherefore>  because  all  emissions  excited  at  different  wavelengths 

^  H  “  "”Pfe  T  'he  light  emissio"! fSJ 

'”„d  “d  s?ttered  potots  Si  samples.  This  is  in  consistence  with  conclusions  of 

Si  havejstmn  nZ  !Z  Studies  “*  polariz!“io"  °f  the  luminescence  of  porous 

WoZTTlZ: Z1 ^  lT,ed  m‘Sra,i0n  °f  <°  "<%hbon„8  SI 

J  ...  ““"d-fold  consideration  is  that  in  nanocrystais,  the  phonon  snectra 

vibration  may  also  become  important.  A  large  number  of  activators  are  exDected  to  he 

he  locTlLr^f latt,Ce  C,nvir0nment  near  the  surface  The  radiative  recombination  in 
the  localized  state  on  the  nanociystallite  surface  causes  the  strong  vicihie 

to  Se '.Ue"ee  ^  “  °f  PL  “  ^  -  USZ 


Transient  optical  response  from  PSi 

recnSr!Udy  0"  °ptical  proPerty  of  P°™is  Si,  particularly  the  transient  ontical 

2ES  rr  b?  n0t  yet  rep0rted  t0  0ur  knowledge,  ft  is  duHSy to  the 
fovS  fffiCUty  10  perfonnin8  such  ^d  of  measurements  With  a  green  light  source 

! pted  lasei? the  absorption  by  PSi 

picosecond  lipht  «  .  6  a^e  ab  e  t0  Pefform  this  kind  of  measurement  by  turning  the 

JltiT*  r  th?  yeI,0W  -  ora"Se  "*»■  using  technique  of  picosecond 
ptical  parametric  generation.  A  frequency-tripled  Y3Al5012:Nd  laser  operated  at  355  nm 
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FIGURE  3.  Picosecond  nonlinear  optical 
response  of  PSi  compared  with  CS^  reference. 


- _ -  with  pulse  width  of  45  ps  was  used  for 

12  -  J  optical  parametric  generation.  The  output 

Hi  cs2  wavelength  can  be  continuously  tuned 

a  -  /  I  from  400  to  700  nm.  In  order  to 

^  l  |  determine  which  wavelength  can  be  best 

|  4  -  I  I  used  for  characterization,  the 

.a  j  l  measurement  was  checked  by  turning 

i  o  -  laser  was  within  550  -  650  nm.  Significant 

&  absorption  by  Psi  free  standing  film  in  the 

J  6  -  I  region  of  550  -  590  nm  would  cause  a 

fl  PSI  damage  of  sample.  In  the  region  of  590  - 

•i=  4  “  I  ul  iuM  I  ajl  |  ij|]|  600  nm  the  obtained  response  signal  is 

dominant  .by  the  slow-response 
Q  2 "  .  j  *  component  which  originates  from  the 

0  _  excited  state  population  grating  formed  by 

1  7  7  1  the  two  pump-pulses  across  inside  the 

-200  -100  0  100  200  300  400  *  L  . , 

n  ^  ,  ,  v  sample.  Tuning  the  wavelength  beyond 

Probe  Pulse  Delay  ( psec )  ,,  A  7  *  i  * 

610  nm  caused  the  DFWM  signal  too 

weak  to  record  in  transverse  geometry. 

FIGURE  3.  Picosecond  nonlinear  optical  The  final  measurement  was  therefore 

response  of  PSi  compared  with  CS^  reference.  conducted  at  604  nm.  In  the  experiment 

the  laser  beam  was  split  into  two  pump 

beams  and  a  probe  beam.  The  latter  carried  10%  of  the  total  laser  power.  The  probe  went 

through  an  optical  delay  line,  and  was  then  incident  on  the  sample  from  the  substrate  side. 

The  two  pump  beams  were  crossed  on  the  film  at  an  angle,  interfering  to  form  a  transient 

grating  in  the  film.  The  delay  line  in  one  of  the  pump  beam  paths  was  used  to  ensure 

simultaneity  of  the  pulses  at  the  sample.  The  beam  size  of  the  probe  was  about  50%  that 

of  the  pump.  Thus,  the  probe  could  uniformly  sense  the  susceptibility  within  the  interference 

region  of  the  two  pump  pulses.  DFWM  methods  for  thin  films  have  been  amply  discussed 

in  the  literature  (8,9).  All  the  beams  in  our  experiment  were  in  backward  phase 

conjugation  configuration,  which  can  allow  higher  signal-to-noise  ratio.  The  pump-pulse 

intensity  was  typically  10-20  mJ/cm2.  The  polarization  for  all  measurements  was 

perpendicular  to  the  incidence  plane  for  all  4  beams,  expressed  as  <ssss>.  In  this 

configuration,  the  time-resolved  NLO  spectrum  can  adequately  be  used  to  characterize  the 

amorphous  PSi’s  NLO  properties. 

Fig.  3  shows  the  DFWM  signal  of  free  standing  PSi  wafer  compared  with  CS2  reference. 
For  CS2  reference,  the  intense  coherent  response  signal  with  a  half  width  of  about  55  psec 
located  at  zero  delay  of  the  probe  pulse  is  associated  with  the  third-order  nonlinearity  of  the 
CS2  material.  It  is  an  instantaneous  response.  It  comes  from  the  change  in  polarizability 
of  grating  peak  versus  grating  valley  when  2  pump  pulses  are  cross  inside  the  sample 
writing  a  periodic  bright-dark  interference  pattern.  For  CS2  reference,  the  fast  response  is 
well  know  as  due  to  the  vibrational  character  of  the  intermolecular  motion  in  liquid.  The 
signal  width  is  consistent  with  the  intensity  autocorrelation  of  the  three  laser  pulses.  For 
PSi  sample,  however,  the  signal  consists  of  2  components.  The  one  is  similar  to  CS2,  with 
the  peak  located  at  zero  delay  of  the  probe  pulse.  The  other  one  is  a  slow  response  signal, 
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can  be  identified  as  to  start  at  ~  +30  ps  of  the  probe  pulse  delay.  The  signal  intensity  of  the 
slow  response  (SR)  is  almost  a  half  of  the  instantaneous  response  signal  but  it  lives  for 
much  longer.  By  estimation  it  can  enter  into  microseconds  region. 

The  first  component  at  zero  delay  is  a  direct  measure  of  x(3)  of  the  PSi  material.  It  has 
two  distinct  features.  Firstly,  the  ratio  of  peak  intensity  to  sample  thickness  is  extremely 
large  compared  to  most  bulk  crystals.  In  solid  state  materials  the  contribution  to  the 
instantaneous  signal  component  may  come  from  electronic  cloud  deformation  as  well  as 
the  nucleus  reorientation.  The  latter  contribution  in  solids  is  usually  smaller  The  second 
leat  ure  of  the  coherent  response  component  is  the  squeezed  nature  in  time  domain 
When  the  two  pump  pulses  are  polarized  in  the  same  direction,  the  coherent  signal 
component  is  directly  related  to  the  third-order  susceptibility  by 


where 


rj  -  exp(- ad/cos 0)  sv^idnAn/Acosd), 

a  12  n  ,,, 

n„<E  >  *  ’ 


here  £  is  optical  electric  field,  d  is  the  grating  thickness,  and  2d  is  the  crossing  angle  of 
e  o  pump  pu  ses.  The  intrinsic  nonlinear  optical  response  signal  generated  in  this 
circumstance  can  be  analyzed  as  the  probe  beam  diffracted  by  the  grating  formed  by  the  two 
pump  pulses  A  and  B,  which  intersect  at  an  angle  of  20  inside  the  thin  film.  The  grating 

vector  isq  ±(kA-kB ),  confined  in  the  film  along  the  x  direction.  The  electric  field 
amplitude  A  is 


A-A,,,'*-  ,  (3) 

where  k  represents  the  wavevector.  The  light  intensity  is  then  modulated  as 

/  =  /4  +  4  +  2AI  cos(2kx)  ,  (4) 

here  M=l/2ne^A-Ah.  In  our  experiment,  /,=/*  so  I=2IA(l+cos(qx))  for  identical 
P,°  anza*lons  (ss)  of. Aa  and  AB.  The  diffracted  signal  intensity  is  directly  related  to  the 
third-order  susceptibility  as  expressed  by  Eq.(l).  Compare  the  FWHM  of  PSi  sample  to 

CSi  ;he(S,rm!r's  Sout  7°0/»ofthe  latter  compared  to  the  electronic  part.  For  isotropic 
“ft*  “  .  3  X.  3  X  >**(10).  We  obtained  x'3^  =  x'3^  =  X13^  =  0. 12  x  KT11 
cm  erg,  which  is  significantly  large  compared  to  many  bulk  crystals 

The  “squeezing”  feature  observed  in  the  PSi  wafer  is  quite  interesting.  We  previously 
reported  a  compression  effect  of  the  DFWM  signal  in  bulk  KNb03  crystal.  It  was  due  to 
TPA  interaction  (11)  between  the  pump  and  probe  pulses  when  the  laser  fluence  was 
greater  than  280  mJ/cm  The  narrowed  signal  was  shifted  towards  negative  delay.  In  the 
present  case,  however,  the  laser  fluence  was  lower  and  the  peak  does  not  shift.  Further 
increasing  the  laser  power  would,  unfortunately,  cause  the  sample  to  be  damaged.  For 
ot  er  t  n  amorphous  samples  such  as  amorphous  quartz  wafer,  no  squeezing  effect  could 
e  observed.  The  reason  for  this  is  not  completely  understood  yet.  In  four-wave-mixing. 
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squeezed  light  was  first  observed  by  R  E.  Slusher  et  al  (12)  in  an  optical  cavity  (atomic 
vapor).  It  was  due  to  the  generation  of  “squeezed  states”.  In  this  phenomenon  the 
statistical  variance  in  one  of  the  two  in-quadrature  components  of  the  coherent  electric  field 
is  reduced  below  the  shot  noise  value,  while  variance  of  the  other  component  is  increased 
so  that  the  uncertainty  principle  is  satisfied.  The  squeezing,  which  could  have  applications 
in  communications,  has  also  been  observed  in  semiconductors,  optical  fibers,  as  well  as  in 
waveguides.  It  is  due  in  these  cases  to  a  nonlinear  mechanism  through  the  third-order 
nonlinear  susceptibility,  including  Kerr  effect,  which  can  affect  the  transverse  spatial  profile 
of  the  generated  light  (13,14).  In  our  case  it  possibly  relates  to  a  similar  mechanism,  except 
for  the  small  optical  interaction  length  in  the  film. 

The  second  component  in  Fig.3  is  a  SR  signal.  This  signal  component  is  attributed  to  the 
population  grating  produced  by  the  interference  of  the  two  crossed  laser  pulses.  In  the 
bright  region  of  the  grating,  the  intrinsic  interband  transition  occurs,  forming  highly 
localized  surface  state  as  discussed  in  the  above  section.  It  in  turn  alters  the  polarizability. 
In  many  crystal  samples,  the  SR  signal,  following  the  coherent  component,  used  to  be  built 
up  for  a  risetime  of  up  to  tens  picoseconds  (14),  which  corresponds  to  relaxation  process 
of  the  excited  state  dynamics.  In  PSi  sample,  however,  no  risetime  was  observed,  which 
further  proves  that  the  excitation  at  604  nm  was  highly  localized  without  further  migration 
or  other  dynamical  processes. 
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ABSTRACT 

We  have  carried  out  a  resonance  Raman  study  of  line-shape  in  silicon  quantum  dots 
synthesized  on  a  quartz  substrate  by  co-sputtering  bulk  Si  and  Si02.  Optical  transmission 
measurements  are  used  to  evaluate  dot  size  distribution.  The  size  distribution  shows  peaks 
around  1.0  and  1.4  nm.  The  Si  dots  exhibit  photoluminescence  in  the  visible  region,  which 
shifts  to  higher  energy  with  decreasing  size.  The  size  dependence  of  Raman  scattering  shows 
phonon  softening  and  increasing  asymmetrical  broadening  for  small  dots  (<  2nm).  The 
observed  spectra  are  compared  with  calculations  considering  electron-hole  interactions  at  a 
quasi-direct  gap  of  a  spherical  quantum  dot. 

INTRODUCTION 

Nanometer  size  Si  crystals  show  marked  improvement  in  photoluminescence  (PL) 
quantum  efficiency  compared  to  bulk  silicon  and  are,  therefore,  attractive  for  silicon  based 
optoelectronic  applications  [1].  The  PL  observed  from  Si  nanocrystals  is  attributed  to 
quantum  size  effects.  The  effect  of  carrier  confinement  and  band  gap  upshift  in  Si 
nanocrystals  has  been  discussed  intensively  in  the  literature  [1].  The  photoluminescence 
energy  for  Si  nanocrystals  have  been  calculated  by  using  the  empirical  psudopotentials  [2], 
the  third  nearest  neighbor  tight  binding  [3]  and  time  dependent  tight  binding  techniques  [4]. 
Despite  the  apparent  disagreement  in  predicting  the  peak  PL  energy,  these  calculations 
demonstrate  that  the  PL  energy  increases  as  the  nanocrystal  size  decreases  and  the  dominant 
contribution  to  the  visible  light  emission  comes  from  nanocrystals  smaller  than  2  nm. 

The  resonantly  excited  PL  spectrum  has  shown  long  radiative  lifetimes  in  Si 
nanocrystals  and  provides  important  evidence  that  phonons  are  directly  involved  in  the 
radiative  recombination  process  [5].  Clearly,  this  suggests  that  the  Si  nanocrystal  energy  band 
remains  indirect  type.  Theoretical  calculations,  however,  indicate  that  a  quasi-direct  gap  in 
nanocrystal  Si  is  possible  for  sizes  1 .0- 1 .5  nm  in  diameter  [3]. 

Raman  scattering  techniques  are  commonly  used  for  semiquantitative  determination  of 
size  effects  on  vibrational  modes,  such  as  confined  optical  phonon,  surface  phonon  and 
confined  acoustic  phonon,  in  nanocrystal  or  quantum  dots  (QD)  [6-11].  An  adequate 
knowledge  of  the  bulk  phonon  dispersion  is  a  prerequisite  for  understanding  lattice 
dynamical  properties  and  electron-phonon  interactions  in  QD.  The  first-order  resonant 
Raman  scattering  via  the  Frolich  interaction  has  been  investigated  for  spherical  QD  [9].  The 
electronic  intermediate  state  was  considered  as  uncorrelated  electron-hole  pairs.  Electron- 
hole  correlation  effects  via  Coulomb  interaction  on  Raman  amplitude  for  spherical  QD  has 
also  been  reported  [10,11].  It  was  shown  that  the  resonant  scattering  efficiencies  depend 
strongly  on  the  Coloumb  correlation  and  only  excitonic  states  and  vibrational  modes  with 
zero  angular  momentum  were  allowed  in  the  Raman  process. 

The  effect  of  carrier  confinement  on  optical  phonons  in  Si  QD  has  been  widely  studied 
by  Raman  scattering  [6-8].  The  size  effect  is  reflected  in  a  finite  wave  vector  of  QD  and 
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So11!r:rrftPhr0n  frequency  fr0m  bu,k-  A  quantum  confinement  model  that 
takes  into  account  the  finite  phonon  wave  vector  and  gaussian  size  distribution  for  spherical 

quantum  dots  is  generally  invoked  to  explain  the  observed  Raman  line  shape  In  thfrreport 
we  investigate  the  Raman  line  shape  from  Si  dots  dispersed  in  a  Si02  matrix  grown  on  a 

iLcominfnr^'i'  Tphasiz\the  role  of  Si  dots  that  favor  the  resonant  interaction  with 
incoming  or  outgoing  photon  in  the  Raman  process. 

EXPERIMENTAL 

_  Th/  Si;Si°2  f,lm  Wasgrown  on  the  substrate  by  rf  co-sputtering.  The  substrate 
Z  J  6'f!n.ch  i°ng  quaif  P'a,e- 0n  a  Si°2  ^get  (6  inch  in  diameter)  small  pieces  of  Si  tips 

?ndth’p,VnnT>  ^dand  ,hey  Were  c°-sputtered  for  3  hours.  The  number  of  Si  tips 
and  their  position  on  the  Si02  were  used  to  control  the  size  of  silicon  dots.  The  film  was 

subsequently  annealed  in  nitrogen  ambient  at  1 100°C  for  30  minutes.  The  deposited  film  has 

vaiTmnnnt  fr  ^  ‘^'T5  °f  lhe  °rder  °f  17  ^  1)01  size  and  *«icon  concentration 

5riZ  !!tn  f  ,  fh  0f  23he  71,6  estimated  Si  dot  density  varies  from 
1.5x10  /m  at  S.02-nch  end  to  3.5xlOu/m3  at  Si-rich  end  of  the  film.  The  transmission 

measurements  were  made  in  the  range  190-820  nm  using  a  HP-8452A  spectrophotometer 

The  photoluminescence  and  Raman  measurements  were  performed  in  a  quasi  back-scatterine 

T8  11,6  5145  nm  Hne  0f  a"  Ar+  ,aser-  ^'atpo“8 

20  ^  and  he  aser  beam  was  focussed  to  a  spot  of  25  pm  size  on  the  film  The 

320  ~  -  —  *  * 


Wavelength  (nm) 


Fig.  1(a).  Transmission  spectra  from  two 
regions  (A  and  B)  of  Si-Si02 
film. 

RESULTS 


Dot  diameter  (nm) 


Fig.  1  (b).  Size  distribution  at  the  region  A 
and  B  of  the  film  calculated 
from  transmission  spectra. 


Figure  1(a)  shows  typical  transmission  spectra  taken  from  two  regions  of  Si-SiO,  film 
to  our  case  the  Si  dot  size  is  dispersed  along  the  length  of  the  film  on  Z  ouaiS  p£e  te 

sfre  R  rmbedded^he^O116  ^  t0  ^  abS°rbed  by  ,he  Si  ^antum  *** 

R  imbedded  in  the  Si02  matnx.  From  the  optical  transmission  we  calculated  the  dot  size 
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.  ,7_!g  dependence  of  energy  gap  on  the  size  of  spherical  dot  as  1  l-nfis 
K  [*2],  size  distribution  can  be  approximated  with  the  , J  ,  °otas  11  +3-65 

the  distributicm^comesbroader!  "m  “  “  «■"»  With  --easing  sfze. 


1  40  1.60  1.80  2.00  2.20  2.40 
Energy  (eV) 

Fig.  2.  Normalized  PL  spectra  from  (a) 
1.4  nm,  and  (b)  2.2  nm  diameter  Si 
dots  excited  with  514.  5  nm  laser 
line. 
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Fig.3.  Size  dependence  Raman  spectra 
from  Si  dots  excited  with  514,  5 
nm  laser  line. 
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decreasing  dot  size.  The  peaks  are  observed  at  1  69  eVd  T?  ‘°  h'8her  energy  with 
and  1.4  nm,  respectively  The  averaee TtM  aund.2  06  eV  for  *>t  sizes  of  2.2  nm 
with  an  error  of  ±  0.2  nm.  We  attribute  the  fpJt,  <^>tal" ed  from  tlle  transmission  data 

recombination  in  Si  QD  and  fluctuation  in  OD  sYzes^FoMhe  ^4  spectrum  t0  band-to-band 
structure  towards  the  high  energy  (2  2  eV)  due  tn  ?  th  !'4  ^  dot  SIZe’  we  also  see 
According  to  Delerue  et  al.  [14]  ^PS  or  tinh/hLr  na"f  excitatlon  wi>h  laser  line, 
density  approximation  (LDA)  techniques  give  reliah]118  3S  .WC  as  corrected  ab  initio  local 
gap.  It  has  been  pointed  out  that  the  ouasinaitirle' ^  predlca,lons  for  Si  nanocrystal  band 
nanocrystals  require  a  large  correctiol  to  the  ha  ande*clton  coulomb  energies  of  Si 

electron  theory.  The  optical  gap  is  defined  as  E  E  -  7  n0t5!"cluded  in  the  on« 

g  -  ^lda-  2c,  where  the  quasi 
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particle  gap  EgQP  is  the  difference  in  energies  resulting  from  the  separate  addition  of  an 
electron  and  hole  to  the  system,  ILDA  is  the  self  energy  corrections  and  Ec  is  due  to  exchange 
correlation  effects.  The  SLDA  practically  compensates  for  the  screened  electron-hole 
interaction  Ecoui  with  full  dielectric  constant.  Though  our  PL  results  are  in  qualitative 
agreement,  they  fall  below  the  calculated  values  of  optical  band  gap  of  silicon  nanociystals 
with  the  excitonic  correction  [14].  It  should  be  emphasized  that  PL  energy  is  lower  than  the 
absorption  edge  for  a  given  dot  size,  and  this  difference  is  considerably  large  for  smaller  size 
dots.  Self  trapped  exciton  at  the  interface  between  Si  and  Si02  or  a  defect  state  is  considered 
to  be  responsible  for  the  observed  Stokes  shift  of  PL  emission  energy  [14]. 

The  Raman  spectra  from  Si-Si02  film  were  excited  with  the  514.5  nm  (2.41  eV)  laser 
line  of  an  argon  ion  laser  from  the  region  that  emits  green-red  luminescence.  Figure  3  shows 
representative  spectra  from  three  different  positions  of  the  sample  characterized  by  the 
average  Si  dot  diameter.  The  observed  spectra  show  phonon  softening  and  line  shape 
broadening  with  decreasing  dot  size.  As  the  average  dot  size  decreases  from  2  2  nm  to  14 
mn,  the  Raman  peak  shifts  from  51 1  cm*1  to  502  cm'1.  The  Raman  lines  are  asymmetrically 
broadened  with  a  tail  towards  the  low  energy  side.  In  the  uncorrected  data,  the  integrated 
intensity  remains  unchanged.  The  broad  line  shape  arises  from  size  distribution  in  an 
ensemble  of  dots,  which  leads  to  a  large  fluctuation  in  wave  vectors,  particularly  for  dot 
sizes  smaller  than  2  nm. 


Fig.  4.  The  radial  dependence  of  1  =  0 
optical  phonon  mode  for  a  Si 
quantum  dot  embedded  in  a  Si02 
matrix. 


Fig.  5.  Calculated  resonance  contribution 
to  the  experimentally  observed 
Raman  line  shape  of  a  1.4  nm  Si 
dot. 


In  order  to  explain  the  observed  resonance  line  shape,  we  assume  that  the  optical 
vibrations  are  confined  to  the  spherical  dot  with  effective  wave  vector  q„  =  ^/R,  where  un  is 
n  node  of  the  spherical  Bessel  function  j,,  and  R  is  the  dot  radius.  A  quadratic  bulk  Si 
dispersion  gives  the  eigen  frequencies  o)„(R).  The  eigen-frequencies  of  optical  vibrations  of 
zero  angular  momentum  (1  =  0)  mode  as  a  function  of  R  in  a  spherical  Si  dot  embedded  in 
Si02  are  shown  in  Fig.  4.  It  should  be  pointed  out  that  a  precise  determination  of  dot  size. 
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particularly  for  R  <  1  nm,  is  necessary  for  the  estimation  of  equivalent  wave  vectors  and 
corresponding  eigen-frequencies. 

The  Raman  line  shape  can  be  written  as  I  (to,  R)  [11], 


/(to,,/?)  =  Io 

_ -  1  _ ___ 

•  [ -  EnMMoh  -  EM- rvTvj’+frvOtor  EJR))+  T»,(h ©,  -  E^W). f 
_ TM _ 

^CO /  -  ft©  ,  “  /ito  n(R)f  +  Y'2 


(1) 

where  Io  is  the  product  of  scattering  efficiency  and  eigen-values  of  the  appropriate  scattering 
matrix  elements  over  all  the  intermediate  states,  R  is  the  mean  QD  radius,  EM,(R)  and  Em2(R) 
are  the  eigen  energy  for  R  in  the  \i\  and  p2  intermediate  states,  respectively,  TM  ( T)  is  the 
intermediate  state  (phonon)  life  time  broadening,  coi  and  cos  are  the  excited  laser  energy  and 
scattered  photon  energy,  respectively.  In  the  incoming  resonance  process,  coj  coincides  with 
one  of  the  electronic  states  Emj  (R)  while  outgoing  resonance  occurs  when  o)|  approaches  Ep2- 
The  resonance  conditions  are  described  as  /i(0|=  Emi(R)  for  incoming  resonance,  and  h(0\  - 
/icon(R)  =  E^(R)  for  outgoing  resonance.  We  assume  that  the  intermediate  electronic  state  is 
an  elecfron  hole  pair  state  and  the  electron  phonon  coupling  is  via  deformation  potential 
interaction.  The  energy  of  the  QD  band  edge  as  a  function  of  dot  size  is  estimated  from  the 
two-particle  exciton  calculation  [4].  The  resonant  dots  at  the  laser  photons  at  2.412  eV  have 
diameter  1.3- 1.4  nm.  Figure  5  shows  the  calculated  line  shape  for  1.4  nm  dot.  It  is  evident 
that  the  main  contribution  to  the  line  shape  comes  from  the  resonant  radii,  though  non¬ 
resonant  dots  give  broad  constant  background. 

The  Raman  line  shape  in  QD  is  indeed  a  convolution  of  dot  size  distribution,  size 
selected  vibrational  modes  in  resonance  with  incoming  photon  or  scattered  photon.  The 
resonant  radii  correspond  to  the  resonant  dots  and  give  the  main  contribution  to  the  line 
shape.  Phonons  in  different  quantum  confinement  dots  are  also  expected  to  play  a  significant 
role  in  the  scattering  process,  especially  when  the  resonance  condition  favors  smaller  radii 
The  broadening  and  asymmetry  of  the  Raman  spectrum  follows  from  the  contribution  of 
several  phonon  modes  of  dots  with  different  radii  in  the  resonance  region.  Appreciable 
contribution  comes  from  |  Re  -  R  |  <  SR,  where,  R*  =  radius  for  which  either  incoming  or 
outgoing  resonance  is  reached.  In  our  case  8R  =  0.04  nm.  The  behavior  of  Ep(R)  as  a 
function  of  R  is  taken  from  the  calculations  of  Ref.  4.  A  lifetime  broadening  of  25  meV  for 
all  intermediate  electronic  states  involved  and  phonon  half  width  of  9  cm'1  was  used  in  the 
calculation. 

CONCLUSIONS 

We  have  synthesized  Si  nanocrystals  in  a  Si02  matrix  on  quartz  plate  by  rf  co-sputtering 
The  Si  dots  are  dispersed  in  size  along  the  length  of  the  plate.  Optical  transmission 
measurements  were  used  to  estimate  the  dot  size  distribution  which  shows  peaks  around  1.0 
±  0.2  nm  and  1.4  ±  0.2  nm  with  an  asymmetrical  tail  towards  the  large  particle  size.  The  PL 
emission  from  the  dots  exhibits  broadband  in  the  red  and  green  region,  and  their  peak 
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p  onon  softening  and  large  asymmetrical  broadening  With  decreasing  dot  cj7p 
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Abstract 

We  present  combined  measurements  of  electroluminescence  (EL)  and  photoluminescence  (PL)  in  p-type  porous  silicon.  The  EL 
spectra  were  measured  using  an  electrolyte  contact  for  electron  injection  into  the  porous  face  of  the  sample.  Upon  applying  the 
current,  the  EL  intensity  first  rises  with  time,  reaches  a  maximum,  and  then  decays  to  zero.  (The  whole  process  takes  about  half 
an  hour.)  At  the  same  time,  the  peak  of  the  EL  spectrum  shifts  from  ss  850  nm  in  the  beginning  to  »  600  nm  at  the  end  of  the 
process.  The  PL,  which  was  measured  simultaneously,  peaked  at  «  750  nm  in  the  beginning  and  was  much  wider  than  all  of  the 
EL  spectra.  Towards  the  end  of  the  EL  process,  the  red  part  of  the  PL  spectrum  practically  disappears.  This  shifts  the  PL  peak 
towards  the  blue,  to  about  the  same  wavelength  as  the  EL  peak  (a 600  nm)  and  the  spectrum  becomes  much  narrower, 
comparable  to  the  EL  spectrum.  The  voltage  across  the  sample  during  the  EL  process  shows  a  moderate  increase  up  to  the  point 
where  the  EL  disappears,  and  then  the  voltage  rises  steeply.  This  behavior  is  associated  with  the  build-up  of  a  thin  oxide  layer 
on  the  porous  surface.  The  combined  results  of  EL  and  PL,  and  especially  the  disappearance  of  the  red  part  in  the 
photoluminescence  spectrum  at  the  end  of  the  EL  process,  suggest  that  in  addition  to  quantum  confinement,  localized  surface 
states  play  an  important  role  in  the  luminescence  process,  at  least  in  the  red  part  of  the  spectrum.  Such  states  may  be  associated 
with  adsorbed  species  and  disappear  upon  oxidation.  ©  2000  Elsevier  Science  S.A.  All  rights  reserved. 

Keywords:  Electroluminescence;  Photoluminescence;  Porous  silicon;  Localized  surface  states 


1.  Introduction 

Porous  silicon  [1-4],  obtained  by  electrochemical 
etching  procedures  applied  to  crystalline  silicon  sur¬ 
faces,  exhibits  high  luminescence  efficiencies  in  the  visi¬ 
ble  range.  It  is  generally  accepted  that  the  visible 
luminescence  originates  from  the  band-gap  enlargement 
due  to  quantum  confinement  in  the  porous  silicon 
nanocrystallites  [3,4].  At  the  same  time,  the  reasons  for 
the  high-efficiency  luminescence  are  still  somewhat  un¬ 
der  debate  [3-5].  It  was  suggested  that  it  is  the  amor¬ 
phous  or  microcrystalline  nature  of  the  porous  silicon 
that  is  responsible  for  the  phenomenon,  or  that  the 
formation  of  silicon  compounds  such  as  siloxene 
(Si603H6)  or  Si-H,  Si-0  and  Si-F  bonds  are  involved 
in  the  luminescence  [3,4]. 


*  Corresponding  author. 


In  this  paper  we  present  combined  measurements  of 
electroluminescence  (EL)  and  photoluminescence  (PL) 
in  p-type  porous  silicon.  It  is  shown  that  our  combined 
results  of  EL  and  PL  point  to  the  important  conclusion 
that  in  addition  to  quantum  confinement,  localized 
surface  states  play  an  important  role  in  the  lumines¬ 
cence  process  at  the  red  part  of  the  spectrum.  Such 
states  may  be  associated  with  adsorbed  species  and 
disappear  upon  oxidation. 


2.  Experimental 

The  starting  material  was  high-grade  p-type  (100) 
silicon  of  resistivity  in  the  range  0.5- 1.5  flcm.  Ap  + 
layer  was  formed  by  diffusing  metallic  A1  into  the  back 
faces  of  the  silicon  wafers  to  obtain  an  ohmic  contact. 
The  sample  was  attached  to  a  cylindrical  Teflon  cell  via 
a  Kalrez  O-ring.  The  sample  constituted  the  bottom  of 
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the  cell,  with  its  front  surface  facing  upwards.  Before 
anodization,  the  samples  were  etched  in  20%  HF.  In 
order  to  prepare  the  porous  surface  [4],  a  solution  of 
HF,  ethanol  and  water  (1:2:1)  was  poured  into  the  cell. 
A  platinum  electrode  was  immersed  in  the  solution  and 
a  spring  contact  was  attached  to  the  bottom  p  + 
contact.  The  anodization  of  the  Si  surface  was  carried 
out  with  a  current  density  of  15-50  mA  cm-2. 

The  EL  and  PL  were  measured  in  the  same  cell  using 
a  CCD-based  computer  board  spectrometer.  In  the 
former  case  the  cell  was  filled  with  an  electrolyte  con¬ 
sisting  of  an  aqueous  solution  of  NaCl.  The  EL  was 
obtained  under  constant  current  conditions,  the  current 
densities  ranging  from  1-3  mA  cm-2.  The  PL  was 
excited  by  a  45-mw  argon  laser  beam  of  488  nm 
wavelength. 


3.  Results 

In  order  to  determine  the  EL  characteristics  we  ap¬ 
plied  to  the  cell  a  constant  current  (the  electrolyte 
negative)  and  continuously  measured  the  EL  spectra.  In 
Fig.  1  we  show  some  of  these  spectra  (assorted  sym¬ 
bols)  obtained  at  different  times  of  the  EL  current  flow 
as  marked  on  the  curves.  We  see  that  upon  applying  the 
current  the  EL  intensity  first  rises  with  time,  reaches  a 
maximum,  and  then  decays  to  zero.  The  whole  process 
usually  takes  about  half  an  hour.  At  the  same  time,  the 
peak  of  the  EL  spectrum  shifts  from  ss  850  nm  in  the 
beginning  to  as  600  nm  towards  the  end  of  the  process. 


The  solid  curve  represents  the  total  luminescence  power 
as  a  function  of  wavelength,  obtained  by  integrating  the 
EL  emission  intensity.  It  is  seen  that  this  curve  forms 
an  envelope  of  the  EL  spectra  at  the  different  stages  of 
the  current  flow,  as  expected. 

The  width  of  the  different  EL  spectra  was  found  to 
be  about  120  nm  and  is  essentially  the  same  for  all 
samples  studied.  This  is  illustrated  also  by  the  data  in 
Fig.  2.  Here  we  plot  the  EL  peak  amplitude  as  a 
function  of  the  area  of  the  spectrum.  The  various 
symbols  (nine  altogether)  correspond  to  samples  pre¬ 
pared  under  different  anodization  conditions,  while  the 
recurrence  of  the  same  symbol  represents  measurements 
at  different  times  of  the  EL  current  flow.  We  note  that 
all  data  lie  on  a  straight  line,  showing  the  proportional¬ 
ity  between  the  emitted  EL  power  (area)  at  any  time 
and  the  peak  amplitude.  This  indicates  that  both  the 
shape  and  the  width  of  the  various  spectra  are  constant. 
Accordingly  we  used  the  peak  amplitude  to  represent 
the  emitted  EL  power. 

The  main  conclusions  from  Fig.  1  are  summarized  in 
Fig.  3.  Here  we  plot  the  emitted  EL  power  and  the 
wavelength  of  the  spectral  peak  as  functions  of  time. 
Also  plotted  is  the  time  variation  of  the  voltage  across 
the  sample.  Note  that  during  the  time  that  the  EL 
persists,  the  voltage  rises  moderately,  but  shortly  after 
the  disappearance  of  the  EL,  it  rises  steeply.  This 
behavior  is  associated  with  the  buildup  during  the  EL 
current  flow  of  an  insulating  layer,  most  probably  an 
oxide,  on  the  porous  surface.  In  fact  we  have  measured 
the  buildup  process  of  such  an  oxide  layer  on  crys- 


Fig.  1 .  Several  EL  spectra  (assorted  symbols)  obtained  at  different  times  of  the  EL  current  flow  (current  density  J  =  3.3  mA  cm  “  2)  as  marked 
on  the  curves.  The  solid  curve  represents  the  total  luminescence  power. 
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Fig.  2.  EL  peak  amplitude  as  a  function  of  the  area  of  the  spectrum. 
The  various  symbols  (nine  altogether)  correspond  to  samples  pre¬ 
pared  under  different  anodization  conditions,  while  the  recurrence  of 
the  same  symbol  represents  measurements  at  different  times  of  the  EL 
current  flow. 

talline  silicon  during  current  flow  under  similar  condi¬ 
tions  to  those  used  for  EL  in  porous  silicon.  Use  was 
made  of  the  space  charge  capacitance  technique  of  the 
silicon -electrolyte  interface  [5].  Fig.  4  shows  the  growth 
of  the  oxide  thickness  with  the  time  of  the  applied 
current.  We  see  that  the  current  flow  causes  an  appre¬ 
ciable  enhancement  in  the  oxide  thickness.  It  is  only  to 


Fig.  4.  Growth  of  the  oxide-layer  thickness  on  crystalline  silicon  with 
time  of  applied  current. 


be  expected  that  the  EL  current  induces  an  oxide 
growth  in  porous  silicon  as  well. 

Fig.  5  displays  the  PL  spectra  (solid  curves)  at  the 
onset  of  the  EL  (curve  1)  and  after  the  EL  disappears 
(curve  2),  We  see  that  due  to  the  EL  current  flow,  the 
red  part  of  the  PL  practically  disappears.  As  a  result, 
the  PL  spectrum  is  appreciably  narrower  and  blue 
shifted.  For  comparison  purposes  we  replotted  from 
Fig.  1  (dotted  curves)  the  total  EL  emitted  power  (curve 
3)  and  one  of  the  last  measured  EL  spectra  (curve  4). 
The  two  curves  were  normalized  to  the  respective  PL 
curves.  It  is  seen  that  each  pair  of  curves  coincide  quite 
closely. 


Fig.  3.  Emitted  EL  power  and  wavelength  of  the  spectral  peak  as  functions  of  time.  Also  plotted  is  the  time  variation  of  the  voltage  across  the 
sample. 
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Fig.  5.  PL  spectra  (solid  curves)  at  the  onset  of  the  EL  (curve  1)  and  after  the  EL  disappears  (curve  2).  For  comparison  purposes,  we  replotted 
from  Fig.  1  (dotted  curves)  the  total  EL  emitted  power  (curve  3)  and  one  of  the  last  measured  EL  spectra  (curve  4). 


4.  Discussion 

The  blue  shift  of  the  EL  peak  with  the  EL  time  and 
the  accompanying  moderate  increase  of  the  voltage 
across  the  sample  are  somewhat  similar  to  results  ob¬ 
tained  by  Bsiesy  and  Vial  [6]  and  are  probably  due  to 
the  same  mechanism.  Because  of  quantum  confinement, 
the  energy  gap  associated  with  the  nanocrystallites  is 
higher,  the  smaller  the  size  of  the  crystallite.  At  the 
beginning,  the  contribution  to  the  EL  is  mainly  from 
the  larger  nanocrystallites.  This  is  because  such  parti¬ 
cles  are  associated  with  smaller  energy  gaps  so  that 
electron  injection  from  the  electrolyte  requires  less  en¬ 
ergy  (lower  voltage).  With  time,  these  nanocrystallites 
become  covered  with  an  oxide  due  to  the  EL  current 
flow.  At  this  point  the  current  is  diverted  to  smaller 
nanocrystallites  and  electron  injection  into  such  crystal¬ 
lites  requires  a  higher  voltage.  As  pointed  out  above, 
the  measurements  were  performed  under  constant  cur¬ 
rent  conditions  and  we  see  indeed  from  Fig.  3  that  the 
voltage  across  the  cell  rises  with  the  time  of  the  EL 
current  flow.  The  EL  disappears  when  all  of  the 
nanocrystallites,  large  and  small,  become  covered  with 
an  oxide.  One  might  conceivably  assume  that  the  EL 
originates  from  band-to-band  transitions  and  thus  ac¬ 
count  for  the  blue  shift  as  the  EL  current  switches  from 
larger  to  smaller  nanocrystallites.  Such  a  model,  how¬ 
ever,  is  inconsistent  with  some  PL  data.  In  the  first 
place  we  reported  previously  [7]  that  the  excitation  of 
the  PL  is  practically  nonexistent  in  the  range  of  wave¬ 
lengths  encompassed  by  the  PL  spectrum  (above  600 
nm).  This  suggests  that  the  PL  involves  some  intermedi¬ 


ate  localized  states,  and  that  the  luminescence  results 
from  transitions  between  such  states  and  either  energy 
band.  As  the  gap  of  the  nanocrystallites  involved  in¬ 
creases,  so  does  the  energy  separation  between  the 
states  and  the  relevant  band.  Further  information 
about  the  nature  of  the  localized  states  is  provided  by 
the  observation  that  the  red  part  of  the  PL  spectrum 
disappears  at  the  end  of  the  EL  process  (Fig.  3).  This 
suggests  that  at  least  the  localized  states  responsible  for 
the  red  part  of  the  spectrum  reside  at  the  surface  of  the 
nanocrystallites.  Such  states  may  be  associated  with 
adsorbed  species  and  disappear  upon  oxidation. 

In  conclusion,  it  appears  that  in  addition  to  quantum 
confinement,  localized  states  play  an  important  role  in 
the  luminescence  process  in  porous  silicon.  Our  com¬ 
bined  results  of  EL  and  PL  suggest  that  the  states 
responsible  for  the  red  part  of  the  spectrum  reside  at 
the  surface. 
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Abstract 

Er3  + -doped  nanocrystalline-Si/Si02  composite  films  were  synthesized  by  RF  co-sputtering  of  bulk-Si,  Si02,  and  Er203  targets. 
The  visible  and  1.54  pm  emission  bands  of  the  samples  were  measured,  as  well  as  their  optical  transmission.  We  varied  the  relative 
concentrations  of  Si,  Er,  and  the  annealing  temperature.  The  nanoparticles  size  distributions  of  the  samples  were  obtained  from 
their  optical  transmission  spectra.  We  analyzed  the  dependence  of  the  characteristic  1.54  pm  Er-emission  intensity  on  the  size  and 
concentration  of  the  nanocrystalline  particles  and  the  dependence  of  the  emission  on  the  preparation  conditions.  ©  2000  Elsevier 
Science  S.A.  All  rights  reserved. 
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nanocrystalline  silicon  due  to  carrier  confinement  ap¬ 
pears  to  favor  energy  transfer  to  Er3  + .  Indeed,  in  Si02 
films  containing  nanocrystalline  silicon  and  Er,  a  corre¬ 
lation  between  visible  emission  from  the  nanocrystal¬ 
lites  and  IR  emission  from  Er3  +  was  observed  [5], 
which  further  corroborate  that  the  IR  emission  of  Er3  + 
is  mediated  by  the  nanocrystallites. 

In  this  paper  we  investigate  the  dependence  of  the 
1.54  pm  emission  intensity  on  the  nanocrystal  size  and 
concentration  in  Er3+  doped  nanocrystalline  Si/Si02 
composite  films.  The  films  were  prepared  by  varying  the 
relative  concentration  of  Si,  Er,  and  the  annealing 
temperature. 


2.  Results  and  discussion 

The  material  was  synthesized  by  RF  co-sputtering. 
The  sputtering  setup  that  was  used  to  prepare  the 
samples  included  a  Varian  system  in  a  target  downside 
—  substrate  upside  vertical  configuration.  The  chamber 
had  62-1  volume  connected  to  a  2500  Is"1  diffusion 
pump.  A  Si02  disk  with  13.2  cm  in  diameter  was  used 
in  the  target  configuration.  The  substrates  were  quartz 
plates  with  dimensions  12.6  x  1.1  x  0.09  cm3.  A  piece  of 
crystalline  silicon  wafer  with  7  cm2  of  area  was 
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j  1.  Introduction 

I  Recently,  efficient  light  emission  has  been  demon¬ 
strated  by  erbium  (Er)  incorporation  into  silicon  [1]. 

|  The  erbium  ion,  in  its  3+  state,  emits  photon  in  the 
infrared  (IR)  region  at  1.54  pm  resulting  from  an  intra 
|  4f  transition  between  the  4I13/2  state  to  the  4I15/2  state. 

;  This  emission  is  particularly  suitable  for  the  low  loss 
i  transmission  through  silica  based  optical  fiber.  The 
j  strong  quenching  of  luminescence  at  room  temperature 
(RT),  through  nonradiative  decay  channels,  is  consid- 
!  ered  to  be  a  principal  limitation  for  the  development  of 
j  Er-doped  silicon  devices.  It  has  been  shown  that  in  the 
presence  of  impurities  such  as  O  and  F,  the  reduction  in 
j  the  quenching  as  well  as  higher  erbium  solubility  can  be 
achieved  [2,3]. 

Er-doped  porous  silicon,  on  the  other  hand,  shows 
strong  1.54  pm  emission  even  at  RT  and  a  considerable 
weak  thermal  quenching  [4].  It  is  argued  that  in  the 
s;  !  porous  silicon,  the  Er3+  excitation  results  from  the 

I  recombination  of  the  photoexcited  carriers  in  the  silicon 
nanocrystallites  and  subsequent  energy  transfer  to  Er3  + . 

|  The  modification  in  the  band  structure  properties  of 
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mounted  on  top  of  the  Si02  target  disk  just  below  one 
end  of  the  substrate  deposition  area.  With  such  a 
configuration  a  density  gradient  of  silicon  along  the 
sputtered  film  can  be  obtained  in  one  sputtering  session. 
The  Er203  material  was  obtained  as  a  powder  and  then 
pressed  and  prepared  as  pellets. 

The  deposited  films  had  an  area  of  1 1  x  0.44  cm2  and 
the  thickness  was  in  the  order  of  ~  3.5  pm  after  6  h  of 
deposition.  The  sputtering  was  made  with  a  residual  gas 
pressure  ~5  x  10“ 7  Torr,  Argon  (99.9995%)  pressure 
of  19  mTorr,  1100  V  substrate-to-target  voltage,  180  W 
forward  power,  and  2.5  W  reflected  power.  The  sub- 
strate-to-target  distance  was  5  cm.  The  substrate  tem¬ 
perature  raised  up  to  -  130°C  during  the  sputtering 
period.  Previously  to  the  deposition,  the  quartz  sub¬ 
strates  were  cleaned  by  immersion  in  a  HF  (20%)  water 
solution  for  30  min  and  then  washed  with  deionized 
water  and  dried  with  a  nitrogen  flux.  After  sputtering, 
the  samples  were  annealed  in  a  N2-rich  atmosphere 
using  a  Lindberg  Heavy  Duty  oven.  The  samples  were 
kept  in  the  oven  during  the  heating-up  and  cooling- 
down  periods. 

Er-doped  samples  were  synthesized  by  adding  to  the 
sputtering  target  configuration  a  line  of  small  ( ~  4  x  9 
mm2)  Er203  pellets  oriented  along  the  substrate.  The 
distance  between  each  pellet  and  the  next  was  —  7  mm 
and  they  were  located  just  on  top  of  the  Si  and  Si02 
targets.  With  this  target  configuration  a  quasi-uniform 
Erbium  concentration  was  sputtered  along  the  film.  The 
samples  were  annealed  at  different  temperatures  from 
1200  to  700°C.  Auger  Electron  Spectroscopy  (AES) 
gave  an  erbium  concentration  of  3.9%  along  the  film.  In 


Wavelength  (nm) 


Fig.  1.  Infrared  emission  spectra  from  uniformly  Er-doped  Si/Si02 
films  annealed  at  700°C  and  excited  with  488  nm  laser  line.  The 
number  on  each  spectrum  designates  the  segment  number  of  the  film. 


order  to  correlate  the  optical  properties  of  the  material  ' 
with  different  density  of  the  constituents,  each  film  was 
divided  in  50  segments  (0.22  cm.  in  length)  and  labeled 
successively  with  the  Si02-rich  end  as  the  No.  1  and  the 
other  end  as  No.  50.  We  measured  the  optical  transmis¬ 
sion  and  the  photoluminescence  properties  of  each  seg¬ 
ment  of  the  films  in  both  the  visible  and  the  IR  (near 
1.54  pm)  regions.  The  transmission  spectra  were  ob¬ 
tained  using  a  HP-8452A  spectrophotometer  and  the 
photoluminescence  was  measured  with  an  ISA-TRIAX  * 
320  spectrometer. 

Fig.  1  shows  the  IR-emission  between  1500  and  1600 
nm  for  a  sample  that  was  annealed  at  700°C.  Each  k 
curve  corresponds  to  a  different  position  along  the  film 
as  written  in  the  graph.  Above  position  30,  Fig.  1  shows 
the  characteristic  emission  of  Er  +  3.  The  intensity  of  the 
IR  emission  increases  as  the  silicon  concentration  does. 

This  gives  a  clear  indication  that  the  Er-emission  is 
enhanced  by  the  presence  of  the  silicon  nanoparticles. 
Positions  30  and  below  that  correspond  to  lower  con¬ 
centration  of  silicon  give  almost  no  emission.  As  men¬ 
tioned  above,  it  is  expected  that  both  the  density  and 
the  size  of  the  nanoparticles  reduce  when  approaching 
the  lower  silicon  concentration  side  of  the  film.  At  the 
same  time  during  the  sputtering,  oxygen  released  from 
Er203  can  react  with  silicon  and  reduce  even  further  the 
available  amount  of  silicon  required  for  the  formation 
of  the  nanoparticles  in  this  part  of  the  film.  The  PL 
measurements  of  this  film  in  the  visible  region  give  no 
direct  emission  from  the  nanoparticles  at  any  of  the 
studied  positions.  This  can  be  explained  as  the  quench¬ 
ing  of  the  visible  PL  because  of  the  energy  transfer  to 
the  Er-emission  for  the  samples  where  the  IR  emission 
is  observed. 

In  order  to  correlate  the  IR  emission  with  the  size  of 
the  nanostructures,  we  determined  the  size  of  the  parti¬ 
cles  at  positions  35  and  45  by  optical  means.  Fig.  2(a) 
shows  the  optical  transmission  spectra  in  the  visible 
region  of  the  film  at  those  positions.  The  optical  ab¬ 
sorption  edge  shifts  are  attributed  to  the  quantum 
enlargement  of  the  effective  energy  gap  of  the  particle 
when  the  size  is  reduced.  From  the  optical  transmission, 
we  calculated  the  particle  size  distribution  by  assuming 
that,  due  to  the  electron  quantum  confinement,  the 
energy  dependence  of  the  absorption  coefficient  for  a 
silicon  nanoparticle  shows  a  sudden  increase  from  zero 
to  a  finite  value  at  the  energy  corresponding  to  the 
enlarged  energy  gap  [6].  Under  this  assumption  the  size 
distribution  is  proportional  to  the  derivative  of  the 
transmission  versus  energy  curve.  Fig.  2(b)  shows  the 
normalized  size  distributions  obtained  from  the  optical 
transmission  data  at  the  same  points,  35  and  45.  The 
maximum  of  the  distribution  changes  from  1.18  nm  at 
point  35  to  1.8  nm  at  point  45. 

The  as-sputtered  silicon  concentration  at  each  site 
was  obtained  from  the  determination  of  the  film  thick- 
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Fig.  2.  (a)  Optical  transmission  spectra  from  segment  35  and  45  of  a 
uniformly  Er-doped  Si/Si02  film  annealed  at  700°C.  (b)  Calculated 
normalized  Si-nanocrystal  size  distribution  at  the  positions  where 
transmission  spectra  were  taken. 

ness,  measured  with  a  Sloan  Dektak  Surface  Profilome- 
ter,  and  the  Hanak  et  al.  algorithm  [7].  With  the  silicon 
concentration  data  and  the  particle  average  size  as¬ 
signed  to  each  site  from  Fig.  2(b),  the  density  of 
particles  can  be  estimated  at  each  position  assuming 
that  all  the  silicon  excess  has  been  grown  into  particles. 
In  that  case  the  silicon  density  at  position  35  is  0.93% 
and  the  particle  density  ~  1.08  x  1019  part  per  cc  and 
at  position  45  the  silicon  density  is  2.3%  and  the 
particle  density  —  7.53  x  1018  part  per  cc.  These  num¬ 
bers  and  the  fact  that  the  maximum  of  the  IR-emission 
at  position  45  is  larger  than  the  IR-emission  at  position 
35,  brings  us  to  the  conclusion  that  the  energy  transfer 
in  the  sample  with  particles  of  the  size  of  ~  1.8  nm  is 
enhanced  as  compared  with  the  sample  with  particles  of 
~  1.18  nm  in  size.  It  is  worth  noting  that,  according  to 
the  statement  that  oxygen  coming  from  the  sputtering 
of  the  Er203  pellets  can  oxidize  the  silicon  during  the 
film  formation,  the  silicon  concentration  used  on  the 
above  calculations  can  be  overestimated  specially  at 
position  35,  where  the  silicon  concentration  is  lower. 

To  study  how  the  variation  of  the  erbium  concentra¬ 
tion  changes  the  optical  properties  of  the  system,  the 
sputtering-target  configuration  was  changed  by  replac¬ 


ing  the  line  of  small  pellets  oriented  along  the  film  by 
only  one  Er203  pellet  with  13  mm  in  diameter  located 
below  the  center  of  the  substrate.  In  such  a  way  the 
erbium  concentration  will  change  down  from  close  to 
the  center  to  both  ends.  Fig.  3(a)  shows  the  IR  emission 
at  positions  20  and  28  of  a  film  that  was  prepared  with 
the  new  target  configuration  and  annealed  at  900°C. 
The  figure  shows  a  reduction  of  the  IR  emission  at 
position  20,  where  the  amount  of  erbium  is  lower  as 
compared  with  position  28.  Fig.  3(b)  shows  the  corre¬ 
spondent  visible  emission  that  is  related  with  the  PL  of 
the  silicon  nanoparticles.  The  intensity  of  the  visible  PL 
increases  when  the  IR  emission  intensity  reduces. 

Another  sample  was  prepared  using  the  new  target 
configuration  and  annealing  the  film  at  700°C.  This 
sample  showed  strong  IR  emission  and  non-detectable 
visible  emission  as  in  the  case  of  the  first  sample  that 
was  also  annealed  at  the  same  temperature.  Fig.  4(a) 
shows  the  dependence  of  the  IR  emission  intensity  as  a 
function  of  the  position.  The  curve  shows  that  the  PL 
decreases  when  the  position  moves  to  the  ends  where 
the  erbium  concentration  reduces  and  the  PL  curve  is 
not  symmetric  with  larger  intensity  where  the  silicon 
concentration  and  silicon  nanoparticles  are  larger  as  in 

Wavelength  (nm) 
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Fig.  3.  (a)  Infrared  emission  spectra  from  non-uniformly  Er-doped 
Si/Si02  films  annealed  at  900°C  excited  with  488  nm  laser  line.  The 
number  on  each  spectrum  designates  the  segment  number  of  the  film, 
(b)  Visible  PL  at  the  same  positions,  where  IR  emission  is  also  seen. 
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Position 


Fig.  4.  (a)  Variation  of  1.54  pm  peak  intensity  with  position  on  the 
film,  and  (b)  variation  of  Er  and  oxygen  concentration  measured  by 
Auger  spectroscopy  with  position  on  the  film.  Positions  1  and  50  are 
assigned  the  SiCL,  rich  and  Si  rich  sites,  respectively,  on  the  film. 

the  case  of  the  sample  of  Fig.  1.  An  interesting  feature 
is  observed  at  positions  where  the  erbium  concentration 
is  maximum  (around  the  center  of  the  film),  the  PL  is 
quenched  at  these  sites.  This  observation  also  suggests 
that  the  silicon  oxidized  at  this  region  because  of  the 
increasing  amount  of  oxygen  coming  from  the  Er203. 
This  conclusion  is  supported  by  the  fact  that  there  is  a 
bleaching  of  the  film  at  these  positions.  We  included  in 
Fig.  4(b)  the  erbium  and  oxygen  concentrations  mea¬ 
sured  by  Auger  spectroscopy.  The  curves  support  the 
fact  that  the  increasing  of  the  Er-concentration  is  also 
accompanied  by  large  amounts  of  oxygen  due  to  the 
target  composition  that  will  promote  the  silicon  oxida¬ 
tion  and  the  quenching  of  the  IR  emission. 


3.  Conclusions 

We  have  synthesized  Er3  +  doped  nanocrystalline 


Si/Si02  composite  films  on  quartz  plates  by  RF  sputter¬ 
ing  using  Er203,  crystalline  silicon  and  Si02  targets. 
During  the  process  we  varied  the  relative  concentration 
of  Si,  Er  and  the  annealing  temperature.  The  character¬ 
istic  1.54  pm  emission  is  found  to  be  sensitive  to  the  Si 
nanocrystallite  density,  size,  Er  concentration  and  oxy¬ 
gen  incorporation  in  the  film.  At  low  temperature  an¬ 
nealing  (700°C),  both  uniform  and  non-uniformly 
Er-doped  samples  exhibit  pronounced  IR  emission  in 
the  high  Si  nanocrystallite  density  region.  This  clearly 
suggests  that  the  IR  emission  is  being  enhanced  by  the 
presence  of  the  Si  nanocrystallites  and  that  high  Si 
nanocrystallite  density  is  essential.  At  higher  annealing 
temperature  (900°C),  both  visible  and  IR  emission  were 
measured  and  we  observed  a  correlation  between  the 
1.54  pm  emission  and  the  energy  transfer  between  the 
Er3+  and  the  Si  nanoparticle.  On  the  other  hand,  IR 
emission  efficiency  depends  not  only  on  the  silicon 
particles  density  and  Er  concentration  but  also  on  the 
size  of  the  nanocrystallites.  We  observed  that,  at  the 
particles  size  range  we  were  working  (~  1.2  to  ~2 
nm),  larger  particles  produce  larger  enhancement  of  the 
IR  emission.  The  use  of  Er203  as  sputtering  target  can 
also  introduce  large  amounts  of  oxygen  and  prevent  the 
silicon  particles  formation  when  increasing  density  of 
erbium  is  required. 
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ABSTRACT 

We  prepared  Eu-doped  films  of  Si  nanoparticles  embedded  in  Si02  using  pellets  of 
Eu203  by  sputtering.  We  studied  their  photoemission,  transmission  and  fluorescence  to  obtain 
data  about  their  composition  and  particle  size  and  the  Eu  interaction  characteristics.  We  were 
able  to  incorporate  Eu(III)  into  the  Si  nanoparticle  /  Si02  host.  We  also  found  we  obtained  Eu(II) 
in  the  process.  We  found  a  lowering  of  photoluminescence  intensity  with  lowering  of 
temperature.  An  as  yet  unanswered  question  is  the  reason  for  the  intense  whitish  luminescence 
found  in  some  regions  of  the  samples.  Some  involvement  with  Eu(II)  is  suspected.  Eu(IQ)  related 
peaks  were  only  observed  where  the  size  distribution  peak  of  the  nanoparticles  was  lower  than 
1.3nm.  Whitish  luminescence  was  related  to  peak  sizes  ranging  from  1.1  nm  to  1.4nm.  Annealing 
the  samples  had  clear  effects  upon  their  photoluminescence,  but  did  not  necessarily  involve 
changes  in  particle  sizes,  nor  were  these  size  changes  necessary  to  increase  luminescence.  The 
Eu  doping  has  a  tendency  to  halt  the  annealing  effects  on  size  and,  when  changes  did  occur,  the 
particles  generally  became  smaller. 

INTRODUCTION 

Europium  is  an  interesting  material  when  incorporated  as  a  dopant  into  solid  state 
materials.  This  is  because  trivalent  Eu(III)  ions  can  be  used  as  a  structure  probe,  and  due  to  its 
large  quadrupole  splitting  and  spin-lattice  relaxation  times,  Eu(III)  is  particularly  attractive  in 
high-resolution  laser  spectroscopy  [1],  One  proposed  application  is  in  a  frequency  domain 
optical  data  memory,  based  on  the  large  ratio  of  inhomogeneous-to-homogeneous  line  widths 
present.  In  some  solid  state  materials  the  Eu(III)  dopant  ions  exhibit  the  longest  ever  dephasing 
time  of  822 ps  and  the  narrowest  homogeneous  line  width  at  low  temperature  [2],  This  has 
stimulated  interest  in  searching  for  appropriate  materials  to  host  Eu(IH)  ions,  which  could  be 
optimized  for  operation  at  room  temperature. 

We  have  produced  and  studied  Si/Si02  material  in  which  visible  light  emitting  Si 
nanoparticles  are  embedded  in  a  Si02  matrix.  We  made  subsequent  studies  incorporating  another 
rare-earth  metal,  Er,  into  the  Si  nanoparticle  embedded  Si02  matrix  and  verified  that  the 
luminescence  of  Er  is  enhanced  due  to  energy  transfer  between  the  nanoparticles  and  the  Er  [3], 

These  considerations  led  us  to  attempt  to  ascertain  the  interactions  between  Eu  and  this 
new  host  material.  Due  to  the  amorphous  nature  of  the  Si02  material  prepared  by  this  particular 
technique,  the  dopant  europium  ions  must  experience  large  variations  of  their  environment. 
Additionally,  some  Eu(III)  ions  may  be  placed  at  the  Si  nanoparticle  /  Si02  interface  while  the 
others  are  inside  the  Si02  host.  This  would  produce  fluorescence  line  broadening.  Also,  our  use 
of  Eu2C>3  as  the  source  of  our  dopant  could  lead  to  interactions  between  the  oxygen  thus 
incorporated  into  our  samples  and  the  Si,  as  had  apparently  happened  in  our  previous  Er 
incorporation  experiments. 
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EXPERIMENT 

We  combined  E112O3,  Si  and  Si02  producing  a  series  of  samples  by  Ar  sputtering 
deposition  on  quartz  substrates.  We  produced  samples  with  Si  and  Si02;  Eu203  and  Si02;  Eu2C>3 
and  Si;  and  Eu2C>3,  Si  and  Si02  .  The  samples  were  made  by  a  12hr  deposition.  The  plates 
holding  the  sample  and  the  target  materials  reached  temperatures  in  the  vicinity  of  110°C.  The 
target  voltage  was  1 100V  with  a  0. 12A  plate  current. 

The  largest  portion  of  the  target  was  a  13.2cm  diameter  disk  of  Si02  (except  for  the 
Eu203  with  Si  sample,  in  which  case  it  was  a  disk  of  Si).  When  including  Si  in  the  sample  (for 
the  nanoparticles)  we  placed  a  3.00cm  by  2.30cm  Si  rectangle  on  the  border  of  the  Si02  disk. 
The  14cm  x  1cm  x  0.1cm  rectangular  quartz  substrate  was  placed  with  one  end  above  the  Si 
rectangle  and  oriented  almost  along  the  diameter  of  the  Si02  disk.  The  target  to  substrate 
distance  was  about  5cm.  This  produced  a  gradient  in  the  Si  concentration  along  the  sample.  The 
Eu203  pellet  was  a  half-circle  placed  with  its  1cm  diameter  perpendicular  to  the  length  of  the 
substrate.  For  different  samples,  the  pellet’s  position  was  varied  from  being  next  to  the  Si 
rectangle,  to  beneath  approximately  halfway  along  the  substrate,  to  a  bit  beyond  the  substrate  on 
the  side  opposite  to  the  Si  rectangle.  This  way  we  could  see  the  effects  of  a  variety  of  Si-Eu 
combination  proportions.  The  Eu2C>3  /  Si02  sample  had  the  pellet  near  the  center.  The  Eu2C>3  /  Si 
sample  had  it  near  one  edge  of  the  sample.  The  samples  were  marked  lengthwise.  Position  50  is 
closest  to  the  Si . 

We  measured  the  luminescence  and  transmission  of  the  samples,  and  then  annealed  them 
at  700°C,  except  in  the  case  of  pure  Si  with  Si02,  which  was  annealed  at  1100°C.  Annealing 
took  place  in  a  Lindberg  Heavy  Duty  Oven  in  a  N2  or  Ar  rich  atmosphere.  Annealing  the  as- 
sputtered  Si/Si02  at  1100°C  results  in  Si  nanoparticles  with  sizes  ranging  from  about  1  to  greater 
than  3nm,  with  the  dominant  sizes  ranging  from  about  1.0  to  2. 1  nm  as  we  go  from  position  1  to 
position  50.  We  observed  that  higher  annealing  temperatures  result  in  smaller  particles  and  less 
variation  in  size  at  each  position  and  overall.  Annealing  should  help  crystallization  and  other 
bond  formation  processes,  and  perhaps  nanoparticle  growth,  but  some  oxidation  could  also  take 
place  with  02  from  the  atmosphere  (or,  in  the  Eu-doped  samples,  from  the  Eu2(>3).  For  Er  doped 
Si/Si02  700°C  was  sufficient  [4],  and  we  tried  it  for  Eu  . 

The  photoluminescence  was  measured  with  an  Argon  laser  514.5nm  beam.  We  used  an 
ISA  Triax  320  spectrometer.  We  measured  the  luminescence  of  the  Eu2C>3  pellet  and  a  EuCfj 
pellet  to  characterize  the  Eu(III)  luminescence  peaks  for  comparison  with  our  samples.  The 
results  were  comparable  to  others  found  in  the  literature  [4].  A  CW  Ar  laser  at  488nm  and  a  high 
power  pulsed  YAG  laser  were  used  for  further  studies  to  again  observe  the  general  behavior  and 
also  identify  the  components  of  the  luminescence  via  luminescence  decay  and  time  resolved 
spectroscopy  measurements. 

RESULTS 

The  Eu203  /  Si02  samples  were  transparent  and  colorless.  The  Eu2(>3  /  Si  sample  looked 
like  Si.  The  Si  /  Si02  sample  was  caramel  like  at  the  Si  rich  end,  getting  lighter  in  color  until 
about  three  centimeters  from  the  Si  side  it  was  almost  transparent,  and  finally  completely 
transparent  at  the  Si  poor  end.  The  Eu2C>3  /  Si  /  Si02  samples  were  the  same  as  the  Si  /  Si02 
samples  except  that  in  the  vicinity  of  where  the  Eu  pellet  had  been  the  Si  caramel  color  (when 
present)  gave  way  to  transparency.  If  this  “bleaching  effect  was  due  to  Eu2C>3  donating  some  of 
its  oxygen  atoms  to  Si,  the  oxidizing  of  Si  would  have  to  be  accompanied  by  a  change  in  the  Eu 
bonding.  In  Eu2C>3  we  have  Eu(III)  present.  Upon  donating  the  oxygen  it  could  become  Eu(II), 
which  is  also  a  relatively  stable  valence  state  for  Eu. 
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The  photoluminescence  (PL)  studies  with  the  514.5nm  Ar  laser  showed  two  types  of 
spectra.  On  the  one  hand  we  saw  wide  bands  (from  550nm-800nm)  with  peaks  in  the  vicinity  of 
600-680nm.  On  the  other  hand,  we  saw  a  well  defined  peak  at  619nm.  This  peak  only  occurred 
in  Eu  doped  samples  and  corresponded  to  the  5Do-7F2  Eu(H[)  transition.  Quite  a  few  times  we 
saw  the  Eu(IH)  peak  piggybacked  atop  one  of  the  larger  bands.  In  those  cases  the  bands  tended  to 
be  triangular  in  shape  (Fig.  1). 

The  as-sputtered  samples  were  also  viewed  under  the  light  of  a  366nm  u-v  lamp.  To  the 
naked  eye,  a  strong  whitish  luminescence  was  observed  on  either  side  of  the  EU2O3  pellet 
position  in  the  EU2O3  /  Si  /  Si02  samples,  while  the  position  corresponding  to  the  pellet  itself 
tended  to  be  transparent  and  not  luminescent.  Farther  off  this  luminescence  ceased.  The 
positions  where  the  whitish  luminescence  occurred  were  somewhat  affected  by  annealing  later.  It 
should  be  noted  that  the  whitish  luminescence  occurred  in  non-caramel  areas  in  the  sample  but 
not  in  the  most  Eu  rich  bleached  areas  centered  where  the  Eu  pellet  was  located.  The  lone 
exception  may  have  occurred  where  the  Si  rich  side  was  bleached  when  the  Eu  pellet  was  next  to 
the  Si  rectangle.  This  luminescence  had  not  been  observed  either  in  the  Si,  SiC>2  samples  or  a 
series  of  Er203,  Si,  Si02  samples.  The  source  for  this  luminescence  was  not  clear.  We  considered 
that  the  Eu  could  be  serving  as  a  nucleation  center  for  the  growth  of  Si  nanoparticles  and  that 
they  be  present  even  before  annealing,  but  we  had  never  observed  a  whitish  luminescence  nor 
one  as  strong  as  we  were  seeing  due  to  only  Si  nanoparticles.  Another  option  was  that  the  Eu 
itself  was  emitting. 

The  transmission  measurements  were  used  to  obtain  the  size  distributions  of  the  Si 
nanoparticles  through  the  derivative  of  the  transmission  curves.  From  each  of  these  distributions, 
the  peak  position  was  used  to  obtain  the  size  of  the  most  numerous  particles  [5],  We  assumed  the 
absorption  is  from  the  nanoparticles  and  that  they  are  of  spherical  shape  for  these  calculations. 
Distribution  peak  particle  sizes  were  found  to  range  from  l.Onm  to  2.1nm  among  the  EU2O3  /  Si  / 
Si02  samples  while  sizes  go  up  to  2.2nm  in  the  Si  /  SiC>2  sample.  As  we  continue  it  is  to  these 
distribution  peak  particle  sizes  to  which  we  will  refer. 

All  the  Si  containing  samples  showed  variation  in  size  of  the  particles  with  position  along 
the  sample.  However,  it  was  present  before  annealing  and,  we  did  not  detect  significant  changes 
in  that  particle  size  (Fig.  2a)  upon  annealing  except  in  the  sample  where  the  Eu  pellet  was  on  the 
side  opposite  to  the  Si  rectangle.  In  this  latter  case  the  changes  occurred  only  after  about  the 
midpoint  of  the  sample,  on  the  Si  rich  side.  In  fact,  in  this  case  the  particles  did  not  grow  but 
rather  reduced  in  size.  These  changes  may  have  been  due  to  an  oxidation  process  with  oxygen 
from  the  EU2O3  doping  material.  Interestingly  enough,  also,  the  changes  in  size  began  just  at  the 
location  where  the  whitish  luminescence  occurred.  It  would  seem  that  the  EU2O3  was  interfering 

because  the  particles  sizes  vary  along  the 
whole  length  of  the  Si  /  SiC>2  sample  (Fig. 
2b). 

Upon  annealing  the  luminescence 
peaks  in  the  Eu  opposite  side  to  Si  sample 
shifted  to  lower  wavelengths  values  after 
annealing.  This  also  occurred  in  the  Eu  and 
Si  on  the  same  side  sample.  The  Eu 
opposite  side  to  Si  sample  showed  a  large 
increase  in  absolute  and  relative 

Figure  1  Photoluminescence  with 
514.5nm  excitation.  Sample  with:  EU2O3 
(middle)  /  Si  /  SiC>2,  positions  20,  25,  30. 
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luminescence  in  the  Si  rich  Eu  poor  side,  beginning  at  the  whitish  luminescence  location  and 
correlated  with  the  increasing  particle  size  region  where  the  change  in  size  also  occurred.  No 
changes  in  size  occurred  in  the  Eu  in  the  middle  sample.  In  all  Eu  containing  samples,  other 
luminescence  changes  also  occurred  (Fig.  2a). 

The  very  bleached  Eu  rich  locations  tended  to  show  the  Eu(III)  peak  both  before  and  after 
annealing.  They  did  not  occur  at  the  high  Si  concentration  areas.  Particle  sizes  ranged  from 
1.1  nm  to  1.3nm  in  these  regions.  To  the  sides  of  these  regions  we  noticed  an  intense  band  of 
luminescence  similar  in  shape  and  peak  position  (640-650nm)  to  that  of  the  Si  /  SiC>2  bands. 
These  correspond  to  the  positions  at  which  we  saw  the  whitish  luminescence  present  before  and 
after  annealing.  The  particles  ranged  from  1.1  nm  to  1.4nm  there.  Farther  from  these  areas,  we 
reach  the  positions  at  which  we  could  find  Eu(III)  peaks  piggybacked  atop  bands.  These  were 
seen  after  annealing  in  locations  where  previously  there  was  a  band.  The  bands  luminescence 
was  reduced  relative  to  the  other  locations  after  annealing.  They  occurred  in  the  lower  Si 
concentration  areas.  Particle  sizes  here  were  lower  than  1.3nm.  Only  one  piggyback  case 
occurred  between  the  Eu(III)  peak  and  the  whitish  luminescence  position.  It  occurred  in  a 
relatively  large  Si  concentration  (for  these  peaks)  and  also  large  Eu  concentration  area  (right  next 
to  the  Eu(III)  peak)  in  the  Eu  at  the  middle  sample.  It  almost  disappeared  completely  both 
relative  to  the  rest  and  in  absolute  terms.  Particle  size  here  and  in  the  whitish  luminescence 
position  next  to  it  was  1.3nm,  and  in  the  pure  Eu(III)  peak  position  particle  size  was  1.1  nm. 

With  respect  to  the  cause  for  the  strong  whitish  luminescence,  various  sources  indicate 
that  Eu(II)  can  luminesce  in  the  visible  range,  with  a  very  wide  band  going  from  IR  to  UV  with  a 
maxima  in  the  green  region  of  the  spectrum  [6,  7,  8,  9],  Given  that  it  was  possibly  due  to  Eu(II) 
but  our  deposition  source  had  Eu(III)  we  had  to  check  for  a  clear  Eu(H)  signal.  The  problem  was 
that  in  the  whitish  luminescence  regions  of  our  samples  we  would  observe  a  Si  /  SiC>2  like  curve 
upon  illuminating  with  our  514.5nm  laser  (albeit  a  very  strong  one  relative  to  the  rest).  Finally, 
to  excite  Eu(II)  we  might  need  photons  in  the  blue  to  u-v  region  [9],  Since  we  observe  the 
whitish  luminescence  under  u-v  light,  this  might  indicate  that  Eu(II)  is  the  actual  source,  even  if 
we  only  observe  the  nanoparticles  with  our  5 14.5nm  or  other  blue  source. 

Time-resolved  fluorescence  measurement  was  employed  to  identify  the  origin  of  the 
observed  photoluminescence.  Using  a  CW  Ar  laser  tuned  at  488  nm  the  general  behavior  of  the 
observed  integrated  fluorescence  was  characterized  by  a  broad  band  feature  with  the  peak 
varying  in  the  region  of  600  -  750  nm,  recognizable  as  the  luminescence  from  Si/SiC>2  host 
material  [10-12],  The  measured  luminescence  decay  was  very  fast  with  a  lifetime  typically 


Figure  2  Sizes  of  particles  (in  nm)  and  maximum  photoluminescence  values  (in  arb.  units) 
versus  position.  Solid  line  and  crosses:  photoluminescence  and  particle  size,  respectively,  after 
annealing.  Dotted  line  and  circles:  photoluminescence  and  particle  size,  respectively,  before 
annealing 
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being  <  10  ns  [10],  For  some  samples  the  fine  PL  features  from  Eu(III)  emission  which  is 
overlapped  with  host  emission  could  be  identified. 

Using  a  high-power  pulsed  YAG  laser  with  a  pulse  width  of  30  ps  as  a  light  source  the 
PL  decay  at  ~  600  -  620  nm  was  carefully  examined.  We  identified  2  components.  Right  after 
excitation  a  fast  decay  component  occurs  with  a  lifetime  of  800  ns  followed  by  a  buildup  of  a 
long-lived  signal  which  has  a  lifetime  of  2  ms.  Time-resolved  spectroscopy  was  thereafter  used 
to  separate  the  two  PL  components. 

Setting  a  gate  width  of  400ns  with  the  delay  at  300ns  (to  eliminate  the  host  emission)  a 
broadband  luminescence  was  obtained,  it  was  PL  from  Eu(II)  ions.  In  order  to  verify  this,  a 
sample  without  the  Europium  dopant  was  measured.  With  the  same  setup,  in  this  other  sample 
no  fast  PL  component  nor  long-lived  component  could  be  observed.  This  pointed  to  the 
possible  origin  of  the  observed  fast  PL  in  Eu-doped  sample  as  divalent  europium  ions. 

Eu(H)  has  an  f  electron  configuration.  The  ground  state,  3.8  eV  below  the  conduction 
band,  is  the  8S7/2  level  and  the  first  excited  state  is  6P  level  of  the  4f65d  electron  configuration. 
The  characteristic  optical  absorption  band  of  ground  state  8S7/2  to  excited  state  6P  transition  is  at 
around  270  -  300  nm.  Using  the  3rd  harmonics  of  an  intense  YAG  laser  operating  at  355  nm,  this 
PL  was  just  observable,  but  overlapped  by  a  strong  Eu(III)  emission.  Using  the  laser  operating  at 
532  nm  with  a  high-power  output  of  -27  GW/cm2,  we  increased  the  relative  intensity  of  Eu(II) 
w/r  Eu(III).  The  two-photon-absorption  (TP A)  induced  Eu(II)  emission  is  shown  in  Fig.  3.  This 
photoluminescence  has  a  lifetime  of  ~800ns,  the  peak  is  at  -540  nm,  and  FWHM  is  -4200  cm'1. 
This  emission  profile  looks  similar  to  the  PL  from  the  Si  nanoparticles  and  from  defect  centers  in 
SiC>2  as  well  [13],  But  the  lifetime  measurement  as  well  as  the  experimental  results  with  the  non- 
doped  sample  convinced  us  in  making  the  judgement. 

The  luminescence  from  Eu(III)  ions  showed  a  well  resolved  feature  in  the  range  of  575 
nm  to  705  nm.  The  PL  intensity  versus  temperature  showed  an  interesting  behavior.  It  was 
found  to  decrease  with  lower  temperature.  Presumably  it  is  associated  with  the  phonon-assisted 
photodynamical  processes  in  the  host  conduction  band.  Its  emission  decay  rate  is  3  orders  of 
magnitude  greater  than  that  of  Eu(II).  With  a  large  delay  of  the  sampling  gate  the  fluorescence 
spectra  of  Eu(HI)  was  obtained  as  in  Fig.4.  It  shows  a  nature  typical  of  Eu(III)  in  an  amorphous 
host.  Its  PL  profile  consists  of  5  emission  groups.  In  the  spectral  range  from  584  to  609  nm 
there  is  a  triplet  structure.  At  first  glance  it  looks  like  the  triple  components  of  the  transition 
from  metastable  5Do  to  7Fi  state  with  j  =  0,  ±  1.  If  the  first  peak  at  580  nm  of  the  triplet  were 
assigned  to  5Dq  -  7Fi  transition,  the  5Do  -  7Fq  transition  would  occur  at  an  even  shorter 


500  550  600  650  700 

Wavelength  (nm) 


Figure  3.  TPA-induced  Eu(II)  PL  in  Si/Si02.  Figure  4.  PL  of  Eu(III)  dopant  in  Si/SiC>2 

Data  points  are  missing  at  laser  wavelength.  at  25°C.  The  transitions  are  labeled. 
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wavelength.  This  is  inconsistent  with  most  data  obtained  for  Eu(III)  in  crystals.  Actually,  the 
component  in  the  middle  at  —590  nm  is  broadened,  can  be  seen  as  a  twin  peak  due  to  the 
superposition  of  two  subcomponents.  Therefore  the  spectral  features  within  584  -  609  nm  can  be 
assigned  to  the  5Do  -  7Fi  transition  while  the  first  peak  at  580  nm  is  assigned  to  5Do  -  7Fo  .  An 
intense  emission  peaked  at  622  nm  is  due  to  the  5D0  -  7F2  transition.  The  other  two  weak  groups 
of  emission  at  longer  wavelengths  of  659  and  703  nm  belong  to  the  transitions  of  5Do  -  7F3  and 
5Do  -  7F4  respectively.  All  components  were  well  resolved  at  room  temperature. 

CONCLUSIONS 

We  were  able  to  incorporate  Eu(IU)  into  the  Si  nanoparticle  /  Si02  host.  We  also  found 
we  obtained  Eu(II)  in  the  process.  We  found  a  lowering  of  photoluminescence  intensity  with 
lowering  of  temperature.  An  as  yet  unaswered  question  is  the  reason  for  the  intense  whitish 
luminescence  found  in  some  regions  of  the  samples.  Some  involvement  with  Eu(II)  or  Eu(III)  is 
suspected.  Eu(III)  related  peaks  were  only  observed  where  the  peak  distribution  size  of  the 
nanoparticles  was  lower  than  1.3nm.  Whitish  luminescence  was  related  to  peak  sizes  ranging 
from  l.lnm  to  1.4nm.  Annealing  the  samples  had  clear  effects  upon  their  photoluminescence, 
but  did  not  necessarily  involve  changes  in  particle  sizes,  nor  were  these  size  changes  necessary  to 
increase  luminescence.  The  Eu  doping  has  a  tendency  to  halt  the  annealing  effects  on  size  and, 
when  changes  did  occur,  the  particles  generally  became  smaller. 
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